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Introduction 



Ab initio nuclear many-body approaches are based on the premise that the dynamics 
can be modeled studying exactly solvable few-body systems. This is a most important 
feature since, due to the complexity of strong interactions and to the prohibitive difficulties 
associated with the solution of the quantum mechanical many-body problem, theoretical 
calculations of nuclear observables generally involve a number of approximations. Hence, 
models of nuclear dynamics extracted from analyses of the properties of complex nuclei are 
plagued by the systematic uncertainty associated with the use of a specific approximation 
scheme. 

Highly realistic two-nucleon potentials, either purely phenomenological [H El [3l H] or 
based on chiral perturbation theory (ChPT) [5l |6], have been obtained from accurate 
fits of the properties of the bound and scattering states of the two-nucleon system [71 
El m EDI [m 1121 [13 ■ Unfortunately, however, the extension to the case of the three- 
nucleon potential, the inclusion of which is needed to account for the properties of the 
three-nucleon systems, is not straightforward. 

The definition of the potential describing three-nucleon interactions is a central issue 
of nuclear few- and many-body theory. Three-nucleon forces (TNF) are long known to 
provide a sizable contribution to the energies of the ground and low-lying excited states of 
light-nuclei, and play a critical role in determining the equilibrium properties of isospin- 
symmetric nuclear matter (SNM). In addition, their effect is expected to become large, 
or even dominant, in high density neutron matter, the understanding of which is required 
for the theoretical description of compact stars. 

Phenomenological models, such as the Urbana IX (UIX) potential, that reproduce the 
observed binding energy of '^H by construction, fail to explain the measured nd doublet 
scattering length, ^a„rf [14j, as well as the proton analyzing power in p-^He scattering. Ay 

The investigation of uniform nuclear matter may shed light on both the nature and 
the parametrization of the TNF. The equation of State (EoS) of SNM is constrained by 
the available empirical information on saturation density, po, binding energy per nucleon 
at equilibrium, Eq, and compressibility, K. Furthermore, the recent observation of a 
neutron star of about two solar masses [16] puts a constraint on the stiffness of the EoS 
of beta-stable matter, closely related to that of pure neutron matter (PNM). 

Nuclear matter calculations are carried out using a variety of many-body approaches. 
The scheme referred to as Fermi-Hyper- Netted-Chain/Single-Operator-Chain (FHNC/SOC), 
based on correlated basis functions and the cluster expansion technique, has been first 
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used to perform accurate nuclear matter calculations with realistic three body potentials 
in Ref. |T^. This analysis included early versions of both the Urbana (UIV, UV) and 
Tucson Melbourne (TM) three body interactions with the set of parameters reported in 
Ref. [18j. The results indicate that the UV model, the only one featuring a phenomeno- 
logical repulsive term, provides a reasonable nuclear matter saturation density, while the 
UIV and TM potentials fail to predict saturation. In addition, none of the considered 
models yields reasonable values of the SNM binding energy and compressibility. 

The findings of Ref. jT7] are similar to those obtained in Ref. [TP], whose authors 
took into account additional diagrams of the cluster expansion and used the UVII model. 
The state-of-the-art variational calculations discussed in Ref. [20], carried out using the 
Argonne fig ||3| and UIX [21j potentials, also sizably underbinds SNM. 

While the authors of Ref. j2Q] ascribed this underbinding to deficiencies of the vari- 
ational wave function, a signal of the limitations of the UIX three-nucleon potential has 
been provided by the authors of Ref. [2^, who carried out a study of symmetric nuclear 
matter within the Auxiliary Field Diffusion Monte Carlo (AFDMC) approach. Their 
results, obtained using the Argonne v'q NN interaction, show that AFDMC simulations 
do not lead to an increase of the binding energy predicted by Fermi-Hyper-Netted-Chain 
(FHNC) and Brueckner-Hartree-Fock (BHF) calculations |23| . 

Different three-nucleon potential models fllX [21] when applied to the calculation of 
nuclear properties, predict sizably different equations of state (EoS) of pure neutron mat- 
ter at zero temperature and densities exceeding the nuclear matter saturation density, 
Po = 0.16 fm~'^ [25j. In this region, the three-nucleon force contribution to the binding 
energy becomes very large, the ratio between the potential energies associated with two- 
and the three-body interactions being ~ 20% at density p ~ 2po (see, e.g. Ref. [26j). 

In Refs. [27[ [28] . we have followed a strategy which is somewhat along the line of the 
Three-Nucleon-Interaction (TNI) model proposed by Lagaris and Pandharipande [29j and 
Friedman and Pandharipande [30] in the 1980s. 

The authors of Refs. f2U\ I3U] suggested that the main effects of three- and many- 
nucleon forces can be taken into account through an effective, density- dependent two- 
nucleon potential. However, they adopted a purely phenomenological procedure, lacking 
a clearcut interpretation based on the the analysis of many-nucleon interactions at micro- 
scopic level. 

The TNI potential consists of two density- dependent functions involving three free 
parameters, whose values were determined through a fit of the saturation density, binding 
energy per nucleon and compressibility of SNM, obtained from FHNC variational calcu- 
lations. The numerical values of the three model parameters resulting from recent studies 
performed by using AFDMC simulations turn out to be only marginally different from 
those of the original TNI potential |31) . 

The TNI potential has been successfully applied to obtain a variety of nuclear matter 
properties, such as the nucleon momentum distribution [32j, the linear response [55[ IM] . 
and the Green's function [35| [36] . 

The strategy based on the development of two-body density-dependent potentials has 
been later abandoned, because their application to the study of finite nuclei involves a 
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great deal of complication, mainly stemming from the breakdown of translation invariance. 
While in uniform matter the density is constant and the expansion of the effective potential 
in powers of p is straightforward, in nuclei different powers of the density correspond to 
different operators, whose treatment is highly non trivial. 

However, the recent developments in numerical methods for light nuclei seem to indi- 
cate that the above difficulties may turn out to be much less severe then those implied 
in the modeling of explicit many-body forces and, even more, in their use in ab initio 
nuclear calculations. 

The approach described in this Thesis, described in Ref. [27t [28] is based on the 
assumption that n-hodj potentials (n > 3) can be replaced by an effective two-nucleon 
potential, obtained through an average over the degrees of freedom of n — 2 particles. 
Hence, the effective potential can be written as a sum of contributions ordered according 
to powers of density, the p-th order term being associated with {p + 2)-nucleon forces. 

Obviously, such an approach requires that the average be carried out using a formalism 
suitable to account for the full complexity of nuclear dynamics. Our results show that, 
in doing such reduction, of great importance is the proper inclusion of both dynamical 
and statistical NN correlations. Therefore, we have used the Correlated Basis Function 
(CBF) approach and the Fantoni Rosati (FR) cluster expansion formalism to perform 
the calculation of the terms linear in density of the effective potential, arising from the 
irreducible three-nucleon interactions modeled by the UIX potential. 

It should be noticed that our approach significantly improves on the TNI model, as 
the resulting potential is obtained from a realistic microscopic three-nucleon force, which 
provides an accurate description of the properties of light nuclei. 

While being the first step on a long road, the effective potential we have derived 
is valuable in its own right, as it can be used to include the effects of three-nucleon 
interactions in the calculation of the nucleon-nucleon scattering cross section in the nuclear 
medium. The knowledge of this quantity is required to obtain a number of nuclear matter 
properties of astrophysical interest, ranging from the transport coefficients to the neutrino 
emission rates [371 EB] • 

Moreover, the density dependent potential has been implemented in the AFDMC 
computational scheme to obtain the EoS of SNM. Similar calculations using the UIX 
potential are not yet possible, due to the complexities arising from the commutator term. 
Monte Carlo calculations does not show an increase of the binding energy with respect 
to the variational FHNC/SOC calculations. This suggests that UIX potential model may 
not be adequate to properly describe three-nucleon interactions. 

In recent years, the scheme based on ChPT has been extensively employed to ob- 
tain three-nucleon potential models [391 HQ]. The main advantage of this approach is the 
possibility of treating the nucleon-nucleon (NN) potential and the TNF in a more consis- 
tent fashion, as some of the parameters fixed by NN and vrA^ data, are also used in the 
definition of the TNF. In fact, the next-to-next-to-leading-order (NNLO) three-nucleon 
interaction only involves two parameters, namely cd and Ce, that do not appear in the NN 
potential and have to be determined fitting low-energy three-nucleon (NNN) observables. 
Unfortunately, however, ttN and A^A^ data still leave some uncertainties that can not be 
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completely determined by NNN observables. 

A comprehensive comparison between purely phenomenological and chiral inspired 
TNF, which must necessarily involve the analysis of both pure neutron matter and sym- 
metric nuclear matter, is made difficult by the fact that chiral TNF are derived in mo- 
mentum space, while many theoretical formalisms are based on the coordinate space 
representation. 

The local, coordinate space, form of the chiral NNLO three nucleon potential, hereafter 
referred to as NNLOL, can be found in Ref. [41]. However, establishing a connection 
between momentum and coordinate space representations involves some subtleties. 

The authors of Ref. |39] have shown that the NNLO (momentum space) three body 
potential obtained from the chiral Lagrangian, when operating on a antisymmetric wave 
function, gives rise to contributions that are not all independent of one another. To 
obtain a local potential in coordinate space one has to regularize using the momenta 
transferred among the nucleons. This regularization procedure makes all the terms of 
the chiral potential independent, so that, in principle, all of them have to be taken into 
account. The potential would otherwise be somewhat inconsistent, as it becomes apparent 
in nuclear matter calculations, which involve larger momenta. 

Momentum space chiral three-body interaction have been also employed in nuclear 
matter H5| HI]. In these studies, the NNNLO chiral two-body potential has been 
evolved to low momentum interaction Viowk, suitable for use in standard perturbation 
theory in the Fermi gas basis. The results, showing that the TNF is essential to obtain 
saturation and realistic equilibrium properties of SNM H3] . exhibit a sizable cutoff 
dependence. At densities around the saturation point this effect is ~ 4 MeV. In addition, 
different values of the constants Cj lead to different Equations of State for SNM |13] and 

PNM m]. 

A comparative study of different three-nucleon local interactions (Urbana UIX (UIX), 
chiral inspired revision of Tucson-Melbourne (TM') and chiral NNLOL three body po- 
tential), used in conjunction with the local Argonne fig NN potential, has been recently 
performed by the authors of Ref. jlS]. They used the hyperspherical harmonics formal- 
ism to compute the binding energies of and ^He, as well as the nd doublet scattering 
length, and found that the three body potentials do not simultaneously reproduce these 
quantities. Selecting different sets of parameters for each TNF, they were able to obtain 
results compatible with experimental data, although a unique parametrization for each 
potential has not been found. This problem is a consequence of the fact that the three 
low-energy observables considered are not enough to completely fix the set of parameters 
entering the definition of the potentials. 

In Ref. [46] we have analyzed the different parametrizations of the TM' and chiral 
NNLOL three body potential in nuclear matter discussed in Ref. |45] . carrying out nuclear 
matter calculations within both the FHNC/SOC variational scheme and the AFDMC 
formalism. The unphysical cutoff dependence of the contact terms of the NNLOL three 
body potential has been discussed. This feature, that clearly emerges already in the 
illustrative example of noninteracting FG, leads to sizable effects in the more refined 
calculations of the Equation of State of both SNM and PNM. 
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We have found that none of the parametrizations simultaneously reproduces the equi- 
librium properties of nuclear matter. Nevertheless, one of the TM' three-body potentials 
and one of the chiral NNLOL potentials provide values of SNM saturation density close 
to the experimental one. This is a remarkable feature of these potentials, as, unlike the 
UIX model, they do not involve any parameter adjusted to reproduce po- 

Over the past few years, the CBF approach and the cluster expansion formalism have 
been also used to develop well behaved effective interactions, which take into account the 
main effects of NN correlations and are suitable for use in standard perturbation theory 
in the Fermi gas basis, thus allowing for a consistent treatment of equilibrium and non- 
equilibrium properties of nuclear matter [rf, ^STJ. In view of the critical role played by 
interactions involving more than two nucleons, the implementation of the results discussed 
in this Thesis in the CBF effective interaction should be regarded as one of the most 
interesting applications of our analysis. 

As a first step in this direction, we have employed the effective interaction approach 
to compute the weak response of symmetric nuclear matter including three-body clus- 
ter contribution. Although the density, transverse and longitudinal responses of nuclear 
matter at momentum transfer around ~ 1 fm~^ have been already computed within CBF 
using the chain summation scheme EH SB], our approach, based on effective weak 
operators and effective potentials, is more general, as it allows for a consistent description 
of the nuclear response in the regions of both low and high momentum transfer, where 
long- and short-range correlations are known to be dominant. 

In Chapter [H we introduce the concept of nuclear matter, emphasizing its features and 
relations with physical systems. The derivation of two- and three- nucleon interactions 
and their capability of reproducing experimental data are discussed and the formalism of 
chiral perturbation theory is outlined. 

In Chapter [21 after discussing the limitations of the independent particle model, we 
describe CBF theory and the FHNC/SOC summation scheme, as well as the Monte Carlo 
many-body formalisms, pointing out that these approaches are able to encompass the 
correlation structure of nuclear matter, originating from the nuclear interactions. 

The first part of Chapter [3] is devoted to the derivation of the density dependent 
potential, obtained from an average of the UIX three-nucleon force, while in the rest of 
the Chapter we discuss a comparative analysis of the chiral inspired three-nuclear forces 
in nuclear matter. 

Chapter mis focused on the inclusion of the effects of three-nucleon interactions in the 
CBF effective interaction, and the application of this approach to the calculation of the 
weak response of nuclear matter. 



Chapter 1 

Nuclear matter and nuclear interactions 



1.1 Bulk properties of nuclear matter 

Nulcear matter is uniform system of nucleons interacting through strong interactions 
only. While being a theoretical construct, it provides an extremely useful model to inves- 
tigate the properties of both atomic nuclei and neutron star matter. Note that in these 
systems, with the exception of neutron stars in the early stages of their life, the tem- 
perature can be safely set to zero, as thermal energies are negligible compared to Fermi 
energies. 

Two quantities that characterize nuclear matter are the density p and the proton 
fraction Xp 

P = Pp + Pn 

Xp^^ (1.1) 

where Pp and Pn indicate the proton and neutron densities, respectively. PNM is the 
limiting case in which Xp = 0, while for SNM Xp = 1/2. Strong interactions do not bind 
PNM, which in neutron stars is packed by gravitational attraction. SNM on the other 
hand is bound, and its equilibrium properties can be deduced deduced from the analysis 
of nuclear data. 

The nuclear charge distribution, Pcfe(^) is almost constant within the nuclear volume 
and its central value is basically the same for all stable nuclei. It can be parametrized by 

= 1 + Sr-R)/D ■ (1-2) 

Elastic electron-nucleus scattering experiments have shown that the nuclear charge radius, 
R, is proportional to A^^^ 

R = roA-^ , (1.3) 

implying that the volume increases linearly with the mass number. The parameters 
ro = 1.15 fm and D = 0.54 fm have been extracted from experimental data, see for 
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instance [49]. Equation (11 .Sp and tlie nuclear mass formula, to be discussed below, imply 
that the equilibrium 

po = =- 0.16 ±0.02 fm~^ (1.4) 

4:71 

In addition, one can observe that the central charge density of atomic nuclei, measured 
by elastic electron-nucleon scattering, does not depend upon A for large A. As shown in 
Fig- II. 1[ the limiting value does not differ from the one resulting from tq. 



0.20 




Figure 1.1: Saturation of central nuclear densities of medium- heavy nuclei as measured 
by electron- nucleus scattering, from |50] . 



The curves of Fig. II. 1^ parametrized by Eq. (II. 2p . show that the charge density drops 
from 90% to 10% of its value over a distance Rt ~ 2.5 fm, independent on A, called 

surface thickness. 

The (positive) binding energy per nucleon is defined as the difference between the 
mass of the bound nucleus and that of its constituents 



B{Z, A) = -[Znip + iA- Z)mn + Zrue - M{Z, A)] . 



;i.5) 



In Table 11.11 the masses and the binding energies for ^^O, ^^Fe, ^^Ni and ^^°Sn are 
shown. The dependence of the binding energy on the atomic and mass numbers can be 
parametrized according to the semiempirical mass formula [5T| l52] . based on the liquid 
drop model and the shell model 



BiZ, A) = 1 \avA - asA^/' - acZ{Z - \)A^I'' - a^^^^^ + ^p-^^'^' 



1.6) 
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The first term in square brackets, proportional to A is the volume term and describes the 
bulk energy of nuclear matter. It is due to the strong nuclear interaction, that does not 
distinguish between neutrons and protons. Because strong interactions are short ranged, 
a given nucleon may only interact strongly with its nearest neighbors and next nearest 
neighbors, this explaining the scaling with A instead of the A{A — 1) characteristic of 
the long-ranged interaction. The term proportional to A^^^, denoted as surface term is 
also due to the strong interactions. It is actually a correction to the volume term, arising 
from the the fact that nucleons close to the surface have fewer neighbors than the inner 
ones. The third term accounts for the Coulomb repulsion between protons. Since the 
electrostatic interaction is long ranged the scaling is given by Z{Z — 1). The fourth term, 
proportional to [{A — Z) — Z]"^, goes under the name of symmetry energy. Its origin can be 
justified on the basis of the Pauli exclusion principle explaining the experimental evidence 
that stable nuclei tend to have the same number of protons and neutrons. The last term, 
pairing term, which captures the effect of the spin coupling, can be exhaustively explained 
in the framework of the the shell model. It accounts for the fact that even-even nuclei 
(i. e. nuclei having even Z and even A — Z) are likely to be more stable with respect 
to even-odd or odd-odd nuclei. Hence, the value of the constant A is —1, and +1 for 
even-even, even-odd and odd-odd nuclei, respectively. 



Table 1.1: Mass and binding energies of some stable nuclei. 





Z 


M{Z,A) (amu) 


B{Z,A) (MeV) 


160 


8 


15.9949 


7.9765 


56pg 


26 


55.9349 


8.7906 


62Ni 


28 


61.9283 


8.7948 


i20Sn 


50 


119.9022 


8.5048 


208pb 


82 


207.9767 


7.8677 



SNM is described by the semiempirical mass formula by putting Z = A/2 and taking 
the limit for A — >■ oo. Neglecting the Coulomb repulsion, the volume term is the only one 
surviving in this limit. Therefore, the coefficient ay can be identified with the binding 
energy per particle of SNM. Typical fits of Eq. (11. 6p gives |53| 

Eo = E{po) = av = -15.6±0.2MeV (1.7) 

In the vicinity of the equilibrium density, the energy per particle of SNM can be 
expanded according to 

^^ = E, + ^{^^)\0(,-,,r. (1.8) 

The coefficient 
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is the (in) compressibility module that can be extracted from isoscalar breathing modes 
[Ml IS5] , and isotopic differences in charge densities of large nuclei [55] 

fsTo = 220 ± 30 MeV . (1.10) 

1.2 Nuclear hamiltonian 

Although impressive steps have been done in this directions [STj [58| [59| [60] . the de- 
scription of nuclear matter properties at finite density and zero temperature within the 
framework of quantum crhomo-dynamics (QCD) still seems to be out of reach for the 
present computational techniques. 

For this reason, in this work we rely on dynamical models in which non relativistic 
nucleons interacts by means of instantaneous potentials, describing both the short and 
the long range interactions, the latter given by meson (mainly pions) exchange. 

Within this picture, nuclei can be described in terms of point like nucleons of mass m, 
whose dynamics are dictated by the hamiltonian 

i j>i k>j>i 

where Vij and Vijk are the two- and three- body potentials, respectively. In principle four- 
and more- body potential could be included in the hamiltonian, but there are convincing 
indications that their contribution is negligibly small. 

1.2.1 Two-body interaction 

Highly realistic two-nucleon potentials, either purely phenomenological or based on 
chiral perturbation theory (ChPT) have been obtained from accurate fits of the properties 
of the bound and scattering states of the two-nucleon system [3 El El UHl Ell [121 113] ■ 

Besides the technical complexities involved in the fits, from the analysis of nuclear 
experimental data, the following main features of the NN interaction can be inferred. 

• The saturation of nuclear densities, see Fig. II. 1^ showing that the density in the 
inner part of nuclei is almost constant and independent on A, indicates that nucleons 
cannot be packed together too tightly. In a non relativistic picture, a coordinate 
space potential with a repulsive core is able to reproduce this experimental evidence. 
Denoting by ri2 the inter particle distance, one has 

z;(ri2) > , ri2<rc. (1.12) 

It is worth mentioning that the authors of Ref. |42], using the renormalization 
group equations to smear off the repulsive core, were able to develop a class of low- 
momentum potentials. In this framework, the saturation properties is explained in 
terms of many-nucleons interactions. 
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• The nuclear binding energy is almost the same for nuclei with A > 12. Together 
with the proton neutron scattering data, this is an indication for the NN interaction 
to be short-ranged 

v{ru) = , ru>ro. (1.13) 

• Combining the neutron-proton cross section data with the properties of the deuteron 
(^H) it is found that the singlet state does not have bound states, as opposite to 
the triplet state. In particular, the deuteron is the only NN bound state, consisting 
of a neutron and a proton coupled to total spin S = 1 and isospin T = 0. This is a 
clear manifestation of the spin-dependence of the NN. 

• From the observation that the deuteron exhibits a non vanishing quadrupole mo- 
ment, it can be deduced that the angular momentum does not commute with the 
hamiltonian. Consequently the NN potential cannot be invariant under the rotation 
of the spatial coordinates alone. 

• The mirror nuclei are pairs of nuclei such that the proton number in one equals the 
neutron number in the other and vice versa. This obviously implies that mirror nu- 
clei have the same mass number but atomic number differing by one unit. Examples 
of mirror nuclei are ^^N and g'^C or ^^N and g^O. The spectra of mirror nuclei show 
striking similarities, as the energy of the levels with the same parity and angular 
momentum are the same, beside small electromagnetic corrections, showing that 
the nuclear interactions are charge symmetric. This is a manifestation of a more 
general symmetry of the underlying theory, the isospin invariance. 

The first attempt of a theoretical description of NN scattering data is due to Yukawa 
[61j. He made the hypothesis of nucleons interacting through the exchange of a particle 
of mass /i, related with the range of the interaction tq through 

ro~-. (1.14) 

For To ~ l.Ofm the above equation gives fi ~ 200 MeV, which is of the same order of 
magnitude of the pion mass m^r ~ 140 MeV. The most simple parity conserving vertex 
between the nucleons and the pseudo scalar pion has the form ig'j^T, where g is a coupling 
constant and r accounts for the isospin of the nucleons. Hence, the non relativistic limit 
of the scattering amplitude described by the Feynman diagram of Fig. II. 2^ leads to the 
definition of a NN potential, whose expression in coordinate space reads 



2 ^3 



vAri2) = ^7^^n2b(ri2)5i2 + (v;(ri2) - ^5(ri2))aJ , (1.15) 

In the above equation, aij = a,; ■ ffj and Tij = Ti ■ Tj, where (Tj and Tj are Pauli matrices 
acting on the spin or isospin of the i-th, while 

= i^rtf^, - S'-'KcT^ , (1.16) 
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with a, /3 = 1, 2, 3, is tlie tensor operator. Tlie radial functions associated with the spin 
and tensor components read 

Y{x) = —iY{x) (1.17) 

X 

T{x)={l + -^^)Y{x)iT{x). (1.18) 

The phase shift analysis of the high angular momentum neutron-neutron [nn] and neutron- 
proton [np] scattering states shows that for g"^ / (Itt) ^ 14 the one pion exchange potential, 
Vtt, provides an accurate description of the long range part (ri2 > 1.4 fm) of the NN inter- 
action. High angular momentum in fact gives rise to a high centrifugal barrier, preventing 
the nucleons from coming very close to each other. In order to describe NN interactions at 
intermediate range one should consider processes in which more two and more pions, pos- 
sibly interacting among themselves, are exchanged. In the short range region, exchange 
of heavier mesons and more complicated processes, that can be modeled through, e.g., 
contact interactions, are expected to be dominant. 



vr 



Figure 1.2: Feynman diagram describing the one-pion exchange process between two 
nucleons. 



Chiral perturbation theory (ChPT) 

A framework in which the above processes are taken into account in a systematic way 
is provided by the chiral perturbation theory. 

By analyzing the spectra of hadrons having up and down valence quarks, it can be 
deduced that the approximate chiral symmetry of QCD is spontaneously broken. The 
corresponding quasi-Goldstone bosons can be identified with pions, that in turn are much 
lighter than all other hadrons. They would be exactly massless if the masses of u and 
d quarks were vanishing. Goldstone's theorem states that the interactions of Goldstone 
bosons become weak for small momenta, of the order of the pion mass. 

The natural cutoff of ChPT is the rho-meson mass, providing the high-energy scale 
of the theory, ^ SOOMeV. Therefore it is possible to expand the scattering amplitude 
in powers of the small ratios between either the external momenta or the pion mass and 
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A^. Pion loops are naturally incorporated and all corresponding ultraviolet divergences 
can be absorbed at each fixed order in the chiral expansion by counter terms of the most 
general Lagrangian, involving pions and nucleons, consistent with spontaneously broken 
chiral symmetry and other known symmetries. Is it worth noting that the Lagrangian also 
contains contact interaction terms among nucleons, needed to renormalize loop integrals, 
make results fairly independent of the regulators, and parametrize the unresolved short- 
distance dynamics of the nuclear force [62]. 

In the pion-nucleon sector ChPT work well as the interaction vanishes at vanishing 
external momenta in the chiral limit. In the pure nucleonic sector the situation is more 
complicated, since the strong interaction does not become perturbative even in the chiral 
limit at vanishing three-momenta of the external nucleons. In his seminal papers [631 
IM] Weinberg proposed to apply ChPT to the "effective NN potential", defined as the 
sum of connected diagrams for the scattering matrix, generated by old fashioned time- 
ordered perturbation theory. Following this idea Weinberg was able to demonstrate the 
validity of the well-established intuitive hierarchy of the few-nucleon forces: the two- 
nucleon interactions are more important than the three-nucleon ones, which are more 
important than the four-nucleon interactions and so on. 

Within ChPT it also possible to explain why the first order diagram of Fig. I1.2l is able 
to provide an accurate description of the long-range part of the NN potential, although the 
coupling constant g is much larger than one. The one-pion exchange is in fact the leading 
contribution in the chiral parameter |q|/Ai^,, where |q| ~ m^r is the spatial momentum of 
the nucleons. 

Two-body potential at next-to-next-to-next leading order (NNNLO) in the chiral ex- 
pansion has been derived independently by Entem and Machleidt |5J and by the Julich 
group [6]. Both these potentials are able to reproduce the Nijmegen phase shifts — 1- 
Unfortunately, a coordinate space expression for such potentials is not available. There- 
fore, they can be employed neither in FHNC/SOC nor in AFDMC formalism yet. 

On the other hand, a local version of the three-body potential at NNLO does exist 
and will be discussed later in this Chapter 

The Argonne potential 

Phenomenological NN potential [H [21 |3l H] are generally written as 



where '07r(ri2) is given by Eq. fll.lSp stripped of the delta function contribution, Vj{ri2) 
describes the intermediate range attraction attributed to two-pion exchange, and Vsiru) 
accounts for the short range repulsion, which may be due to the exchange of heavier 
mesons and/or to the overlap of the quarks distributions of the nucleons. Comparison 
with ChPT suggests that Vsiru) is strictly related to the contact terms in the chiral 
Lagrangian. The highly realistic Argonne vis potential (AV18) [3j can be written in the 
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form 

18 

^^I8(ri2) = $^t;"(ri2)6?2- (1-20) 

p=i 

The static part of AV18, given by tlie first six operators 

6^=^-^ = (l,a,„S',,)®(l,r,,-). (1.21) 

sufficient to describe deuteron properties and the phase shifts corresponding to S and 
D states. In order to explain the P-wave phase shifts, the spin-orbit term has to be 
introduced 

Of='-' = L,,-S,,®(l,r,,). (1.22) 
In the above equation Lj^ is the relative angular momentum 

= i(r, - r,) X (V, - V,) (1.23) 
and Sij is the total spin of the pair 

S., = i(o-, + o-,). (1.24) 

The remaining 10 operators are required to achieve the description of the Nijmegen scat- 
tering data with — 1. They are given by 

^p=9-i4^(L^LV,,,L-S)®(l,r,,) 

df=''-'' = (1, S.,) ® (T,,, t: + r/) , (1.25) 

where 

= (3r:^4 - '^"Vfr/ . (1.26) 

The last four operators account for the charge symmetry breaking effect, due to the 
different masses and coupling constants of charged and neutral pions. 

Instead of the full AV18, we will be using the so called Argonne fg and Argonne Vq po- 
tentials, which are not simple truncations of the original model, but rather "reprojections" 

m- 

For the purpose of the following discussion, the ^5*1, ^Di, ^Pq and phase shift 
calculated using the Argonne fig potential and its reprojected versions fg and fg are 
displayed in Fig. 11.31 

The Argonne v'g potential is obtained by refitting the scattering data in such a way that 
all 5* and P partial waves as well as the ^Di wave and its coupling to ^5*1 are reproduced 
equally well as in Argonne fig. The differences with the full AVI 8 starts appearing in 
higher partial waves' phase shift, like the plotted of Fig. 11.31 In all light nuclei 
and nuclear matter calculations the results obtained with the fg are very close to those 
obtained with the full fig, and the difference fig — fg can be safely treated perturbatively. 
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The Argonne v'q is not just a truncation of fg, as the radial functions associated with 
the first six operators are adjusted to preserve the deuteron binding energy. Our interest 
in this potential is mostly due to the fact that AFDMC simulations of nuclei and nuclear 
matter can be performed most accurately with Ug-type of two-body interactions. Work 
to include the spin-orbit terms in AFDMC calculations is in progress. On the other 
hand we need to check the accuracy of our proposed density-dependent reduction with 
both FHNC and AFDMC many-body methods before proceeding to the construction of a 
realistic two-body density-dependent model potential and comparing with experimental 
data. 




100 200 300 400 100 200 300 400 

£,„i(MeV) -BioUMeV) 

Figure 1.3: ^^i, ^Di, ^Pq and nucleon-nucleon phase shifts for Argonne Vis (solid 
line), Argonne v'g (dotted line) and Argonne fg (dotted dashed line) potentials, from [66] . 



1.2.2 Three-body interaction 

It is well known that using a nuclear Hamiltonian including only two-nucleon interac- 
tions leads to the underbinding of light nuclei and overestimating the equilibrium density 
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of nuclear matter. Hence, the contribution of three-nucleon interactions must necessarily 
be taken into account. 

In order for the three-body potential to be symmetric under the exchange of particles 
1, 2 and 3 (remember that the sum of Eq. (11. lip has the constraint k > j > i), it has to 
written as a cyclic sum. For all the potentials we are considering in this Thesis, it turns 
out that there are only three independent cyclic permutations 



with V{i : jk) = V{i : kj). 

UIX three-body potential 

One of the most widely used three-body potential is the Urbana IX (UIX) [21], that 
consists of two terms. The attractive two-pion {2tt) exchange interaction V'^'" turns out to 
be helpful in fixing the problem of the underbinding in light nuclei, but makes the nuclear 
matter energy worse. The purely phenomenological repulsive term prevents nuclear 
matter from being overbound at large density. 



Figure 1.4: Feynman Diagram associated with the Fujita Miyazawa three-nucleon poten- 
tial term V^^{3 : 12). 

The V^"" term was first introduced by Fujita and Miyazawa |67| to describe the process 
whereby two pions are exchanged among nucleons and a A resonance is excited in the 
intermediate state, as shown in the Feynman diagram of Fig. 11.41 It can be conveniently 
written as a sum of an anticommutator and a commutator term 



Vi23 = V{1 : 23) + V{2 : 13) + 1^(3 : 12) , 



(1.27) 



A 



71 



f2-(3 : 12) = A2.{Xi3,X23}{ri3,r23} + C2.[Xi3,X23][ri3, rss] 



(1.28) 




(1.29) 



^y{x) = ^t{x) = 1 - e 



— CX' 



.2 



(1.30) 
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In the UIX model, the cutoff parameter is kept fixed at c = 2.1 fm~^, the same value as in 
the cutoff functions appearing in the one-pion exchange term of the Argonne fig two-body 
potential. On the other hand, A2Tr is varied to fit the observed binding energies of '^H. 
The three-nucleon interaction depends on the choice of the NN potential; for example, 
using the Argonne Vis model one gets = —0.0293 MeV. 

The repulsive term is spin-isospin independent and can be written in the simple 
form 

r^(3 : 12) = UoT\m^n3)T^im^r23) , (1.31) 

with T{x) defined in Eq. (ll.lSp . The strength Uq, adjusted to reproduce the empirical 
nuclear matter saturation density, is Uq = 0.0048 MeV with vis- 




Figure 1.5: Energies per particle of ground and low-lying excited states of light nuclei, 
resulting from GFMC calculations, computed with the AV18 and AV18+UIX interactions, 
compared to experiment [24J. The Monte Carlo statistical error are represented by the 
light shaded region. The red dashed lines indicate the breakup thresholds for each model 
or experiment. 



The two parameters ^2,^ and Uq have different values for v'^ and v'^. We disregard 
such small differences in this analysis, mostly aimed at testing the quality of the density- 
dependent reduction of the UIX three-body potential, rather than reproducing empirical 
data. 

As displayed in Fig. II. 5[ showing the results of the GFMC calculations of Ref. [24], 
when the UIX potential is used the binding energy of is exactly reproduced by construc- 
tion, and that of ^He turns out to be very close to the experimental value. A significant 
improvement is also observed for the binding of the p-shell nuclei. However, more and 
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more underbinding is provided by tlie AV18+UIX for increasing A and A — Z. In par- 
ticular a problem with the isospin dependence of the interaction model is revealed by the 
fact that ®He is more underbound than ^Be. With respect to the pure AV18 case, the 
relative stability of the lithium nuclei is improved, but the Borromean helium nuclei are 
still unbound. Additional GFMC calculations of higher-lying excited states, not shown 
in Fig. II. 5[ indicate that the AV18+U1X model underestimate the spin-orbit splittings 
among spin-orbit partners such as the 3/2~ and 1/2^ states in ^He. 

Moreover, the phenomenological UIX model fails to explain the measured nd doublet 
scattering length, ^a„d as well as the proton analyzing power in p-^He scattering. Ay 
[15J. 

Chiral inspired models of three nucleon forces 

In recent years, the scheme based on ChPT has been extensively employed to obtain 
three-nucleon potential models HHl ESI ES]- The main advantage of this approach is 
the possibility of treating the nucleon-nucleon (NN) potential and the TNF in a more 
consistent fashion, as the parameters ci, C3 and C4, fixed by NN and irN data, are also 
used in the definition of the TNF. In fact, the next-to-next-to-leading-order (NNLO) three- 
nucleon interaction only involves two parameters, namely cjj and Ce, that do not appear 
in the NN potential and have to be determined fitting low-energy three-nucleon (NNN) 
observables. Unfortunately, however, ttN and NN data still leave some uncertainties on 
the q's, that can not be completely determined by NNN observables. 

A comprehensive comparison between purely phenomenological and chiral inspired 
TNF, which must necessarily involve the analysis of both pure neutron matter and sym- 
metric nuclear matter, is made difficult by the fact that chiral TNF are derived in mo- 
mentum space, while many theoretical formalisms are based on the coordinate space 
representation. 

The local, coordinate space, form of the chiral NNLO three nucleon potential, hereafter 
referred to as NNLOL, can be found in Ref. [41]. However, establishing a connection 
between momentum and coordinate space representations involves some subtleties. 

The authors of Ref. |3S] have shown that the NNLO (momentum space) three body 
potential obtained from the chiral Lagrangian, when operating on a antisymmetric wave 
function, gives rise to contributions that are not all independent of one another. To 
obtain a local potential in coordinate space one has to regularize using the momenta 
transferred among the nucleons. This regularization procedure makes all the terms of 
the chiral potential independent, so that, in principle, all of them have to be taken into 
account. The potential would otherwise be somewhat inconsistent, as it becomes apparent 
in nuclear matter calculations, which involve larger momenta. 

A comparative study of different three-nucleon local interactions (Urbana UIX (UIX), 
chiral inspired revision of Tucson-Melbourne (TM') and chiral NNLOL three body po- 
tential), used in conjunction with the local Argonne fig NN potential, has been recently 
performed [45J . The authors of Ref. [45J used the hyperspherical harmonics formalism to 
compute the binding energies of '^H and ^He, as well as the nd doublet scattering length. 
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and found that the three body potentials do not simultaneously reproduce these quanti- 
ties. Selecting different sets of parameters for each TNF they were able to obtain results 
compatible with experimental data, although a unique parametrization for each potential 
has not been found. This problem is a consequence of the fact that the three low-energy 
observables considered are not enough to completely fix the set of parameters entering 
the definition of the potentials. 

In a chiral theory without A degrees of freedom, the first nonvanishing three-nuclon 
interactions appear at NNLO in the Weinberg power counting scheme IM] . The 
interaction is described by three different physical mechanisms, corresponding to three 
different topologies of Feynman diagrams, drawn in Fig. 11.61 [39j. The first two diagrams 
correspond to two-pion exchange (TPE) and one-pion exchange (OPE) with the pion 
emitted (or absorbed) by a contact NN interaction. The third diagram represents a 
contact three-nucleon interaction. 




132 132 132 

TPE OPE NNN contact 



Figure 1.6: TPE, OPE and NNN contact interactions of the chiral three-body force at 
NNLO. For the OPE, the diagrams in which 1^2 also needs to be considered. 



As shown in Eq. fll.27p . the full expression for the TNF is obtained by summing all 
possible permutations of the three nucleons. The Feynman diagrams of Fig. 11.61 refer to 
the permutation (3 : 12) of the chiral potential V^, that can be written as 

V^{3 : 12) = CiVi{3 : 12) + C3V3(3 : 12) + CiVi{3 : 12) 

+ CdVd{^ : 12) + CeVe{^ : 12) . (1.32) 

The first three terms Vi, V3 and V4 come from the TPE diagram and are related to 
ttN scattering. In particular, Vi describes the S-wave contribution, while V3 and V4 are 
associated with the P-wave. The other terms, Vd and Ve, are the OPE and contact 
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contributions, respectively. Their momentum space expressions are [39j 

° ^ {qf + ml) {ql + m2 ) 
K) ((Ti ■ qi) (^2 • q2) 

2 (gf + ml) {ql + m2 ) 

K) ^ - \ - / ^ (<?! ■ qi) (0^2 ■ q2) 

— Ta • (ri X T2)a'i ■ (qi x q2)- 



V^i(3 


12) 


V^3(3 


12) 


V^4(3 


12) 




12) 


Ve(3 


12) 



4 (gf + ml) [qi + ml) 

■ q2) ^ , i^i- qi) r 

7^- ^(fTi ■ q2) + ^((T2 ■ qi) 

(g| + m^) [qt + m^) 



K)''n2. (1.33) 
The strengths of the TPE, OPE and contact terms, Vq, V(f and Vq^ are given by 

^"-i%r -'"-s^. -'-jk 

where qa = 1-29 is the axial-vector coupling constant, F^^ = 92.4 MeV is the weak pion 
decay constant and is the chiral symmetry-breaking scale, of the order of the p meson 
mass. 

The low energy constants (LEG) Ci, C3 and C4 also appear in the sub- leading two-pion 
exchange term of the chiral NN potential and are fixed by irN [7D| [7T] and/or A^A^ |^ 
data. The parameters cd and Ce are specific to the three-nucleon interaction and have 
to be fixed using NNN low energy observables, such as the binding energy and the nd 
doublet scattering length "^and |39] . 

The many-body methods employed in our work, namely FHNC/SOC and AFDMC, 
require a local expression of the three-body potential in coordinate space, that can be 
obtained performing the Fourier transform [5T] 



V^{3 ■■m = j ■■ l2)FM)F^{ql)e^^-^^^e^^-^-^ , (1.35) 

where the cutoff functions Fa, defined as 

4 

FA(gf)=exp(-^) , (1.36) 

can depend on the momenta transferred among the nucleons, g^, only. This feature has 
important consequences for the OPE and contact terms, that will be discussed at a later 
stage. 

The cutoff A in the previous equation, while not being required to be the same as 
A^, is of the same order of magnitude. Choosing the fourth power of the momentum in 
Eq. fll.36p is therefore convenient, as the regulator generates powers oi q/A which are 
beyond NNLO in the chiral expansion. 



1.2 Nuclear hamiltonian 



21 



The Fourier transform can be readily computed, and provides the following coordinate- 
space representation of the chiral three-body potential: 

Vi{3 : 12) = Womiai ■ ri3)(a2 ■ 23) y{r 13) y{r 23) 

¥3(3 : 12) = WoTi2[ai2yiri3)y{r23) + (^i ■ r23)(52 ■ r23)^('"23)2/('^i3) 

+ (ai ■ ri3)((T2 ■ ri3)t(ri3)?/(r23) + (ri3 ■ r23)(<?i ■ ri3)(CT2 • r23)t(ri3)t(r23)] 
¥4(3 : 12) = Wo (fg ■ fi X fa) [((73 • (T2 x ai)i/(ri3)?/(r23) 

+ i^3 ■ r23 X ^l)i^2 ■ r23)tir23)yiri3) + (^2 ■ ri3 X 53) ((Ti ■ ri3)t(ri3)?/(r23) 

+ (0^3 ■ r23 X ri3)(ai • ri3)(CT2 • r23)t(ri3)t(r23)] 

Vd(3 : 12) = W^o^Ti2[o-i2y(r23)2:o('"i3) + {^1 ■ r23){^2 ■ r23)t{r23)zo{ri3) 
+ cri2y{ri3)zoir23) + (0^2 ■ ri3)(o'i • ri3)t(ri3)2;o(r23)] 

Ve{S : 12) = iyo''n2^o(ri3)^o(r23) , (1-37) 

where Wq, and are obtained multiplying the corresponding Vq, V^f and V^f by a 
factor m^/(47r)^. The radial functions appearing in the above equations are defined as 




(1.38) 



while Zn, proportional to Z„ introduced in Ref. jT8], is given by 

Zn{.r) = — / X3/ 2 I 2^n ^ = §" / ^ 2 i 2\J o{qr) , (1-39) 

with jo (2;) = sin(a;)/x. Note that, due to the form of the cutoff function of Eq. fll.36p . the 
radial functions are not known in analytic form, and must be obtained from a numerical 
integration. 

Recently, the authors of Ref. |45] have studied the low energy NNN observables using 
the hyperspherical harmonics formalism and a nuclear hamiltonian including the NNLOL 
potential and the Argonne fig [3j two-body interaction. 

This mixed approach requires a fit of all the LECs appearing in the chiral three-body 
interaction, not cd and Ce only. Hence, consistency in the treatment of two- and three- 
nucleon interactions, that would be achievable using a hamiltonian in which all potentials 
are derived from chiral effective theory, is lost. Nevertheless, it is possible to exploit 
chiral perturbation theory to assess the importance of the different terms contributing to 
the TNF. This procedure allows one to select the most relevant spin-isospin structures 
entering the three-nucleon potential, as well as the shape of the corresponding radial 
functions,. 

Within the chiral approach, to obtain a potential yielding a fit to the experimental 
data of accuracy comparable to that achieved by the Argonne fig model, one has to 
include terms up to NNNLO [9, 10]. As a consequence, a fully consistent calculation 



22 



Nuclear matter and nuclear interactions 



in principle requires a NNNLO three-body interaction, the expression of which has been 
only recently derived in Ref. |^ ES]. It turns out that some of the terms appearing at 
NNNLO can be taken into account shifting the constants q of about 20-30% with respect 
to their original values [40j. This procedure has been followed in precision studies of TNF. 
By fitting all the LECs of the NNLOL interaction, the authors of Ref. [45j have improved 
upon the NNLO approximation, as they have effectively included the corrections to the 
Cj appearing at NNNLO level. 

The best fit parameters for the and ^He binding energies and for the nd scattering 
length, "^andi are listed in Table 11.21 For all the different parametrizations, denoted by 
NNLOLi, ci and have been fixed to their original values 0.00081 MeV"^ and 700 MeV, 
respectively [391 . The momentum cutoff of Eq. (11.361) has been set to 500 MeV. 

Table 1.2: Parameters of the NNLOL interactions of Ref. |45] . 



Potential 


C3 (MeV-^) 


C4 (MeV"') 


CD 


ce 


NNLOLi 


-0.00448 


-0.001963 


-0.5 


0.100 


NNLOL2 


-0.00448 


-0.002044 


-1.0 


0.000 


NNLOL3 


-0.00480 


-0.002017 


-1.0 


-0.030 


NNLOL4 


-0.00544 


-0.004860 


-2.0 


-0.500 



As noticed in Ref. |72] . despite the different underlying physical mechanisms, both 
TM and UIX three-nucleon interactions can be written as a sum of terms of the same 
form as those appearing in Eq. f ll.37p . The differences among NNLOL, TM and UIX lie 
in the constants and in the radial functions. 

The TM' potential only involves the Vi, V3 and V4 contributions |73j. The cutoff 
function for this potential is not the same as in Eq. fll.36p . but 

'A2 -m2\2 



;i.40) 



A2 + g2 

The above form allows for the analytical integration of Eq. fll.39p . yielding the radial 
functions 



y[r) 



t{r) 



TV 

„-rA 



m 



Ir"^ - 2(1 + m^r)e"(^-^'"-) + rA(2 + rA) 



— 6 + 2(3 + 3m^r -|- m^r )e 



2^2x r-(A-m^) 



2m^r^ 

+ m^r^(l + rA) - rA[6 + rA(3 + rA)] 



;i.4i) 



The TM' potential corresponds to the following choice of the strength constants (com- 
pare to Eq. fOT]) ) 



Wa 



;i.42) 
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and 

ci = Az , C3 = 26 , C4 = -4rf , (1.43) 

a, b and c being the parameters entering the definition of the TM' potential [73]. The 
authors of Ref.[45j have determined the parameters of the TM' potential fitting the same 
set of low energy NNN observables employed for the NNLOL potential. In order to get 
a better description of the experimental data, they introduced a repulsive three-nucleon 
contact term, similar to the chiral Ve, but with T12 omitted 

Ve{3 : 12) = iyo''^o(ri3)^o(r23) , (1.44) 

where 

W,-=(^Y^-^. (1.45) 
The corresponding radial function can be computed analytically from Eq. fll.39p 

zo{r) = p^iml - Ar . (1.46) 

As in the original paper [TS], in Ref. [45j the value of the pion-nucleon coupling constant 
is set to g"^ = 179.7 MeV, the pion mass is = 139.6 MeV and the nucleon mass is defined 
through the ratio 171^/771.,^ = 6.726. The symmetry breaking scale A^,. of Eq. f ll.45p has 
the same value, 700 MeV, used for the NNLOL potential. 

The parameters of the TM' potentials, TM'^, that according to Ref.[45j reproduce the 
binding energies of and ^He and ^a„d, are listed in Table [TT^ It turns out that Vi, gives 
a very small contribution to the low energy NNN observables. Therefore, the parameter 
a has been kept to its original value — 0.87m~^. 



Table 1.3: Parameters of the TM' potential reproducing low energy the NNN experimental 
data with a = —0.87m~^ 1451. 



Potential 






ce 


A(m^) 


TM[ 


-8.256 


-4.690 


1.0 


4.0 


TM'2 


-3.870 


-3.375 


1.6 


4.8 


TM'3 


-2.064 


-2.279 


2.0 


5.6 



It can be shown that the anticommutator and commutator terms of the UIX potential, 
displayed in Eq. (I1.28p . correspond to V3 and V4 of Eq. f ll.37p . provided the following 
relations between the constants 



C3W0 = 4^2- 



(1.47) 
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and the radial functions 




(1.48) 



are satisfied. 

On the other hand, the repulsive term of the UIX potential of Eq. fll.3ip is equivalent 
to the Ve term appearing in the TM' potential and (aside from the T12 factor) in the 
NNLOL chiral potential if the following relations hold 

T\m^r) = Zo{r) , = CeW^ . (1.49) 

In Ref. [45j it has been found that the original parametrization of the UIX potential 
underestimates ^a^^ and slightly overbinds of ^He. 

The authors of Ref. have calculated the differential cross section and the vec- 
tor and tensor analyzing powers oi p — d scattering at Eiab = 3 MeV for the different 
parametrizations of NNLOL and TM' potentials. They found that all of them lead to 
underestimating Ay (the so-called Ay puzzle remains unsolved) and Tn, while the cen- 
tral minimum in T21 is always overestimated. However, NNLOL model provides a slight 
improvement with respect to the UIX potential in the description of the polarization ob- 
servables. On the other hand, no substantial modifications from the UIX results are given 
by the TM' interactions. 



Chapter 2 



Many body description of nuclear 
matter 



Ab initio nuclear many-body approaches are based on the premise that nuclear dy- 
namics can be modeled studying exactly solvable systems, having mass number A < 3. 
This is a most important feature since, due to the complexity of strong interactions and to 
the prohibitive difficulties associated with the solution of the quantum mechanical many- 
body problem, theoretical calculations of nuclear observables generally involve a number 
of approximations. Hence, models of nuclear dynamics extracted from analyses of the 
properties of complex nuclei are plagued by the systematic uncertainty associated with 
the use of a specific approximation scheme. 

In ab initio approaches, the hamiltonian entering the time-independent many-body 
Schoredinger equation 



is the one defined in Section II. 2^ without any additional adjustable parameters. 

In the first Section of this Chapter we discuss the independent particle model, and 
argue that it is not suitable to encompass correlation structure induced by the nuclear 
hamiltonian. The following Sections are devoted to more advanced approaches, allowing 
one to take into account correlation effects. We will focus on the variational method, based 
on correlated basis function (CBF) theory, and the diffusion Monte Carlo technique. 

2.1 Mean field approach: the Hartree-Fock method 

D. R. Hartree [H], V. A. Fock [75j and J. C. Slater [76j, proposed to use as a start- 
ing point toward the solution of the many-body Schroedinger equation describing atomic 
electrons, the central field approximation. Within this approximation, based on the inde- 
pendent particle model, each nucleon moves in a single-particle effective potential repre- 
senting the average effect of the interactions with the other A — 1 nucleons. Each nucleon 
is described by its own wave function, ipnX^i) eigenfunction of the hermitian operator 




(2.1) 
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h (xi)V'n,(a;i) = e„^Vn,(a;i) , (2.2) 

where the generahzed coordinate Xi = {r^, cr,, Tj} represents both the position and the 
spin- isospin variables of the i-th nucleon. The operator h^^, called the one-particle Fock 
hamiltonian, is given by 

h''^{x,) = -^+v''^{x,), (2.3) 
2m 

where v^^{r) is the single particle Hartree-Fock potential, that is built from the states 
ipmixi) using a self- consistent iterative procedure, based on the variational principle, de- 
scribed in Appendix 1X1 

Within this approximation, the many particle ground-state for a system made of A 
nucleons is a single Slater determinant \1' of one-nucleon states 

^ = A[^n,{x^)...^IJnSXa)]. (2.4) 

where A is the antisymmetrization operator. 

As shown in Appendix 13 the single particle energy is given by 

en. = ^ + / dXj^l)l^{Xi)^^l^{Xj)Vij[ll)nXXi)^n,{Xj) - V^„^. (Xi) (Xj )] , (2.5) 

where J dxj stands for integration over the coordinate and trace over the spin and 
isospin variables of the j-th nucleon. 

The total energy of the system, -E[\l'], is not the sum of the single particle energies, 
but rather 

^[^] = 5Z ^ ^ 5Z / '^^j^n^ i^i)^*n, (^i)^U [^n. {Xi)i^n, {Xj) - ^Pn, (x^^n, (Xj)] . (2.6) 

A physical meaning to the single particle energies can be given through Koopmans' 
theorem. Assuming that the spin orbitals of the A — 1 system are the same as those of the 
A system, from the previous equation it can be shown that e„. is the separation energy of 
the nucleon in the state ipn. 

= Ea - EA-i{t) . (2.7) 

As explained before, the self-consistent field method allows for the determination of 
the spin-orbitals of the A occupied states, {ipi, . . . ,iIja}, with single-particle energies 
ei, . . . , eA, Ea being the Fermi energy of the system. The remaining eigenfunctions of 
h^^, which satisfy Eq. (lA.9p . are associated with unoccupied (virtual) states having 
single particle energies larger than the Fermi energy. Unlike {ipi, . . . ,iPa}, they are not 
determined in a self-consistent fashion, as they do not enter the definition of the Fock 
hamiltonian. 
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The key point of the Hartree Fock approach is that occupied and virtual states provide 
a natural basis to describe the many-body system [77]. While the many-body ground state 
is the Slater determinant of occupied single-particle states, Eq. fl2.4p . excited many-body 
states are constructed by removing n occupied states from the Slater determinant and 
replacing them with n virtual states. Such excited states are called n — particle n — hole 
{np nh) states and are eigenstates of the Hartree Fock hamiltonian, also known as "Fokian" 

i=l 

A ri 

H^^\'^hi,...,h„;pu-,Pn) = [ ^ri, + '^{(^p^ - (^h,) \'^hi,...,h„;pi,...,p„) ■ (2.8) 

ni=l i=l 

The Hartree-Fock procedure is the basis, for instance, of the nuclear shell model, 
that has been successfully applied to explain many nuclear properties. [THl ED]- As 
far as nuclear matter is concerned, the single particle wave functions are known to be 
plane waves, as dictated by translation invariance. Therefore, a uniform system can be 
conveniently described within a box of volume V with periodic boundary conditions |52] . 
using the wave functions 

^nA^i) = —l^Va, , (2.9) 



where r/^. = Xa^Xn represents the product of Pauli spinors describing the spin and the 
isospin of particle i. In order to satisfy the periodic boundary conditions, the wave vector 
k is discretized; for a cubic box of side L, it turns out that 

ki = —ni i = x,y,z rij = 0, ±1, ±2, . . . . (2.10) 

The momentum of the occupied states is smaller than the Fermi momentum hp, which is 
related to the density of the system, p, through kp = (67r^p/z/)^/^, and u is the spin-isospin 
degeneracy (z/ = 2 for PNM, z/ = 4 for SNM). 

The plane waves of Eq. (12. 9p are already solutions to the Hartree Fock equations; 
in other words they are the best single-particle wave functions for uniform systems. A 
remarkable feature of nuclear matter is that the starting single particle wave functions 
are known and simple, unlike what happens, for instance, in finite nuclei. Due to the lack 
of translation invariance, even generating the single particle wave functions is a difficult 
task, as it requires the solution of the Hartree-Fock equations [52]. 

The single particle energy of nuclear matter can be easily derived substituting the 
wave function of Eq. (12. 9p in Eq. (lA.lSp . In the case of SNM (z/ = 4) for potentials of the 



form of Argonne Vis, carrying out the summation over the occupied states with m | < kp 
yields 



2m 



+ pj dr,j [v^^ - l£(A;^r,,)e-^'^-^- (v^^ + 3vl^ + 3v'[j + 9v^;)] . (2.11) 
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where the Slater function is given by 



iikpr) = 3 



(2.12) 



Summing over spin-isospin states of Eq. (I2.5p amounts to tracing over the spin-isospin 
variables of the nucleon j. Such a trace is normalized, as it incorporates the factor l/u 
coming from the summation over the momentum hj. The factor A arising from the same 
sum, divided by wave function normalization factor V produces the factor p, appearing 
in Eq. fl2ll]) . 

Standard perturbation theory performed in the basis of the Hartree-Fock solutions can 
not cope with the repulsive core of the nuclear force, which cause individual terms of the 
perturbative expansion to diverge |20j . 

As an example|81]. consider the scalar repulsive potential 



The single particle energy computed from Eq. fl2.1ip using this potential is seen to be of 
order pro|vo|; if the potential approaches the hard sphere interaction, similar to the strong 
repulsive core of the nuclear interaction, the single particle energy keeps increasing. 

In other words, since the eigenfunctions of the Fock hamiltonian are the same as 
those of the non interacting Fermi gas, the many-body wave function largely differs from 
the exact ground state associated with the nuclear hamiltonian. Standard perturbation 
theory in such a basis can not be expected to be convergent as the matrix elements of the 
nuclear hamiltonian between np nh states are not perturbative corrections to the ground 
state expectation value. 

To circumvent this problem, one can follow two different strategies, leading to either 
G-matrix or correlated basis function (CBF) perturbation theory. 

Within the former approach proposed by Brueckner [82 | [831I8H 185] . the bare potential 
fjj, is replaced by a well behaved effective interaction, the G-matrix, which is obtained 
by summing up the series of particle-particle ladder diagrams. The physical basis of 
this theory was elucidated by Bethe [86j, while Goldstone introduced the linked clus- 
ter expansion [ST]. For a more recent review of the G-matrix approach, also known as 
Brueckner-Bethe-Goldstone expansion, see Refs. |88 [ 

In this Thesis we have been using the CBF approach, to which the following Section 
is devoted. 

2.2 Correlated basis functions theory 

Theories of Fermi liquids based on correlated basis functions are a natural extension 
of variational approaches in which the trial ground state wave function is written in the 
form 






(2.14) 
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where F is suitable many-body correlation operator. The simplest choice suitable for 
dealing with the strong short-range repulsions is the scalar correlator of the form 

A 

^= U fin,) , (2.15) 
j>i=l 

known as Jastrow correlator [81j. However, this choice for the correlation operator is only 
suitable for purely central potentials, such as those describing the interaction between '^He 
atoms. For state-dependent potentials, like the Argonne nuclear interaction, spin-isospin 
dependent correlations, to be introduced at a later stage, are needed. 

The variational approach consists in the minimization of the expectation value of the 
hamiltonian 

Ev = i^o\H\^o) (2.16) 

which is an upper bound to the true ground-state energy Eq. For instance, in the pure 
central Jastrow case, minimizing Ey allows for finding the radial function f{rij). Apart 
from the technical difficulties involved in finding the optimal radial function, it is clear 
that the resulting correlation function is small within the repulsive region of the NN 
potential. 

As noted in the review of Clark ^U\, historically, the development of the variational 
approach has been somewhat discouraged not only by the difficulties involved in the 
calculation of Ey, potentially leading to violations of the variational principle, but also 
by a psychological obstacle: the embarrassing conceptual simplicity of the method, in 
other words, its lack of "snob appeal". 

Nevertheless, the variational approach succeeded in treating the atomic helium in 
both the liquid and solid phases (HH Although nowadays the numerical problem of 
solving the many-body Schrodinger equation for the ground state has been resolved, to 
a large extent, by the Green's function Monte Carlo method [93], this approach does not 
provide a quantitative understanding of the ground-state wave function [9l]. However, the 
knowledge of the analytic form of the ground-state function would be particularly useful to 
extend the microscopic theory to treat the elementary excitations and finite-temperature 
properties of helium liquids. A successful approach in this direction has been provided by 
the variational theory |i95j, including also the back flow correlation, proposed by Feynman 
and Cohen [96]: a velocity dependent correlation, arising from the fiow induced by a 
moving atom. 

As far as the nuclear many body problem in concerned, this method is supported 
by a variety of experimental evidence [97. 98j showing that short range NN correlations 
are a fundamental feature of nuclear structure. The description of nuclear dynamics in 
terms of interactions derived in coordinate space appears to be the most appropriate, 
for both conceptual and technical reasons. First of all, correlations between nucleons 
are predominantly of spatial nature, in analogy with what happens in all known strongly 
correlated systems, like liquid ^He. In addition, one needs to clearly distinguish the effects 
due to the short-range repulsion from those due to relativity. 
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The correlated basis theories of Fermi liquids are a natural extension of the variational 
approach. A non orthogonal but complete set of correlated basis states can be defined as 

= , (2.17) 

(M/„,|^t^|M/„)' 

where {"^n) is the n — particle n — hole state of Eq. ( 12. 8p . The correlation operator, J-", 
is determined by the variational calculation of the ground state energy. The variational 
energies E^, although only has been variationally estimated, are given by the diagonal 
matrix elements of the hamiltonian between correlated states 

El = {mr^\H\^n). (2.18) 

The energies E"^ are extensive quantities, as they are of order A, while excitation energies 
E^ — Eq are of order 1. 

In order to compute the perturbative corrections to the variational energies, the hamil- 
tonian H is decomposed in two terms 

H = H^ + HK (2.19) 

where, as will became clear in the following, neither nor Hi are hermitian operators. 
The "unperturbed" hamiltonian is defined through the correlated basis states and the 
variational energies, in such a way that 

^°|^„) = E^|^„) (2.20) 

Notice that, since the correlated states are not orthogonal, H'^ is not diagonal in this basis 

(^„|i^o|^j = E;;iv„^. (2.21) 

The metric matrix A^„m is defined by 

Nnm = ^nm + Snm = (^n|^m) ; (2.22) 

where Snm is the overlap matrix, with 

5'„„ = 0. (2.23) 

It is convenient to distinguish the diagonal part from the not diagonal part of the 
hamiltonian 

Hnm = (^nl^l^™) = E::^6nm + Km , (2-24) 

where H'^^^ = 0. The closer the CBF states are to the true eigenstates of the hamiltonian, 
the smaller H'^^ becomes. Using Eq. fl2.2ip and (I2.22p it turns out that the off-diagonal 
part of the hamiltonian matrix element is 

^nm ~ E^Snm + -f^nm ) • (2.25) 
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This is consistent with the fact that H^^ = 0, resulting from Eqs. (1238]) . fl219|) and ([222]) • 
Assuming that the nondiagonal elements of both the metric and the hamiltonian be 
small, there have been two fundamental ways of treating the problem pertubatively. One 
way, consists in diagonalizing Hnm as it is, without bothering to orthogonalize the basis, 
using the nonorthogonal perturbation theory. The other way is employing some procedure 
to orthogonalize the basis first, and then apply the standard perturbation theory. 

Within the former approach, the authors of Ref. |1U1] . introducing the so called 
diagonal metric, were able to show that the perturbative corrections to the variational 
energy E"^ can be casted in a way that is formally identical to standard perturbation 
theory in an orthogonal basis 



where En is the exact eigenvalue of the full hamiltonian the eigenstates of which are 
denoted as \^-n\ 

^|^„,} = E„|^4. (2.27) 

The differences with respect to the orthogonal case are enclosed in the matrix V^m, the 
perturbative expansion of which reads 

k 

SnkSkl{H'irn ~ SlmE^) + . . . . (2.28) 

k,l 

Replacing Vnm in Eq- fl2.26p with its expansion leads to 

Z? Z?" \^ i-^'nk~ EnSnk){H'f^n~ EnSkn) /'o On^ 

En-E^ = '^ + . . . (2.29) 

k^n ^« 

Earlier derivations of the latter result, not involving the diagonal metric formalism, 
can be found in Refs. ||102[ 1103] . Like in ordinary many-body perturbation theory, each 
order of the perturbative expansion diverges with the number of particle, A. However in 
Ref. |104j it has been shown that divergent terms appearing at different orders cancel 
each other. 

A major difference with respect to ordinary many-body perturbation theory is that 
there is an energy dependence in the matrix element Vnm of Eq. fl2.26p . arising from the 
non orthogonality of the CBF state. Another peculiar feature of CBF perturbation theory 
is the fact that Vnm is a many-body operator, as, through F, it incorporates the effect of 
the correlations among all the particles of the system. 
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In the earlier calculations |105[I106] . where the correlator was taken to be of the simple 
Jastrow form of Eq. fl2.15p . the second order term of the perturbative corrections has been 
found to be large. The NN potential has indeed a complicate spin-isospin structure, that 
can not be encompassed considering radial correlations only. In particular, since this 
wave function is spherically symmetric, the expectation value of the tensor component 
of the NN interaction averages to zero. In the pure Jastrow case, the CBF states are 
not sufficiently close to the exact eigenstates of the hamiltonian, and more terms in the 
perturbative series need to be calculated. In liquid ^He or in electron gas, where the 
potential is purely central, Jastrow CBF is much more justified [107] . 

A generalization of the Jastrow correlation operator whose structure refiects the com- 
plexity of the NN interaction has been proposed in Ref. |108[ I109[ 1110] 

A 

P{X) = (5 n 4) , (2-30) 

j>i=i 

with 

6 

4 = E/'(^^^)^^ (2.31) 

p=i 

Note that the symmetrization operator S is needed to fulfill the requirement of anti- 
symmetrization of the state |^n), since, in general, [OlpOlf^ 7^ 0. 

Since the first six operators present in the NN potential form a closed set, this choice for 
the correlation operator has a tremendous advantage in analytic manipulation necessary 
to compute the energy per particle. As will be shown in Section 12. 5[ the product of any 
two of the O^, p < 6, can be reduced to a linear combination of elements from this set. 

In this Thesis we will stick to this choice for although in Ref. |111] the correlation 
operator has been extended including spin-orbit correlations. In fact, the variational 
choice of Eq. f l2.3ip (Fg model) implies that spin orbit correlations are neglected. We 
motivate this choice mainly with the technical difficulties of consistently including spin- 
orbit correlations; in spite of the calculations performed of Ref. |111] . we believe that the 
contribution of the spin-orbit correlation is still an open problem. In several FHNC/SOC 
calculations of the binding energy of SNM the spin orbit terms of the potential have been 
included only pertubatively. Moreover in all the FHNC/SOC calculations of the linear 
response |33| IM] . optical potential and Green's function |35[ 136] of SNM, whose results 
have been used to explain a variety of experimental data, spin-orbit correlations have 
been neglected. 

Before moving to the cluster expansion technique, which has been developed to com- 
pute the matrix elements of the hamiltonian, it is worth spending few words on orthogonal 
CBF theory. As a matter of fart, a clear analysis of the convergence properties of the 
non-orthogonal CBF perturbation theory has not been performed yet. For instance, the 
truncation of the series at some perturbative order leads to non-orthogonality spuriousi- 
ties, whose effects may not be negligible. Moreover, the calculation of quantities other 
than the ground-state energy, like the response function, is made difficult by the fact that 
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properly orthogonalized eigenvectors cannot be easily extracted tlie nonortliogonal CBF 
basis. 

If one attempts to diagonalize the CBF states using tlie standard Lowdin transforma- 
tion |112] . the resulting states are worst than the original one. For instance, the expec- 
tation value of the hamiltonian in the Lowdin orthogonal ground-state is larger than E^. 
To avoid this inconvenient, a two-step orthogonalization procedure which preserves the 
variational diagonal matrix elements of the hamiltonian and allows for using ordinary or- 
thogonal perturbation theory in zero temperature calculations, has been developed [ 113] . 

2.3 Cluster expansion formalism 

Both correlation operators of Eq. f l2.15p and fl2.3ip are defined in such a way to 
possess the cluster property. This means that if the system is split in two (or more) 
subset of particles that are moved far away from each other, the correlation operator 
factorizes into a product of two (or more) correlation operators in such a way that only 
particles belonging to the same subset are correlated. For instance, consider two subsets, 
say im and im+i, ■ ■ ■ iA] the cluster property implies 



The above property allows for expanding the matrix elements of the hamiltonian, 
(or of any other many-body operator), between CBF states in sum of terms involving 
an increasing number of particles, knowns as clusters. In the literature are present both 
analytic |90 |ll06j and diagrammatic cluster expansion formalisms |114[I115]I116] : moreover 
different classification schemes have been adopted, corresponding to different choices for 
the smallness parameters of the perturbative expansion. 

In the calculation of the expectation value of any many-body operator it is convenient 
to perform separate cluster expansion for the numerator and the denominator, the latter 
arising from the normalization of CBF states 



It is a general property of the cluster expansion, to be discussed in detail below, that 
divergent terms coming from the expansion of the numerator and of denominator cancel. 

The two following subsections will be devoted to the description of the original Fantoni 
Rosati (FR) diagrammatic cluster expansion formalism |114[ 1115) and to its generalization 
|117[ 1118"] . developed to deal with spin-isospin dependent correlation operators. 

2.4 Fantoni Rosati cluster expansion and FHNC sum- 
mation scheme 

The FR cluster expansion has been obtained through a generalization of the concepts 
underlying the Mayer expansion scheme, originally developed to describe classical liquids 




(2.32) 




(2.33) 
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[119] . to the case of quantum Bose and Fermi systems. 

Within the FR approach, both the term ThT of the numerator Nnm and the term 
TT of the denominator -D„m associated with the expectation value of the hamiltonian 
are expanded in terms of 

h{r,,)^nn,f-l. (2.34) 

Notice that for scalar correlation operator of Eq. f l2.15p to respect the cluster properties 
one can impose 

f^ir >(f) = l. (2.35) 

The variational parameter rf"^, to be fully explained later on, is the central healing distance 
encompassing the fact that when two-particles are further apart than dc they are not 
anymore correlated. Hence, the quantity h{rij) can be seen as a smallness parameter for 
the cluster expansion, as it is indeed in the case of the "power-series" (PS) expansion 
scheme. 



Two-body distribution function 

In the calculation of the ground-state expectation value of any two-body scalar op- 
erator, it is very useful to employ the scalar two-body distribution function, (yf^(ri,r2) 
defined as 

P- }dx, A%(X)TfTMX) 

In terms of g'^{ri, r2), the expectation value of the two-body potential reads 



{v) ^ 5^(0|^;.,|0) = ^^^^{0\vu\0) = ^ [ dr,,,g'^{r,,r,)v{ri2). (2.37) 

i<j 

For the first equality we have exploited the symmetry property of the wave function, 
which is due to the fact that J-" is symmetric and \E'o antisymmetric in the generalized 
particle coordinates. Since nuclear matter is uniform, g'^{ri,r2) = g{ri2), implying that 
(v) diverges with the number of particles. However, the potential energy per particle is 
finite and reads 

^ = f 1 dn^g^nMru) ■ (2.38) 

Particles 1 and 2 are denoted as interacting particles and are distinguished from the other 
particles in the medium. 

The scalar two-body distribution function obeys the following sum rule 

p J rfri2[/(ri2)-l] = -l, (2.39) 



that can be easily derived by integrating Eq. fl2.36p and using translation invariance of 
(yf'^(ri2). Note that the latter results is a consequence of the fact that the scalar two-body 



2.4 Fantoni Rosati cluster expansion and FHNC summation scheme 



35 



distribution function can interpreted as the joint probability of finding two particles with 
coordinates ri and r2. 

Following |100] . in the following subsections we will provide a detailed description of 
the FR cluster expansion of the scalar two-body distribution function. 

Cluster decomposition of T'^T 

In the scalar Jastrow case the product of the correlation operator reduces to 

j-t^=n/'M=n[i+Mr.,)]. (2.40) 

i<j i<j 

It is convenient to put aside the correlation between the interacting particles, denoted 
as "active correlation", while the others are called "passive correlations". Without loss of 
generality we can write 



= f(ri2)(l+ J]x(3)(ri,r2;r,)+ J] X(^)(ri, r^; r„ r,) + . . . ) . (2.41) 

The generic cluster term /^(ri2)X*^"^(rj, . . . , r„) correlates the positions of the two inter- 
acting particles and of the n — 2 medium particles and should be considered as an n-body 
operator. For the sake of clarity we give the explicit expression of the first cluster terms 

X^^\r,, rs; r^) = h{ru) + h{r2^) + h{ru)h{r2i) 
X(^^(ri, T2] Vi, Tj) = h{rij) + h{rii)h{r2j) + h{rii)h{rij) + h{r2i)h{r2j) + h{rii)h{rij) 

+ h{r2i)h{rij) + h{ru)h{r2j)h{rij) + ... (2.42) 

Expansion of the numerator in cluster diagrams 

The cluster expansion of the numerator can be performed by substituting the rhs of 
Eq. dMU) in Eq. ([M 

A(A-r 

num =- 



Tri2 f dx3,...,A/'(ri2)(l+ $^X(3)(ri,r2;r,) + ...)|^o(X)|^ (2.43) 



The integration of the cluster term on the squared modulus of the Fermi-gas wave 
function, which is invariant under the exchange of any particles, gives rise to a combinatory 
factor: 



^ f dx3,...,^X(^nri,r2;r„r„...)|vI/o(X)p 



i<j<,-f^l,2 



{A-2)\ 



{A-N)\{N-2)\ 



J da;3,...,AX(^)(ri, rs; rs, r4, . . . , r^)|^°(X)p . (2.44) 



36 



Many body description of nuclear matter 



Using the above result in Eq. fl2.43p leads to 



num =f{ri2) ^ 



N=2 



{N ~2)\ 



rfr3,...,^X(^)(ri,r2;r3,...,r^)(7i^^^'(ri,...,r^) , (2.45) 



where X^'^^=l. Note that we have introduced the mean-field A^-body correlation function, 
defined as: 



{A-N)\p^ 



Tri,...,jv / rfa;^,...,A|*oWP 



(2.46) 



We now proceed by integrating out the variables x^+i, . . . ,xa from g'jf'^^ by using the 
orthogonality of single particle states. As shown in Appendix [Bl extracting particles 
from the Slater determinant of the ground-state yields 



\E'r 



{A-N)\ 



+---+n.jv+l 



ni<---<n]v 

^[^„,(a;0...^n.(a;^)]*^^"^'-'""i 



Xn+1 ■■■ Xa) 



(2.47) 



where \E'[J*'^'^^''"'"^(xAr+i . . . xa) is the minor ipQ, describing a system of A — particles 
with holes ni <,...,< tiat. The minors satisfy the following orthonormality condition 

dXN+l,...,A'^l'^^''^^""'''' {xn+i ■ ■ ■ XA)^'[^^'''-'''^(xAr+i ... Xa)= S^^i^ . . . (5„^,^ . (2.48) 

With the help of the above equation, the mean-field A^-body distribution function can be 
written in the form 

MF/ \ 

(ri,...,r7v) = 



E Tri,...,Ar A[^i^{Xi) . . . Vl^(x^)]^[?/'„i(Xi) . ..^n^ixN)] 



ni<---<n]v 



P 



N 



E Tri,...,^Ut^(xi)...^t^(x^)^[^.,(xi)...7A.,(x^)]j , (2.49) 



ni,...,njv 



implying that if the number of particles, A^, is larger than the number of quantum states, 
A, the mean-field A^-body distribution function vanishes, i. e. 



g%'''{r,,...,rN) = if N>A^ 



(2.50) 



We preliminary remark that this property is crucial for the exact cancellation of the 
unlinked diagrams of the numerator with the denominator to take place. 
The antisymmetrization operator A can be written in the form 



A — 1 — Pij + E] {P'ljPjk + PikPkj) + • • • • 

i<j<k 



(2.51) 
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For a uniform system like nuclear matter the single particle states are normalized plane 
waves, see Eq. ( I2.9p . Thus, the two-particle exchange operator, defined by the relation 

Piji'nA^i)^n^{Xj) = ^nA^j)^nj{Xi) ■ (2.52) 

can be written as 

4- = X , (2.53) 

where 

^7 = + ^^.^(l + ^^.O = E ^'^f. (2-54) 

p=i 

acts on the spin-isospin degrees of freedom of the nucleons' wave function, while 

P[j = exp[-t{k,-kj)-v,j] (2.55) 

exchanges the radial coordinates of particles i and j. 

Because Pauli matrices are traceless, in the pure Jastrow case, when the traces of Eq. 
f l2.49p are carried out, the exchange operator reduces to its central part 

Pij ^ exp[-z(ki - kj) ■ Tij] . (2.56) 

and one is left with 

g%''^{r,,...,r^)=l^-^j 0*^(rO • . • 0^^(r^)^[0k.(ri) • • • 0r.(k^)]] . (2.57) 



Therefore, g^^ can be written in terms of the Slater function, defined by 

- J2 </'k(r.)0k(r,) . (2.58) 



P 



In the limit of infinite volume, the sum over the discrete momentum can be replaced by 
an integral. Hence it can be easily shown that 

The first terms of the mean-field A^-body distribution function reads 

(n, . . . , r„) = 1 - E + E ^^ir^Mr,k)i{n.) " • • • • (2.60) 

i<j i<j<k 

The factors l/p come from the normalization of the exchange operator of Eq. f l2.56p . In 
particular, producing a two-particle loop, £^(rjj), requires one exchange operator Pij with 
the associated factor l/u. For a loop involving n > 2 particles and n-1 exchange operators. 
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the corresponding factor is (l/z/)'^~^. Moreover, tliere are two possible orderings of tlie 
exchange operators producing loops having more than two particles exchanged, bringing 
and additional factor 2. 

We are now ready to give the general structure of the cluster decomposition for the 
numerator of Eq. (I2.45p . It is very useful to do this pictorially, introducing the so called 
"cluster diagrams". The diagrammatic rules are the following: 

• The diagrams consist of dots (vertices) connected by different kinds of correlation 
lines. Open dots represent the active (or interacting) particles (1 and 2), while black 
dots are associated with passive particles, i.e. those in the medium. Integration over 
the coordinates of a passive particle leads to the appearance of a factor p. 

• The dashed lines, representing the correlations h{vij) and denoted as "correlation 
lines", cannot be superimposed. 

• The statistical factor —l[rij)/v, coming from the expansion of (7^^-^(ri, . . . , r„) is 
represented by an oriented solid "exchange line". The exchange lines must form 
closed loops and, as can be readily seen from the expansion of A in terms of the ex- 
change operators of Eq. ( I2.5ip . different loops cannot have common points. Hence, 
the total exchange pattern consists in one or more non touching exchange loops. 

• Each solid point must be reached by at least one correlation line; in fact in Eq. (I2.45P 
each integration over rj is associated with a term X(^)(ri, T2] rs, . . . , r,, . . . , r^r). 




Figure 2.1: Examples of cluster diagrams. 



Fig. 12.11 shows two examples of cluster diagrams whose analytic expressions read 
dSla) =p j dr^h{n^)h{r23) , (2.61) 

X / dr,dudr,{- ^) ( - ^) ( - ^)Mn3)Mr.5) (2.62) 
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In Fig. 12.21 two examples of forbidden diagrams have been drawn. Diagram (12. 2 l a) is 
not allowed because the solid point 3 is not reached by any correlation line (moreover, 
between points 2 and 4 there are two superimposed dashed lines); on the other hand in 
diagram (12. 21 b) there are two touching exchange loops. 




Figure 2.2: Examples of not-allowed cluster diagrams. 



The cluster terms have no specific prefactor, except those coming from the exchange 
rules and a p factor for each integration. One might wonder where the l/{n — 2)\ of 
Eq. (I2.45P ended up. The factor is due to the counting of the permutations of the n — 2 
internal points and it is automatically taken into account by considering only topologically 
different graphs, or, in other words, by the fact that the labels of the solid points in the 
cluster diagrams are dummy indices. The only remnant of that factor is the inverse of 
what is usually called "symmetry factor", s. This counts the permutations of the solid 
points' labels that, without renaming the integration variables, leave the cluster term 
unchanged. For instance diagrams (a) and (b) of Fig. 12. 3[ being topologically identical, 
give the same contributions 



l6 62 



3«t - - ^4 

\ / 

X 

/ \ 

1 ^2 



(b) 



Figure 2.3: Topologically identical clusters diagrams. 



dOla) = j dr3drih{ri3)h{r3i)h{r2i) , 



(2.63) 
(2.64) 
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As a consequence we take into account only one of them and no prefactors appear. 
On the other hand a prefactor s = 1/2 is associated with the diagram (a) of Fig. f l2.4p . 
because the exchange of points 3 and 4 leads to an identical expression, even without 
relabelling the dummy variables 



(1231a) = 




dr3dr4h{ri3)h{ru)h{r34)h{r23)h{r24) . (2.65) 



The reason of this fact lies in the constraint i < j in the expansion of J^^J-" of Eq. (I2.40p , 
implying that a diagram analogous to diagram (12. 41 a) with the points 3 and 4 exchanged 
does not appear in the cluster expansion. 



3fr - - ^4 

I ^ / I 

X 

I / \ I 



Figure 2.4: Graph with a prefactor s = 1/2 associated with. 



Cluster diagrams may be reducible or irreducible. The integrals corresponding to 
irreducible diagrams cannot be factorized. Obviously an irreducible diagram must be 
linked, i.e. each couple of points must be connected by a sequence of lines, which can 
be both correlation and exchange lines. For instance, diagrams (12.1 l a). (12.31 a). (I2.3[ b). 
(12.41 a) are linked and irreducible, while diagram (12. li b) is unlinked (and reducible). If we 
add to the irreducible diagram (12.11 a) the correlation line /i(r34), we obtain the reducible 
diagram of Fig. (12. 5 p , the expression of which reads 

([231a) = p2 y drsduh{r,s)hir23)h{rs^) , (2.66) 
= P J dr 3h{r is) h{r 23) x p j dr^ihir^i) , (2.67) 

which is just the cluster term of (12.11 a) multiplied by the cluster term of the added part. 

Translational invariance gives rise to a factor V for each unlinked part of the diagram 
but for the one containing the external points 1 and 2. Therefore, the order of magnitude 
of a cluster diagram is V^^~^, where iV„ is the number of unlinked parts. For example, 
the order of magnitude of the linked diagram (12. II a) is 1 while the one of diagram (12. li b), 
having two unconnected parts, is V . 
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3^ - - ^4 

/ \ 
/ \ 

l6 h2 

(a) 

Figure 2.5: Linked and reducible diagram. 

Expansion of the denominator in cluster diagrams. 

The same procedure followed for the expansion of the numerator can be used for the 
denominator. However, in this case there are no interacting particles; hence T'^J^ can be 
conveniently expanded as 

PP= (l + J]x(2)(r„r,)+ ^ X(=^)(r„r„r,) + ...) (2.68) 

i<j i<j<k 

where the cluster term X^^^ correlates N particles. The explicit expressions for N = 2 
and = 3 are 

X(2)(r„r, ) = %,,), 

X'-^\ri,rj,rk) = h{rij)h{rik) + h{rik)h{rjk) + h{'^ij)h{rjk) + h{rij)h{rjk)h{rik) . (2.69) 

Substituting the expansion of J^'^J^ of Eq. fl2.68p in the denominator of Eq. f l2.36p and 
exploiting the invariance of I^E^a^fP under any two-particle exchange, one finds 

^ N f 

den = 1 + 5^ ^ / f/n . . . dvMgf'^iv^, r;v)X(^)(ri, . . . , rjv) , (2.70) 

N=2 ' 

with g^^j^ defined in Eq. flZlel) . 

Comparing Eqs. f l2.45p and fl2.70p . we see that the diagrammatic rules for the denom- 
inator are very similar to those of the numerator. The only difference is that the cluster 
diagrams of the denominator have only solid points, hence their order of magnitude is 
yNu Fig. 12.61 two examples of cluster diagrams coming from the expansion of the 
denominator are depicted. 

Two-body distribution function as a sum of cluster diagrams. 

The ratio of Eq. f l2.36p involves two infinite series of cluster terms, corresponding to 
the expansion of the numerator and of the denominator. 
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Figure 2.6: Examples of cluster diagrams of the denominator. 



Let us consider a generic n body linked (reducible or irreducible) cluster diagram, £„, 
of the numerator [120J . where each internal point is connected to the points 1 and 2 by 
at least one continuous path of correlation and/or exchange lines. Each cluster diagram 
of the numerator can be built as a product of times a factor Uq, with q = A — n, 
representing the sum of all the q body unlinked diagrams 



An example of the above equation is depicted in Fig. 12. 7^ where the diagram fl2.1[ a) 
belongs to £3 and the sum of the diagrams enclosed in round parenthesis is Ua-s- 
Considering the expression of the expansion of the denominator, it is readily seen that 



A 





n=2 
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h2 



Figure 2.7: Series of unlinked diagrams associated with 



Ln = Ma). 



A 




(2.72) 



n=l 



Hence, one might naively think that the denominator cancels the disconnected dia- 
grams of the numerator in the thermodynamic limit, A — 00, only. However, by using 
the property of Eq. (I2.50p it is possible to extend up to infinity the upper limits of all the 
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sums appearing in Eqs. fl2.45p and (12. 70 p . yielding 



oo 



oo 



num 



n=2 q=l 



oo 



den 



(2.73) 



n=l 



that proves the fundamental linked cluster property 



oo 



^2(^12) = ^^n- 



(2.74) 



n=2 



Because of this property, which is a common feature of diagrammatic expansion tech- 
niques, the divergent terms of the numerator and of the denominator cancel out, giving 
rise to finite physical quantities. It has to be remarked that for Fermi systems the exact 
cancellation of unlinked diagrams still holds when spin-isospin dependent correlation are 
considered. 

Perturbative schemes 

In this subsection we briefly present two possible choices for the smallness parameter 
of the cluster expansion, namely, the expansion in the number of points and the power 
series (PS) expansion. 

Expansion in the number of points Within this scheme, the order of magnitude 
of a linked diagram is given by the number of its internal points: expanding the two body 
distribution function at n = m. As a factor p""^ is associated with n internal points, the 
smallness parameter is nothing but the density. 

This was the first scheme adopted for classifying diagrams of Fermi liquids; however the 
series is rapidly convergent for low-density regimes only. This is not the case for nuclear 
matter, as found by the authors of Ref. |121] . They have shown that at the equilibrium 
density, the 1-5 body cluster contributions to the energy of SNM are 22.1 MeV, 43.7 
MeV, 10.8 MeV, 3.4 MeV, and 2.6 MeV, while those with n > 5 give 0.8 MeV. The 
computational cost of this approach exponentially increases with the number of internal 
points: actual calculations do not go beyond the three-body cluster contribution in the 
case of spin-isospin dependent operators. 

It has to be noted that the expectation value of the hamiltonian at any finite order 
of the expansion in the number of internal points is unbound. The main reason for this 
lies in the normalization of the wave function that is not properly taken into account. 
Nevertheless, as it will be fully explained at a later stage, a procedure usually employed 
to obtain the correlation function f"^{r) consists in minimizing the two-body cluster con- 
tribution to the energy per particle imposing constraints on the variational parameters of 
the correlation functions. 
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Power series (PS) expansion Suppose to multiply each correlation line h{rij) by a 
parameter a, it follows that the term of order m in the PS expansion scheme |114] of 
92(1^12), denoted as corresponds to the terms of order a™' resulting from the sum of 
Eq. (I2.74p . Of course the zeroth order corresponds to the mean field results, (73° = 92^^- 

A given order of the power series expansion mixes different orders of the expansion 
in the number of points. For instance the first order in PS includes all of the 2-body 
diagrams, ten 3-body diagrams and nine 4-body diagrams. 

A remarkable feature of the PS scheme is that the sum rule of Eq. (I2.39P for the 
two-body distribution function is satisfied at any order. A very simple argument to proof 
this statement, that can be found in [lOOj, is the following. For any choice of a the sum 
rule reads 

p I dri2($^a'^^7r"(n2) - 1) = -1. (2.75) 

Since the mean-field two-body distribution function satisfies the former sum rule, it turns 
out that 

p I rfri2$^a"^7r"(n2) = 0. (2.76) 

n>l 

The integral of the series is zero for any choice of a within the convergence radius, thus 
the integral of each of the coefficients (yfg"" has to be zero. 

The drawback of the PS expansion is that it does not appropriately describe the short- 
range behavior of the pair function, which is crucial for the calculation of the binding 
energy. The factor /(ri2), while being set aside in Eq. (I2.4ip . does not in fact multiply 
all terms of the expansion. 



2.4.1 Fermi Hyper-Netted Chain (FHNC) and RFHNC 
schemes 

To allow for a good description of both the short and the long range behavior of the 
pair function, the development of a scheme in which infinite sets of cluster diagrams are 
summed up is required. This goal is achieved by the FHNC summation procedure, two 
versions of which had been originally proposed: one by Fantoni and Rosati [114[ I122[ 1115] 
and the other by Krotscheck and Ristig |123[ 1116) . 

The two schemes are essentially complementary: the one by Krotscheck and Ristig is 
better suited to the treatment of long-range correlations, while the FR is more convenient 
for the evaluation of the energy in the presence of strong short-range correlations. Since in 
this work have been dealing with the calculation of the expectation value of the potential, 
we have extensively used the FR approach, that we present for the case of a Fermi liquid 
and a wave function with Jastrow-type correlations (for which FHNC has been originally 
developed). 

In what follows, we will not assume that the reducible clusters diagrams cancel out, 
which is strictly true only for a uniform Fermi liquid described by scalar correlations, 
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as rigorously proved in Ref. [114]. In Section [2.5.31 we will show how the cancellation 
mechanism works in the case of the 3-body cluster contributions. 

The first detailed classification of the cluster diagrams, limited to the bosonic case, can 
be found in the fundamental work of J.M.J, van Leeuwen, J. Groeneveld and J. de Boer 
[124] . while the extension to Fermi systems is extensively analyzed in the more recent 

Refs. fTnnifT2n] . 

Vertex Corrected Irreducible (VIC) diagrams 

In this subsection we will show that, when scalar correlations only are considered, 
reducible diagrams can be included in the calculation of the two-body scalar distribution 
function through vertex corrections to the irreducible diagrams. 

For the sake of introducing the reducible diagrams, the correlation line /i(r34) had 
been attached to the point 3 of diagram f l2.1[ a). It is readily seen that to the point 3 
we can add all the possible linked one-body diagrams, or, in other words, all the linked 
diagrams contributing to the one-body correlation function girs). The net result is that 
the sum of all the reducible diagrams having diagram f|2.1[ a) as the irreducible part, can 
be represented by diagram fl2.1[ a) with the vertex 3 renormalized, in the sense that it is 
vertex corrected by C,d, which in the case of translation invariant system is a constant 

= ^ + (linked one-body cluster terms) = g'^^r) . (2.77) 

A pictorial representation of g'^{r) is given in Figure (12. 8p . 




Figure 2.8: Diagrams contributing to g'^{ri), i.e. to the vertex correction ^rf. The same 
diagrams but the first of the second line contribute to C,e as well. 



Note that there are no restrictions on the kind of cluster terms in the sum, because 
in diagram (12.1 l a) no exchange lines reach the point 3, which is dubbed "d" point (the 
label d comes from "dynamical correlations" which is the name originally given to the 
correlation lines). Conversely, if the reducibility point is reached by exchange lines, like, 
for example, the point 3 of diagram (12.11 b). which is an "e" point, then, to avoid the 
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forbidden superposition of excliange lines, the cluster diagrams of the vertex correction 
needs to be attached to point 3 through a correlation line. It follows that there are two 
kinds of vertex corrections: and for reducibility points of type d and e, respectively. 




Actually there is a third type of vertex correction, occurring when the reducibility 
point is an internal point connected to the rest of the diagram containing the external 
points only through exchange lines. This correction cannot be the full ^e, since any solid 
point must be reached by at least one correlation line, hence in this case the vertex 
correction is = — 1- 

The above procedure can be applied to the external point of diagram f l2.1[ a) and to 
all irreducible diagrams. The conclusion is that the sum of reducible and irreducible 
diagrams may be seen as a sum of vertex corrected irreducible diagrams, which are called 
"VIC diagrams", namely irreducible diagrams whose points carry the vertex corrections 
C,d, and ^c- Note that, taking into account vertex corrections, diagram (12.21 a). without 
the two superimposed /i(r34) lines, is allowed, because point 3 is vertex corrected by ^c- 
Consequently, the second diagrammatic rule does not hold for VIC diagrams that may 
have internal points reached by exchange lines only. 

For a uniform system like nuclear matter, the one body distribution function is equal 
to unity, 5'i(r) = 1. This implies the sum rule = 1 that is an useful check for the 
accuracy of the vertex corrections calculation. 

We define a new set of one-body VIC diagrams Ud{ri) and f/e(ri) in terms of the 
following equations 

[l + f/e(ri)]exp[f/,(ri)] (2.78) 
ee = exp[f/rf(ri)] . (2.79) 

The external point 1 of the diagrams belonging to Ud{ri) is not reached by any exchange 
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lines; conversely in those forming Ue{ri) the point 1 must be reached by a loop of exchange 
lines. Moreover diagrams in both Ud{ri) and f/e(ri) cannot be built by pieces connected 
only by means of the point 1. The exponential in the above equations is due to the fact 
that any number of d-structures forming Ud can be attached to the point 1. Since the 
symmetry factor for n topologically identical structures is 1/n!, the global contribution 
of Ud is given by T.ni^/n\)Udiri)'^. 

For the sake of illustration, some of the diagrams belonging to Ud{ri) and to Ue{ri) 
are depicted in Fig. 12.91 

Simple and composite diagrams 

When in an irreducible diagram it is possible to distinguish two or more pieces that 
are connected with the rest of the diagram by means of the point i and j only, we denote 
these parts as suhdiagrams |124j . Two or more subdiagrams are said to form parallel 
connections between the external points 1 and 2 when the whole diagram consists of two 
or more parts which are only connected by means of points 1 and 2. As no integration 
is carried out over the two external points, a factorization of the integral associated with 
the diagram takes place. 

An irreducible diagram consisting of two or more parallel subdiagrams is called com- 
posite or X — diagram. When such division in parallel subdiagrams is not possible the 
diagram is called simple. For example diagram (12.11 a) is simple, as well as those of Figs. 
12.31 and 12.41 Ignoring the unlinked part, diagram (12.11 b) is composite, as shown in Fig. 

EIDl 

As for the vertex corrections, composite diagrams can be classified on the base of the 
kind of their external points 1 and 2. The set of composite diagrams with only correlation 
line reaching points 1 and 2, like diagram (a) of Fig. 12. IH are denoted with Xdd- When 
two exchange lines are attached to point 1 (2) the corresponding composite diagrams 
are labelled with Xed {Xde)i see diagrams (b) and (c) of Fig. 12.111 Analogously, diagram 
( I2.11[ d) belongs to the set X^e since two statistical lines arrive at both external points. To 
build the ee diagrams with the external points in the same statistical loop, it is convenient 
to define composite diagrams where a statistical open loop starts from the external point 
1 and ends to the external point 2, like for instance diagram (e) of Fig. 12.111 and the 




Figure 2.10: A composite diagram decomposed in simple subdiagrams. 
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Figure 2.11: Classification of composite diagrams based on the kind of the external points. 



one in Fig. I2.1UI Note that the set of these diagrams, denoted by X^c, does not directly 
contribute scalar two-body distribution function. 



Nodal diagrams 

An important concept concerning the diagrams classification is the possible occurrence 
of a "node". Such a node is an internal point through which all possible paths joining the 
external points 1 and 2 pass through. A diagram with one ore more nodes is denoted as 
"nodal" diagram. Diagrams fl2.1[ a). (12. 31 a) and (12. 31 b) are nodal. 

A nodal diagram is also necessarily a simple diagram, but not all simple diagrams are 
nodal, as for example diagram (12. 4p . Non nodal simple diagrams are called "elementary" 
diagrams, oi E — diagrams. Nodal diagrams can be classified according to the same 
scheme, based on the kind of external points, adopted for composite diagrams. Examples 
of diagrams belonging to sets Ndd, Ned, Nde, Nee and Nee are depicted in Fig. 12.121 
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2.4.2 FHNC equations 

Consider a nodal diagram contributing to Nxy{r\2) and label 3 the node closest to 
point 1. All the nodal diagrams can then be found by convoluting the sum of all non- 
nodal 1 — 3 subdiagrams, Xxx'{fiz)) with the set of 3 — 2 subdiagrams Xy/y(r23) + A^y'y(r23) 
(with or without nodes). This leads to the following integral equation 

Nxy{ri2) = J2pf dT3Xx.in3)C,,y,[Xy,y{r23) + Ny,y{r23)] , (2.80) 



x'y' 



the indexes x and y running over the types d and e of the external point. The coefficients 
C,x'y'i accounting for the vertex corrections and for the proper treatment of the exchange 
loops, read 

Crfrf = id , Cde = Ced = ie , Cee = . (2-81) 

In order to better understand the effect of the long-range part of the correlation 
function, it is worth rewriting the convolution Eq. (I2.80p in momentum space 

N{k) = piX{k) [X{k) + N{k)] . (2.82) 

where N{k) and X{k) are the Fourier transforms of the nodal and composite functions, 
respectively. Note that we have omitted the subscripts referring to the kind of verices, 
whose presence is irrelevant to the purpose of this discussion. Solving Eq. fl2.82p for N{k) 
yields 



\ 

\ 

1 6 ^2 








Figure 2.12: Classification of nodal diagrams according to the kind of the external points. 
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As a further simplification, consider the nodal diagram A^"(r) formed by n correlation 
lines h{r). The analytic expression of N^{k) can be easily derived iterating Eq. (I2.82p . 
where X{k) has to be replaced by its first order contribution h{k) 



Therefore, the sum of all the A^"(r) is given by 



(2.84) 



1 - pih{k) ' 



(2.85) 



which is in turn a particular case of Eq. (I2.83p . If h{r) is long ranged, say h{r) ajr^ for 
r — )■ oo, then h{]i) — )■ —a/{2'Kk) for /c — )■ 0. Consequently, each N"'{k) is more divergent 
than h{k) in the long wavelength limit, while their sum is well behaved as it diverges like 
h{k). From this fact, we can gather that the expansion in the number of points diverges 
at any finite order, while the chain summations leads to a well behaved long range limit. 

The iterative equation for Ncc-, cannot be written in a form analogous to that of Eq. 
fl2.80p . It differs more from those given in ref., where because the cyclic nodal diagrams do 
not show not all the cancellations occurring when the cancellation of reducible diagrams 
had been assumed from the very beginning |125] . We need to distinguish two different 
kinds of external points: the point x which is reached by an exchange line and at least 
one correlation line and the point p which is reached by an exchange line only. Hence four 
type of nodal cyclic functions, namely N^^, N^^, N^^ and N^^ and correspondingly three 
type of composite functions, X^f, have to be properly taken into account. For instance, 
the convolution of N^^ with any other of the cyclic functions brings a vertex correction 
^e, whereas the convolution of two pp cyclic functions requires a C,c vertex correction. The 
four nodal equations are given by 

N^:{ru) = P [ dr,^Mri3)[XUr32) + K^i^r,^) + m^^T,^)] 



Xlt{r,2) 



rfr3ee^(ri3)[X,,(r32) + iV-(r32)] + iJ.{r^zWc:{r32) 



iVfi'(ri2) = / rfr3^e^(ri3)iV:,^'(r32) +ec^(ri3)[iV,T(r32) - ^^(r32) 



with 



The total cyclic nodal functions is given by 



NUru) = K^{n2) + N:!{n2) + m:{n2) + iVf^(ri2) . 



(2.86) 



(2.87) 



(2.88) 



2.4 Fantoni Rosati cluster expansion and FHNC summation scheme 



51 



Summing the Eqs. fl2.86p . one can solve the following integral equation for the total Nc, 



Ncciru) = P dr3Xcc{ri3)i 
Jv 

^(ris) 



V 



V 

[rx3 



V 



Xec(r32) + P(r32) 

l{r32 



+ X,,{r32) - V{r32) 



(2.89) 



where V is given by 



^(ri2) = P \ dr^^.X^ris) ^..(raa) + NAr32) - 



m2) 



(2.90) 



At this stage, it is worth introducing the expressions for the partial two-body scalar 
distribution functions 

9M{ri2) = /'(ri2) exp[iV,,(ri2) + Eddi^n^)] , 
9de{'r-i2) = gtdiru) = Nde{ri2) + Xde{ri2) , 
^7ef (ri2) = ^7S,(ri2){iVee(ri2) + E^^ {r^) + [Nde{n2) + Edeiru)?} , 



cexch 
iJee 



(n2) = -ug'id{ri2) NAri2) - -/{ru) + EAri2)\ ' + ElT\n2)g'id{ri2) 



gcciru) = N^ciru) + X^^r^) - -iir^) ■ 



(2.91) 



where Exy{r 12) represent the sum of the xy elementary diagrams. The explanation for 
the presence of the exponential in the above equations is analogous to the one given 
for the vertex corrections of Eq. (I2.79p . Note that, as for the vertex corrections, there 
are no exponentials associated with the e— structures, as two exchange loops cannot be 
superimposed. 

The composite functions, which in turn can be seen as generalized links, are defined 

as 



Xdd{r 12) = gddin2) - Ndd{r 12) - 1 

Xdeiru) = Xediru) = gddin2)[Ndeiru) + Edeiru)] - Ndeiru) , 
Xeeiru) = gddiri2)\Neeiru) + Eeeiru) + [Ndeiru) + Edeiru)? 



V 



Xeeiru) 



Xecirx2) - ^^(ri2) + Eeciru)\ 
Xecirx2) = g'ddiri2)[Neciru) - l^iru) + Eeeiru)] - iVcc(ri2) + ^^ru) 



(2.92) 



The functions Ud,e appearing in eq. f l2.79p and entering the vertex corrections ^^^e are 
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solutions of the following integral equations 



Ud^ p dri2Cd[Xdd{ri2) - Eddir^) - Sdd{r 12) {gdd{r 12) - 1)] 
Jv 



+ ^e[-^de('^12) - Edeiru) - S dd{r 12) 9 de{r 12) 
- Sde{ri2){9dd{r 12) - 1)] + Ed , 

P dfi2Cd[Xed{ri2) - -^6^(^12)] + Ce[^ee(?^12) - E^eiru)] 

Jv 

-^d [S dd{r 12) ged{ri2) + Sed{ri2){9dd{ri2) - 1)] - ^e[See{ri2){9dd{ri2) - 1) 

+Sed{ri2)9de{ri2) + Sdd{r 12) 9ee{r 12) + Sde{r 12) 9ed{r 12) - '2dScc{r 12) 9cc{r 12)] 
-^(ri2)[Ae(ri2) - ^iiru)] + E, , (2.93) 

with 

Sxyiru) = ^N^yiru) + E^y{ri2) , 

A/2=(ri2) = A^cT(n2) + A^cT(n2) = CeP [ dr:,2XUri3)[Xcc{rs2) + N,,{rs2) - -^(ras)] , 

Jv ^ 

^L{n2) = K^ru) + N^Jiru) = N^dru) - A4:(ri2) , (2.94) 

and Ed and Ee are the one-body elementary diagrams with external point d and e re- 
spectively. 

The FHNC equations are solved numerically by means of iterative procedures. At the 
nth step, the nodal diagrams and the 1-body structures resulting from the step n — 1, 
namely Nxy{n — 1) and Ux{n — 1), are employed to compute the partial scalar two-body 
distribution functions as well as the composite functions Xxy{n). Hence, the new nodal 
functions and the 1-body structures are calculated making use of those quantities. 

The iterative procedure is stopped when the difference between the values of a test 
quantity, e.g. the nodal diagrams or the energy, computed in two successive iterations 
is smaller than a given convergence parameter. In order for the procedure to start, the 
nodal diagrams are initially set to zero, while for vertex corrections one sets = = 1- 

For dense systems, like liquid helium or nuclear matter, the convergence may be dif- 
ficult to reach, and one may want to smooth out the iterative process. One of the most 
used technique consists in multiplying by a "mixing parameter", < amix < 1; the nodal 
diagrams resulting from the step n — 1 of the iteration procedure and by 1 — amix those 
obtained in the current step n. Then all the other quantities at the iteration n, like, for 
instance, the composite diagrams, are obtained using the mixture 

amixNxy{n - 1) (1 - a^ix)Nxy{n) . (2.95) 



In order to close the FHNC scheme, an iterative equation for the elementary diagrams 
would be required. However, because of their topological structure, a consistent treatment 
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of the elementary diagrams based on two-body kernel equations, like those for the nodal 
and composite diagrams, is not feasible. 

The simplest approximation consists in neglecting all the elementary diagrams, set- 
ting Exy = (FHNC/0 approximation). Although elementary diagrams are neglected 
in this approach, FHNC/0 provides a very good description of the long-range part of 
the two-body distribution function, and it is also accurate enough at short interparticle 
distances (although Exy{r) are short-range functions). For instance, in nuclear matter 
calculation the elementary diagrams' contribution is likely to be small, provided that 
accurate minimization procedures, like the one described in Section I2.5.4[ are performed. 




Figure 2.13: Some 4-body elementary diagrams. 

A brute force approach based on the direct calculation of the ra— body elementary 
diagrams E^, like the 4— body diagrams of Fig. 12. 13^ and the inclusion of their contribu- 
tions in the FHNC equations, (FHNC/n approximation) does not seem to be a promising 
strategy. In fact, in the cases where the elementary diagrams give appreciable effects, like 
in liquid helium, the FHNC/n series shows poor convergence. 

A more reliable procedure for bosonic systems, denote as scaling approximation (HNC/s) 
[94] amounts to approximating the full set of elementary diagram by a;£'4. The scaling 
parameter, a, is determined by matching the expectation values of the kinetic energy ob- 
tained following the Pandharipande-Bethe and the Jackson-Feenberg prescriptions, pre- 
scriptions, to be discussed in the next Section, which in an exact calculation would give 
the same results. The existence of a scaling property for fermionic systems is not at all 
clear, as there are different classes of elementary diagrams, depending on the kind of the 
external points. Nevertheless, the fermionic generalization of the scaling approximation 
has been used in liquid ^He with some success |126] . 
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Preliminary results obtained using iterative equations with four-body kernels are quite 
encouraging, as the difference with respect to Variational Monte Carlo results turns out 
to be much smaller with respect to the FHNC/0 case |127] . 

The total scalar two-body distribution function is given by the sum of the four different 
kinds of the partial two-body distribution functions, multiplied by the appropriate vertex 
corrections 

9%ru) = Cg^ddiru) + Ueig'^eiru) + 9:,{ri2)] + C'[^/ef ^(^12) + 9:r\n,)] . (2.96) 

The numerical solution of the coupled FHNC integral equations is not trivial at all. 
Moreover, the numerical convergence of the solution does not ensure that this solution 
is acceptable from the physical point of view, as neglecting elementary diagrams could 
in principle lead to a violation of the variational principle. A useful tool to check both 
the numerical and the theoretical accuracy of the calculations is the fulfillment of the 
sum rule of the scalar distribution function of Eq. fl2.39p . Comparing the values of 
the kinetic energy expectation value obtained using the Pandharipande-Bethe and the 
Jackson-Feenberg prescriptions, is also an indicator for the variational principle to be 
respected. 

Using Eqs. fl2.38p and f l2.96p one can evaluate the potential energy contribution to 
the energy per particle. In addition, the knowledge of the two-body distribution function 
allows for the calculation of the kinetic energy per particle, to be discussed in the next 
Section. 



2.4.3 Kinetic energy prescriptions 

In this Section we shall perform a number of manipulations on the kinetic energy 
expectation value. To keep the discussion general, we do not assume that J-" takes the 
simple scalar Jastrow form of Eq. f l2.15p adopted in the FR approach, but only that it 
is a real function of the relative particle coordinates, so that the correlator of Eq. fl2.30p 
satisfies J-"'' = J-". 

The kinetic energy expectation value is given by 

i 

(V?) = den-1 j t/xi,...,^** J-t Vl , (2.97) 

where den= J dxi ,,,^A'^0'^J^'^o accounts for the normalization of the CBF wave function 
In the first line we have exploited the symmetry properties of the correlated wave function, 
as in Eq. (I2.37p . It is worth noting that a crucial role in this context is played by the 
symmetryzation operator S appearing in Eq. (12.30 p . 

The most straightforward form for the kinetic energy is obtained by applying the 
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laplacian to the right 



f2.98) 



This expression of the kinetic energy is called "Pandharipande-Bethe" (PB) |128] . although 
it was first implemented by Iwamoto and Yamada |129| . 

The first term in square brackets generates the Fermi gas energy 



h^e , 3 h^kl 

- A-- 



2m 5 2m 
The full PB kinetic energy is given by 

{T)pB = Tp + W^"^ + Wf + U + Uf. 



(2.99) 



(2.100) 



Integrals involving VfFij are included in W and are completely analogous to those 
arising from the two-body potential. Three body terms with derivatives acting on the 
correlation only, ViFijViFik, are contained in U. On the other hand, Wp and Up accounts 
for the contributions coming from V iFijV iiij and ViFij ■ Viiik, respectively. The different 
contributions can be readily illustrated in the case of pure Jastrow correlations [126] 



W 



kin 



dri2g{ri2) 



Wiring) 



U = p' 



m f'iru) 
C?ri2(iri35(3(ri2,ri3,r23) 



Vi/1ri2) ■ Vi/'^(ri3) 
2m nru)nn,) 



Wf = p drugcciru) 



Uf 



dr^dr 



13 



2m firu) 



ri2 ■ riz 



X 



2m /'=(ri2) 

^{gcc{.ri^)gdy{.ri2)[gdy'ir2^) - 1] + gddiri^)gcc{r2^)gcc{.ri2)} 



(2.101) 



exch 



In order to simplify the notation, cancellation among irreducible diagrams has been as- 
sumed; moreover, in the expression for Up the Abe diagrams appearing in the three-body 
distribution function have been neglected 



5-3(^2, ^^13,^^23) = ^gxx'{ri2)gyy>{r2z)gzz'{. 



(2.102) 



exch 



Note that terms containing Vi\I'o gives zero contribution in direct diagrams after sum- 
mation over kj. The only terms contributing to Up are exchange diagrams in which 
X;i^kiexp(ikij ■ Yij) gives riji{rij). 
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Integrating by parts the last term in square brackets of Eq. fl2.97p and using the 
identity 

Y,[iVi%)T^iV,T)^o - = (2.103) 

i 

yields the "Clark- Westhaus" form of the kinetic energy 

{V^,)cw = num-i j dx^_A%[J'^ J'iV^.'^o) - (V,;^^) ■ (V,;^)^o] 



{T)cw = Tf + W^'^-U, 



In the case of pure scalar Jastrow correlation it turns out that 



(2.104) 



(2.105) 



Integrating once by parts the first line of Eq. fl2.97p leads to — (Vi\E'o'^"'") ■ (ViJ^'^o), 
while yet another integration produces an expression in which the laplacian acts on the 
left, (Vf\E'Qj-'^)J-'\E'o- The Jackson-Feenberg form of the kinetic energy is obtained by 
averaging these contributions 



(Vi)jF = num ^ 



{T)jF = TF + Wr + W^ + U^. 



(2.106) 



The Wb two-body integral contains the kinetic contributions involving derivatives on 
correlations only (V^-Fjj — ViFij ■ ViFij). In the case of central correlations it has the 
form of a two-body integral for a Bose liquid 



W, 



B 



dri2g{ri2) 



/ nri2)vinri2)-r'{n2f 



2m V 



[2.m) 



The W(i, and Ucf, have a fermionic origin as they result from V?\l/o^'o and for a scalar 
Jastrow correlation read 



2m/ 

l^Qdd (?^12 ) -Ecc (?^12 ) (r 12 ) 



{qm - l){[^(ri2) - ^{N,, + E,,)]V?£(ri2) + d'^n^f] 



Am 



[(^7dd(ri2) - l)f (ri2)][(^?..(ri3) - l)f (ria)] 



(2.108) 



Although the three forms of the kinetic energy that we have derived are formally 
equivalent, each has its own distinctive advantages and disadvantages in actual cluster 
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expansion calculations. The CW form has the remarkable feature of not involving second 
derivatives and there are no additional terms arising when one goes from a Bose system 
to a Fermi system. 

Because of the large cancellation among the two-body potential contribution and 
yj/fc*" _l_ Wp, the PB kinetic energy is rather unsensitive to the short-range uncertain- 
ties of gc{r 12)- The three body term of the JF procedure is smaller than the three 
body terms U and Up of the PB and CW prescriptions, making the JF kinetic energy 
essentially unaffected by the approximations involved in the three-body distribution func- 
tion. The drawback of the JF prescription mainly resides in the deficient cancellation 
occurring between Wfp + and the two-body potential contribution. Hence, the JF ki- 
netic energy is more affected by the poor knowledge of the two-body distribution function 
at short distances. 



2.5 Extension to operators: FHNC/SOC 

In this Section we will summarize the extension of the FR cluster expansion scheme, 
extensively used to deal with spin-isospin dependent correlation operators, introduced in 
Eq. flZaOj) 

HX) =[S\[ 4) = (5 n E/'(^^.)Of,) • (2.109) 

j>i=l j>i=l p=l 

The operatorial structure of the correlations and their non commutativity make the 
development of a full FHNC summation scheme for diagrams containing spin-dependent 
correlation prohibitive. Here, we will briefly discuss the so called Single Operator Chain 
(FHNC/SOC) summation scheme, a detailed description of which can be found in |118[ 

In addition to the function h{rij) of Eq. (I2.34p . one has to also consider the products 

2nr,,)r>\r,,) , F''\r,,)r>\r,,). (2.110) 

where the factor 2 of the first quantity accounts for the term in which the central corre- 
lation is on the right of while the operatorial one is on the left and for the reversed 
arrangement. 

For the calculation of the expectation value of the NN potential depending on spin- 
isospin operators, like the AV18 of Eq. fll.20p . it is worth introducing the two-body state 
dependent distribution functions, defined analogously to (7^(ri,r2) of Eq. fl2.36p 



The expectation value of the two-body potential can be conveniently rewritten in terms 
of the two-body state dependent distribution functions 

{v) ^ ^(0|{).,|0) = {0\vu\0) = ^J2 I dv,^,g^{v,,r,)v^{r,2). (2.112) 
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Because of translation invariance, the state dependent distribution functions, like the 
scalar one, depends on the magnitude of the relative distance, (7^(ri,r2) = g^iru)- Thus, 
like for the scalar case, the expectation value of the two-body potential diverges with 
number of particles, while 

^ = §E / drug'iruViru) (2.113) 



p 



is a finite quantity. 

Since the total spin and the total isospin of SNM are both vanishing, the following 
sum rules are satisfied by g^{ri2) 



p J dvi2g''{ri2) = -3 
p j c/ri2fi'^(ri2) = -3 

pjdri2g''^ir,2) = 9. (2.114) 

As far as the expansion of the numerator is concerned, Eq. (I2.4ip can be easily 
generalized 

Y,X^'\xi.X2.Xi)+ X^^\xi,X2;x„xj) + ... . (2.115) 
Note that the cluster terms X" are operators; for example 

X^^\xu X2; X,) = {SFl2FlFli)Ol2{SF^2FuF2i) - F/2 • (2.116) 
The numerator of g^{ri2) can be expanded analogously to Eq. fl2.45p 



,Af-2 
\ ^ P 

num = 

N=2 



p f 

(N -2)1 ^'^^^ / dX3,...,NX^^\xi,X2;X3,...,XN)gN^{xi,...,XN), (2.117) 

N=2 ^ 

where the operatorial A^— body mean field distribution function, defined as 

g^^{x,,...,XM) = p4lv)!;^ / ^^^,...,A(X)*vl>o(X) . (2.118) 
can be readily shown to be (see Eq. (I2.46P and (I2.49P ) 

g^^^{Xi,...,XN) = ^ Y d(a;i) • • •C^(a;7v)^[^/'ni(Xl) . . . Vn^(x7v)] . (2.119) 



2.5 Extension to operators: FHNC/SOC 



59 



P 




Figure 2.14: Operatorial correlation bonds. 



The property of Eq. fl2.5Up holds for , allowing for the cancellation between the 
unlinked diagrams of the numerator and the denominator to take place even in the case 
of operatorial correlations. 

Writing the antisymmetrization operator as in Eq. fl2.5ip and summing over the plane 
wave momenta, for the first terms of we get 



g^ l^Ji, . . . , xatj ^ 



j_(TT 1)2, 

i<j 



+ E ^P^^^!^ + A7^Xn^(r.,)^(r,.)£(rfc,) - . . . ] r/,, . . . r^,, . (2.120) 

i<j<k 

This expression is very similar to the one of Eq. (12.601) . however the sum over the spin- 
isospin states of the latter equation cannot be directly performed, as Tri ... jv? is not 
embodied in the definition of g^^ . 

Substituting back the expression of g^^ in Eq. (I2.117p . one finds 



num 



E 

N=2 



^N-2 



(iV-2)! 



rfr3,...,ArCTri 



X2]X3, 



i<j<k 



(2.121) 



i<j 



The symbol "CTr" denotes the normalized trace of the spin-isospin operators, originating 
from the sum over the spin-isospin states of Eq. fl2.120p and the sum over the spin- 
isospin degrees of freedom, Tri ^r, of Eq. (I2.117p . The factor l/u^ accounts for the 
normalization of the trace, such that CTr(l) = 1. 



Diagrammatic rules 

The diagrammatic rules given in Section 12.41 need to be extended to account for 
2/'^(rjj)/^'*^(rjj) and f^^^{rij)f'^^^{rij). The former is represented by a single wavy line, 
the latter by a double wavy line; in both letter indicating the kind of the operator 

involved in the correlation is placed close to the bond itself, see Fig. 12.141 
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A thick solid line, displayed in Fig. 12.151 has been introduced to represent the inter- 
action term, F12O12F12, of the operatorial two-body distribution function. 

1, p, q 

Figure 2.15: Graphical representation of an operator interaction line. 

Note that the value of the diagrams in general depends on the ordering of the operators, 
operators; hence all permutations need to be considered. 

The diagrammatic classification is analogous to the one used for the FR cluster expan- 
sion technique. As already said, disconnected diagrams of the numerator exactly simplify 
with the denominator and only connected diagrams have to be calculated. 

Unlike the Jastrow case, reducible diagrams do not completely cancel out. At a later 
stage in this Thesis, we will describe the calculation of the two- and three- body cluster 
contributions to the energy per particle. In those simple examples, we will show how to 
deal with reducible diagrams. 



2.5.1 Traces 

As can be realized from Eq. fl2.117p . the calculation of {v'^)/A requires the evaluation 
of the traces of spin-isospin dependent operators present in both the potential and the 
correlations. Since these operators are scalar in the Fock space formed by the product of 
configuration, spin and isospin spaces the Pauli identity can be written as 

{a-Oi){h- (Ji) = a-h + iOi - {a Ab) , (2.122) 

— * 

where a and b are generic vector operators not containing cxj. The latter equation, which 
applies for aj — 7^ also, can be used to express a generic operator product as 

Y[dij = C + iest (2.123) 

where C does not contain any spin-isospin dependent operators while the rest contains 
terms in which each and Tk occurs at most once. Owing to the fact that Pauli matrices 
are traceless 

CTr(afc) = CTr(r,) = 0, (2.124) 

the only contribution of Yl Oij is C. In general C depends on the ordering of the operators 
appearing in n^«i' hence all the possible orderings arising from (iS]^F) needs to be 
properly taken into account. 

The authors of Ref. [118J distinguished three different operatorial structures of the 
cluster term, to the analysis of which we devote the following Sections. 
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Product of operators acting on the same pair 

As a first example, consider a cluster term in which the points i and j are joined by 
two operators, hence 

CTr(6f^.6^^.)=^%., (2.125) 
with = 1,3,3,9,6,18 for p = 1,6. The CTr of diagrams in which more than two 
operators insist on the same pair ij can be easily evaluated with the aid of the 
matrices, defined as 

r 

The values of K^'^^ are given in Table 1 of Ref . [118j . Comparing the last two equations 
it is readily seen that K^i^ = 6pqAp. Using Eq. fl2.126p it turns out that 

CTiiO^p'iPlj) = J2 KP'^'CTiiOtpij) = K^^'A' . (2.127) 

s 

Note that, since operators acting on the same pair of points commute, the order of operator 
in the previous equation is immaterial, thus 

J^pgrj^r ^ j^gpvj^r ^ ^grp^p _ _ _ ^2.128) 



Single operator rings (SOR) 

Single operator rings, like the one showed in Fig l2.16[ are characterized by having at 
most two operators on a given point. The normalized trace of a SOR does not depend on 
the ordering of the operators having the point i in common. Because of the Pauli identity, 
the non commuting terms are indeed linear in either (fi ot fl, thus their trace vanishes. 




Figure 2.16: Four-body SOR diagram. 

Let O^j and Oj^ the only two operators arriving at the point j. Making use of the 
Pauli identity it is possible to completely eliminate the operatorial dependence on point 
j. Integrating over the azimuthal angle (pj and tracing over the spin-isospin degrees of 
freedom of particle j yields 

CTr, J d4>Pld% = J2j dct>,Qldl^ . (2.129) 
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The coefficients ^f^^ depends on the internal angles of the triangle rij, rjk, r-ik 

^ifk = ^<Tr (3 COS^ dj 6tr ^ [-9 COS 9i COS Oj COS 9k 

- 3(cos^ Oi + cos^ Oj + cos^ 0k) + '2] 

prrr _ r 

d?^^^^^^^^ = e!;T, P,g,r = a,t. (2.130) 

The evaluation of SOR diagrams is rather simple: once the operators with one point 
in common are placed next to each other, e. g., O^OJ^O^^ . . . , successive contractions 
over the common points can be made by means of Eq. (I2.129p . Every contraction gives a 
^ factor until at the end one is left with two operators acting on the same pair, resulting 
in a factor AP . 



Multiple operator diagrams 

Consider the normalized trace of the diagram (a) of Fig. I2.17[ where more than two 
operators arrive at both points i and j. In principle, all possible orderings of the operators 
have to be considered in the evaluation of the normalized trace. However, invariance under 
cyclic permutations is a general property of the traces. As a consequence, there are only 
two different orderings of the operators: a "successive" order, in which Of^ and 0\j can be 

placed next to each other, and an "alternate" order, in wich either or O^^^ is placed 
between them. For the successive order, using Eq. fl2.126p . f|2.127p and f|2.129p it turns 
out that 




Figure 2.17: Multiple operator diagram 
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= J2k^''A^' I dct>kC:S''. (2.131) 

On the other hand, for the alternate order one has 

j dct^,CTi{dpidl&;i) = J2 L''- I ■ (2.132) 

To determine the matrix L^'^'" one has to note that either 

CTT{dfp:l,d!^]d%) = (2.133) 

or 

CTr(df^.{d[;,6f^.}q;') = 0. (2.134) 

It can be easily seen that L^'^^ = K^'^'^A'^ and L^'^^ = —K^i^A^ in the former and in the 
latter case, respectively. 

Another possibility that needs to be discussed contemplates two SOR meeting at the 
point i, like in the diagram (b) of Fig. 12.171 Because of the invariance of the trace upon 
cyclic exchanges, again there are only two distinct cases. When the two operators acting 
on the pairs ij and ik are contiguous, it turns out that 

CTT{d0t,Ot,dt^ = 6rr'A^6,s'A^ , (2.135) 



where we have used K^'^^ = 6pqAp. 

In order to deal with the alternate order, we introduce the matrix Drs 

J2 OipiPt, = ^..'^^(1 + Drs)dt, , (2.136) 

where in the case of tensor operators, the above equation implies an integration over the 
azimuthal angle 0. The entries of Drs depend on the kind of the operators and 

g 

4 

Daa = -Drr = -D(<j,r)(f7T,tT) — g • (2.137) 

Thus, for the alternate order trace finds 

GTT{dipiPtPi) = ^rr'A'-il + Drs)Sss'A^ • (2.138) 
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2.5.2 FHNC/SOC approximation 

The technique for summing linked cluster diagrams containing operatorial correlations 
is made technically difficult because of the their non commutativity, which makes a full 
FHNC summation prohibitive. Diagrams having one or more passive operatorial bonds are 
calculated at leading order only. Such an approximation is justified by the observation that 
operatorial correlations are much weaker than the scalar ones. Based on this feature, one 
would be tempted to conclude that the leading order amounts to dressing the interaction 
line with all possible FHNC two-body distribution functions. This is not true as, besides 
the short range behavior, the intermediate range behavior of NN correlations also plays 
an important role that needs to be taken into account. In particular, tensor correlations, 
and to some extent also exchange correlations, have a much longer range than the central 
ones. 

In order to handle this problem, summing the class of chain diagrams turns out to be 
to be of great importance, as remarked in Sec. 12.4.21 for the pure central case (see Eq. 
f l2.83p and the subsequent discussion). 

The above issue is taken care of by summing up the Single Operator Chains (SOC) 
in the corresponding FHNC/SOC approximation |118[ 1130] . SOC are chain diagrams in 
which any single passive bond of the chain has a single operator of the type f'^{rij)f^{rij)0^j 
or —h{rij)i{kFrij) x P^j, with p < 6, or FHNC-dressed versions of them. Note that if a 
single bond of the chain is of the scalar type then the spin trace of the corresponding 
cluster term vanishes, as the Pauli matrices are traceless. Then the SOC is the leading 
order, and at the same time it includes the main features of the long range behavior of 
tensor and exchange correlations. 

The calculation of SOC, as that of FHNC chains, is based upon the convolution integral 
of the functions corresponding to two consecutive bonds. Unlike FHNC chains, however, 
the SOC have operatorial bonds. Therefore, the basic algorithm is the convolution of two 
operatorial correlations having one common point of Eq. fl2.129p . 

The ordering of the operators within an SOC is immaterial, because the commutator 
[Oik, Okj] is linear in ffk and r^, and Pauli matrices are traceless. The only orderings that 
matter are those of passive bonds connected to the interacting points 1 or 2, discussed in 
Eqs. fITTM]) . (Effl), (^rm> and (EM). 

A second important contribution which is included in FHNC/SOC approximation is 
the leading order of the vertex corrections. They sum up the contributions of sets of 
subdiagrams which are joined to the basic diagrammatic structure in a single point, like 
diagram (b) of Fig. 12.171 Therefore, a vertex correction dresses the vertex of all the 
possible reducible subdiagrams joined to it. In the FHNC/SOC approximation they are 
taken into account only at the leading order, i.e. including SOR. Vertex corrections play 
an important role for the fulfillment of the sum rules. 

The full FHNC/SOC equations including the SOR vertex corrections can be found 
in the reference paper |118[ I13U] . For pedagogical purposes, we limit ourselves to the 
equations for SOC diagrams, as this eliminates the problem of the reducible diagrams, as 
all the SOC diagrams are irreducible [120j . 
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We will make the further approximation of neglecting elementary diagrams. 
Regarding the notation, the symbols for nodal and composite diagrams carry an ad- 
ditional index, specifying the operatorial dependence: A^^^, X'^^y. 

The generalization of Eq. (I2.80p accounting for operatorial nodal diagrams reads 

Ky{r^2) = E E CTrsp / dv,XlAri^)eZCx'y'[Xl,y{r,,) + Nl,y{r,^] . (2.139) 

p,q=l x'y' ^ 



Since irreducible diagrams only are present, the factor only selects contributions re- 
specting the Pauli principle 

Cdd = Cde = Ced = I , Cee = 0. (2.140) 

The partial two-body distribution functions, g^y = N^y + X^y, are given by (compare 
to Eq. inm} ) 

g'ddiru) = h^{rumri2) 

+ Nl{r,,) + 2N2,{r,^)Nl{r,,)] , 



(2.141) 



where 



h^{n^) = 2rin2)nn2) + nn2)'Kdiri2) 

C{ri2) = iVcc(ri2) - ^(ri2)/z/ (2.142) 

The composite functions can be effortlessly obtained by subtracting the contribution of 
the nodal diagrams from the partial two-body distribution functions 

Xly{ri2) = gly{ri2) - my{r,2) ■ (2.143) 
The total operator distribution function is given by 

g'iru) = gUru) + 2gUri2) + gLin^) • (2.144) 

In a generic exchange loop all links but one carry an operator dependence. The only 
gap is filled by a dynamical operator to complete the operator chain. Within the SOC 
approximation there are two distinct possibilities: the dynamical operator may be inserted 
to the left (L) or to right (R) of the chain, as in the nodal diagrams (a) and (b) of Fig. 
12. 18^ respectively. To deal with cyclic exchange, we need a new bond 

X'.lA^u) = h%n2Wiru)C{n2) + r(ri2)'iV^^,^(ri2)] - N'^Uru) , (2.145) 
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Figure 2.18: Nodal diagrams with a cyclic exchange loop. 



while the corresponding nodal functions read 

p,q=l 

p,g=l ^ 

Kciri2) = K^iru) + K^in,) . (2.146) 

Finally, the partial two-body distribution functions with circular exchanges is given by 

g:;{n,) = g:,{n,)Ar (2.147) 

The cc nodal functions enter as closed SOR in the generalized equations for X^^ and 
N^^. Moreover, they contribute to the energy expectation value, the full calculation of 
which will not be reported in this Thesis. The interested reader is again referred to the 
Refs. |118[ 1130] . Nevertheless, in the next Section, we do present the calculation of the 
two- and three- body cluster contributions to both the potential and the kinetic energy 
per particle. 



2.5.3 Two- and three- body cluster contribution 

In this Section we analyze in detail the two- and the three- body cluster contributions 
to the energy per particle for the case of a "static" potential, without momentum depen- 
dent terms. This analysis, although being useful for pedagogical reasons, in particular for 
the treatment of the reducible diagrams within the FR summation scheme, will be the 
cornerstone of the effective interaction, that will be developed for the calculation of the 
response. 

With the notation used in Section 12.51 the two-body contribution of T'^Vx^T can be 
expressed as 

Pv^^T = X\xx,x^) = FuVuFu • (2.148) 

2b 
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In what follows, the cluster expansion of the denominator will be disregarded as it is 
understood that, for fermionic systems, the denominator exactly cancels the disconnected 
diagrams of the numerator. 

The two-body contribution of the potential is then 



2b 



[ c?a;i(ia;2</)*^(a;i)(/)*2(x2)A2^^i2A2(l - A2)0nii 

ni,n2 



^"2 



.^2 



(2.149) 



Because of the translation invariance of the system, it is possible to integrate out the 
coordinate of the center of mass R12 = |(ri+r2), so that the two-body cluster contribution 
to the potential energy per particle reads 



A 



^2.150) 



The two-body term of the cluster expansion of the kinetic energy, T, is 



Ki Ki 



2b 



(2.151) 



where the commutator removes the Fermi gas energy, which is a one-body contribution. 
Using the symmetry of the wave functions one gets 



{f)\^^ = Yl J ^a:irfX20:,(Xl)0:,(x2)Fi2[V?,Fi2](l - A2)0ni(Xl)0„2(a^2) • 



rai,n2 



(2.152) 



In order to remove the term with the product of the gradients acting on both the corre- 
lation function F12 and on the plane wave, it is convenient to integrate by part the latter 
expression (see Section r2.4.3p . with the result 



{T)l^ = ^Yl J ^a:irfx2C,(xi)C^(a;2)(ViA2)(ViA2)(l-A2)</)n,(xi 

ni,n2 



1 (Pn2[X2 



[2.153) 



Since V1-F12 = Vi2-Fi2; we can integrate out the coordinate of the center of mass, getting 
the following expression for the kinetic energy per particle 



(T) 
A 



- _^ 
2b 2m' 



p / c/ri2CTri2[(Vi2Fi2)(Vi2Fi2)(l-P/2%)]- 



^2.154) 
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For calculating (Vi-Fi2)(ViFi2) one has to account for the fact that the tensor operator 
depends on fi2, hence 

(ViA2)(ViA2) = J2i(^^f?2)o'i2 + /f2(Vi6?2)][(Vi/f2)o?2 + fU^idl,)] 

= E[(Vi/f2)0?2(Vi/f2)0?2 + f?2i^iOl,)fUV^dt,)] , (2.155) 



where the following property of the gradient of the tensor operator has been used 

Vi/(ri2)Vi^i2 = 0. (2.156) 

The first term of Eq. fl2.155p can be conveniently written in terms of the derivative with 
respect the magnitude of ri2 



(2.157) 



P,9 



Thanks to the relation 



(ViSi2)(Vi5i2) = 4-(6 + 2ai2 + ^12) 

^l2 



(2.158) 



the second term of Eq. fl2.155p turns out to be 



6 



E[/'2(ViO?2)/f2(Vi6?2)] = —ifl2 + fllrune + 2ai2 + ^12) 

'12 



(2.159) 



p,<i 



The three-body cluster contribution appearing in the expansion of J-'^v^J-' is given by 

3b ^-^ 



i>2 



i>2 



(2.160) 



Within the PR diagrammatic scheme, the three-body cluster contribution to ('O12) it is not 
merely the expectation value of the latter result, unlike the two-body case. As a matter of 
fact, the reducible diagrams arising form four-body cluster term of J^^^v^J^, the detailed 
calculations of which can be found in appendix [Cl needs to be taken into account. 

The direct term of the three-body cluster contribution in the PR expansion scheme is 
given by 1121] 



(^^12) 



dir p2 

— / rfri2C?ri3CTri23 



36 



(iSFi2-Fi3F23)-0i2(iSFi2-Fi3F23) 
— -Fl2£'l2-^12(-^l^3 + ^23 ~ 1' 



(2.161^ 
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It includes tlie term 

(5^12^3^23)^)12(5^2^3^23) - ^2^)12^2 (2.162) 

from tlie tliree-body cluster contribution of J-"''fi2J^, whereas the reducible four-body 
diagrams of Fig. 12.191 contribute with the factor 

- FuVi2Fi2{F^s + - 2) . (2.163) 

It is worth remarking that in the pure central Jastrow case, Fij = ffj, the four- and 
three- body reducible diagrams completely cancel, as discussed in Section 12. 4. 1^ and we 
obtain the well-known irreducible contribution 

dir-Jastrow f 

{V12) = ^ / dv^2dv^sr\r^2)vin2){r\ns) - l)(r'(r23) - 1) • (2.164) 




Figure 2.19: Four-body reducible diagrams, v2^_^^i^ and their three-body reduction 



The sum of the diagrams where particles 1 and 2 are exchanged gives 
(^12) 



P12 0^ 

= I rfri2rfri3^'(ri2)CTri23 iSF,2F^sF23)v,2{SF,2F^3F23)P, 

3b Z 



err 
12 



F,2Vl2Fl2{F^3 + ^23-^)Pl 



fJT 

12 



(2.165) 
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Figure 2.20: Four-body reducible diagram, ^45^3^, and its three-body reduction. 



The corresponding four-body diagram producing the term |-Fi2'^^i2-^i2(-^f3 + -^23 ~ 2)-P| 
is drawn in Fig. 12.201 

The diagrams in which particles 1 and 3 are exchanged contributes with 



err 
12 



(^^12) 



Pl3 



3b 



dri2rfri3£2(ri3)CTri23 (5^2^3^23)^)12(5^2^3^23)^'] 



(TT 

13 



Fi2'0l2 -^12 -^1^3 Pl3 



^2.166) 



where the term |-Fi2'£'i2-^i2(-^i3 ~ 1)AT comes from the four-body reducible diagram of 




-> -1 X 




Figure 2.21: Four-body reducible diagram, f|),JL.3b, and its three-body reduction. 



Since the potential is invariant under Xi ^ X2- the diagrams with the exchange 
between particles 2 and 3 give the same contribution reported in Eq. fl2.166p . The 
associated four-body reducible diagram is very similar to the one of Fig. 12.211 but with 
the loop attached to particle 2 instead of particle 1. 

Consider the diagrams with the circular exchange involving particles 1, 2 and 3. In 
this case there are no reducible four-body diagrams that partly cancel the reducible part 
of the three body diagram. In addition, there are no three-body reducible diagrams with 
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circular excliange at all. However, tlie four-body diagram of Fig. \2.21\ witli no correlation 
lines linking particles 1 and 2 to the others, can be reduced to a three-body term, so that 
the three-body diagram with a circular exchange reads 



(^12) 



3b 



p2 / rfri2«(ri2)£(ri3)£(r23)CTri23 (5^2 -^13-^23)'^^12 (5^12-^13^23) P^^P] 



err 
13 



Fl 2 -0 1 2 Pi 2 -^1^3 Pl3 P12 



[2.m) 




Figure 2.22: Four-body diagram, v1^^2,bi ^^^^ contributes to the three-body diagrams 
having a circular exchange between particles 1, 2 and 3. 



As explained in Section [2331 three-body cluster contribution to the PB kinetic energy 
contains terms of the kind Vi(5Pi2Pi3P23). Their explicit expressions can be obtained 
from the corresponding equations for the two-body potential by substituting the first term 
of the normalized traces with 



VuiSFuFiA) ^ -2[5(V?Pi2)Pl3P23] - 2[5(ViPi2) ■ (ViPi3)^23] 



(2.168) 



while 



i)i2Pi2 ^ -2(V?Pi2) (2.169) 
for the second term. Following the notation of Section 12.4.31 terms with (V1P12) are 
denoted by W^^^, those having (V1P12) ■ (V1P13) are included in U. On the other hand, 
the three-body cluster terms belonging to Wp arise from the diagrams where particles 1 
and 2 are exchanged 

P12 



(T) 



■ibWF 



m 



—p' / rfri2rfri3£(ri2)£'(ri2)fi2 ■ CTri23 (5P12P13P23) x 



[5(ViPi2)Pl3i^23)i^] 

and from the ones with circular exchange 



fJT 

12 



A2(ViPi2)(AU ^13-1)^1 



fJT 

12 



(2.170) 



(f) 



3b Wf 



—p' / rfri2rfri3£(ri3)£(r23)f (ri2)fi2 ■ CTri23 {SF^^F.^F^^) 



m 



[5(ViPi2) ^3^23] - A2(ViPi2) ^'3^7^] 



err 
12 



(2.171^ 
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The contributions to U stem from tlie diagrams witli the exchange P13 

Pl3 ff 



(T) 



3b Up m 



p / (iri2(iri3£'(ri3)fi3 • CTri23 (5^12^13^23) x 



and from those having circular exchange 

cir fl^ 
ZbUp 



(f) 



p2 / rfri2rfri3£(ri2)£(ri3)£'(ri3)f23 • CTri23 (5Fi2 -^13-^23) 



[5(ViA2)A3^23]AY^] 



^2.172) 



^2.173) 



Note that in this case there are no subtraction terms arising from reducible diagrams. 



2.5.4 Determination of the correlation functions 

An upperbound to the binding energy per particle, Ey/A, can be obtained by using 
the variational method, which amounts to minimizing the energy expectation value (if) /A 
with respect to the variational parameters included in the model. Its cluster expansion is 
given by 

{H) 



A 



Tp + {AE)2 + higher order terms 



(2.174) 



where Tp is the energy of the non interacting Fermi gas and {AE)2 denotes the contribu- 
tion of two-nucleon clusters 



{AE)2 = Pv^^T + (T) 

26 



2b 



(2.175) 



Neglecting higher order cluster contributions, the functional minimization of {H)/A 
leads to a set of six Euler-Lagrange equations, to be solved with proper constraints that 
force t and fP>^^ to "heal" at one and zero, respectively. That is most efficiently achieved 
through the boundary conditions [221 lll<^ 



Fir > <F) 
dfP{r) , 



5, 



pi 5 



dr 



0. 



^2.176) 



Numerical calculations are generally carried out using only two independent "healing 
distances": dc = d^^^'""^ and df = d^'^. 

Additional and important variational parameters are the quenching factors whose 
introduction simulates modifications of the two-body potentials entering in the Euler- 
Lagrange differential equations arising from the screening induced by the presence of the 
nuclear medium 



6 

E 



(2.177) 
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Figure 2.23: Diagrammatic representation of the two-body cluster contribution {AE)2 of 
Eq. fl2.174p . The thick lines represents both the potential and a kinetic contribution, 
involving derivatives acting only on the correlation functions. The effect of the other 
derivatives is included in Tp. 




The full potential is, of course, used in the energy expectation value. In addition, the 
resulting correlation functions are often rescaled according to 



4 = E/5^'/'(^^^)^'.' (2.178) 



The energy expectation value {H)/A, calculated in full FHNC/SOC approximation is 
minimized with respect to variations of dc, dt, Pp, and ap. 

To determine the best values of the variational parameters we have used a version of 
the "Simulated annealing" algorithm |131] . In metallurgy the annealing procedure consists 
in heating and then slowly cooling a metal, to decrease the defects of its structure. During 
the heating the atoms gain kinetic energy and move away from their initial equilibrium 
positions, passing through states of higher energy. Afterwards, when the metal slowly 
cools, it is possible that the atoms freeze in a different configuration with respect to the 
initial one, corresponding to a lower value of the energy. 

In minimization problems the analog of the position of the atoms are the values of the 
parameters to be optimized, in our case dc, dt, Pp and ap, while the energy of the system 
corresponds to the function that has to be minimized, that in our case is the variational 
energy 

Ev = Ev{dc,dt,Pp,ap). (2.179) 

In the simulated annealing procedure, the parameters dc, dt, Pp, ap are drawn from 
the Boltzmann distribution, exp{—Ev/T), where T is just a parameter of the simulated 
annealing algorithm, having no physical meaning. 

We have used a Metropolis algorithm, with acceptance probability of passing from the 
state s = {dc, dt, (3p, ap} to the proposed state s' = {d'^, d[, (3p, a'p} given by 



Ps,s' = exp 



E{s')-E{ 
f 



^2.180) 



By looking at the distribution of the parameters resulting from the Metropolis random 
walk, it is possible to find the values dc, dt, (3p and ap corresponding to the minimum of 



74 



Many body description of nuclear matter 



Ey, e.g. to the maximum of the Boltzmann distribution. As the fictitious temperature T 
is lowered, the system approaches the equilibrium and the values of the parameters get 
closer and closer to d^, dt, Pp, ap ■ 



2.6 Auxiliary field diffusion Monte Carlo 

A central issue in many-body physics is the evaluation of multidimensional integrals, 
like the one of Eq. fl2.11ip . An alternative to the cluster expansion technique is repre- 
sented by stochastic algorithms using the central limit theorem to compute multidimen- 
sional integrals, known as "Monte Carlo methods" |132[ 1133] . 

Using standard numerical integration methods, like the Simpson rule, the computation 
of a D-dimensional integral requires an exponentially growing number of operations. To 
be definite, in order to estimate the value of a /^-dimensional integral with an accuracy 
e, the quantity of operations scales with e~^. The central limit theorem guarantees that 
Monte Carlo methods scale as e~^, regardless from the dimensionality. 



2.6.1 Variational Monte Carlo 

Variational Monte Carlo (VMC) uses the stochastic integration method for evaluating 
the expectation values for a chosen trial wave function. In the CBF approach, the trial 
wave-function is given by \E'r(X) = (X|J^|$o)- In order to comply with the standard 
notation of Monte Carlo formalism, we do not use the variety of the bra and ket symbols 
introduced for CBF theory and |$o) = l^o)- 

The expectation value of any operator O on the state \E't 



(\E't|0|\E't) 



/ dx,_A'i'UX)O^T{X) 



f ciXi,...,A^?.(X)^T(X) 

can be conveniently written making the spin-isospin sum explicit 



(2.181) 



a/3 



with 



On 



:2.182l 



(2.183) 



V{R)o,p 

The VMC algorithm prescribes to sample the configuration R from the probability 
density 

V{R)ap = ^*Ta{R)'^Tp{R) (2.184) 
and to estimate the integral with the sum 



(2.185) 



{R} 
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where Nc is the number of sampled configurations. 

VMC can be seen as an alternative to the cluster expansion technique, allowing for 
controlling the approximation arising from elementary diagrams and SOC approximation. 
The main drawback of VMC, shared with the CBF, is that the goodness of the result 
entirely depends on the accuracy of the trial wave function. 

In actual facts, neutron matter calculations with for the energy per particle, have 
been limited to 14 nucleons in a box |134j . This is due to the operator structure of the 
correlations, implying a sum over the spin-isospin degrees of freedom of A particles. The 
possible spin states of A nucleons are 2^ and since Z of the A nucleons are protons there 
are A\/Z\{A — Z)\ isospin states. Hence the total number of spin-isospin states is 

2.6.2 Diffusion Monte Carlo 

The diffusion Monte Carlo (DMC) method |135[ I132| I133[ fT36] . overcomes the limita- 
tion of the variational wave-function by using a projection technique to enhance the true 
ground-state component of a starting trial wave function. The trial wave-function can be 
expanded on the complete set of eigenstates of the the full hamiltonian, introduced in Eq. 

(I22Z1) 

l^r) = J]cn|^„). (2.187) 

n 

The evolution in imaginary time r = it/h, that project out the true ground state from 
a trial wave function, provided that it is not orthogonal to the true ground state, i. e. 
co^O 

lim e-^^l^r) = hm V c„e~-^"^|^„) = lim coe~^"^|^o) ■ (2.188) 

T — ^nri T — ^no f ^ t — ^no 



n 



Diffusion Monte Carlo (DMC) is a stochastic projector method for solving the imaginary- 
time many-body Schrodinger equation 

- l-\xij(r)) = {H - Et)\'^{t)) -> |^(r + Ar)) =e-(^-^^)^^|^(r)). (2.189) 
Imposing the initial condition 

|^(r = 0)) = |*t) , (2.190) 

it is readily seen that in terms of the imaginary time wave-functions equation fl2.188p 
implies 

lim |^(r)) = lim coe-(^»-^^)^|^o) • (2.191) 



The energy offset is adjusted to be as close as possible to Eq with the aim of making 
the damping exponential factor constant. 
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For the sake of simplicity, at this point we limit ourselves in considering the 3N spatial 
coordinates only. An entire Section will be dedicated to the inclusion of spin-isospin de- 
grees of freedom through the auxiliary fields. Inserting a completeness on the orthonormal 
basis {\R')} in Eq. (127789]) yields 

|^(r + At)) = J dR'e-^^-^^^^^\R')'^{R', r) . (2.192) 

Projecting on the coordinates {R'\ leads to 

^(i?,r + Ar) = j dR!G{R,R!,AT)^{R!,T) (2.193) 

where the kernel 

G{R, R', At) = (i?|e-(^-^^)^ (2.194) 

is Green's function of the operator H + Notice that no assumptions on the smallness 
of At have been made so far. In terms of the wave function the Schrodinger Eq. fl2.189p 
reads 

- ^^(^, r) = {H- Et)^{R, t) . (2.195) 
If we neglect the interaction terms in the Hamiltonian 

# = f = -|^5^V^ (2.196) 

i 

the imaginary time Schrodinger equation fl2.195p reduces to the master equation of a 
diffusion stochastic process [137] 

The associated Green's function is a 3N-dimensional Gaussian for the spatial coordi- 
nate having variance r in each dimension 



3A 



describing the Brownian diffusion of N particles with a dynamic governed by random 
collisions. 

No DMC algorithm is necessary to solve this problem, as the limit r — > oo leading to 
an uniform wave function can be analytically taken. However, to explain how the DMC 
algorithm works, let us represent represent the distribution \E'(i?, r) by a set of discrete 
Brownian sampling points or random walkers 

^{R,T) = Y,SiR-Rk). (2.199) 
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Evolving this discrete distribution for an imaginary time Ar by means of Eq. fl2.193p we 
obtain a set of Gaussians, centered in the positions Rk 

^{R,T + AT) = Y,Gd{R,Rk,AT) (2.200) 

k 

In order to get a discrete representation for the positions of the walker at r + Ar, each 
Gaussian is sampled by a new delta function. The procedure of propagation/resampling 
is then iterated until convergence is reached. 

When the full Hamiltonian, including both the potential and the kinetic term, is 
considered 

H = f + V, (2.201) 

the analytic solution of the associated Green's function {R\e~^'^~^^~^^^^^\R') is in general 
not known. An approximation to the Green's function can be obtained using the Trotter- 
Suzuki formula 

g(i+B)Ar _ g-iAr/2g-BArg-iAr/2 ^ (^(^^3) (2.202) 

In the limit of small Ar the Green's function can be factorized 

(i?|e-(^+^"^^)^^|i?') = Gb{R, R', AT)Gd{R, R', Ar) . (2.203) 
with the branching factor given by 

f V(R) + V(R') TP \ A 

Gb{R, R', Ar) = e~i ^ ""^i^" . (2.204) 

The wave function at r + Ar can then reads 

/ m \ ^ C , m(R~R')^ V(R) + V{r') a CI a 

^(i?, r + Ar) = ( ) / rfi?'e-^^e-^^^^"e^^^"^(i?', r) . (2.205) 

\2-nn^l\TJ J 

The long imaginary time evolution is obtained by iterating the last equation, which is 
valid for Ar — )• 0, until convergence is reached. Note that since the error in Eq. fl2.202p 
is (9(Ar'^), iterating for a finite amount of imaginary time leads to an overall error of 
the order Ar^. On the other hand, using the less-refined version of the Trotter-Suzuki 
formula 

g(i+B)Ar ^ g-iArg-BAr ^ ^(^^2) _ (2.206) 

provides a still acceptable bias of the order Ar. 

The integral of Eq. ( I2.205P is solved using the Monte Carlo diffusion algorithm, the 
main step of which can be summarized as follows |138] : 

• The initial set of walkers is sampled from the distribution \I'(i?',r = 0) = '^t{R') 
and the starting trial energy Et is chosen, for instance from a variational calculation. 

• The coordinates of the walkers are diffused by means of a Brownian motion 

R = R! + i (2.207) 

where ^s a stochastic variable distributed according to a Gaussian probability 
density with a = hAr/m and zero average so that the walkers are distributed 
according to Gd{R, R', Ar). 
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• The branching or birth/death algorithm is applied: the weight 

w = GbiR, R' , At) (2.208) 

is assigned to each walker and a number of copies of the walker proportional to w 
is generated. Then the convolution theorem governing the composition of random 
variables guarantees that the distribution of the walker is "^{R, t + At). In actual 
facts, the integer number of copies is given by the branching factor 

m = INT(w + ri) , (2.209) 

where INT denotes the integer part of a real number and r/ is a random number 
drawn from the uniform distribution on the interval [0, 1]. The energy offset Et is 
adjusted to keep the the total population of walkers fluctuating around a desired 
value. 

• Once convergence is reached, i.e. for large enough r, the conflgurations are dis- 
tributed with a probability density \E'(-R, r). Therefore, the ground-state expectation 
values of observables that commute with the hamiltonian 

^^^^0^^^,. {^Tmir)) JdR{^T\d\R)^{R,T) 

^ ' (^ol^o) (^Tl^(r)) r^oo j dR^T{R)^{R,r) ^ • ' 

can be computed by 

2.6.3 Importance sampling 

The basic version of the DMC algorithm described in the previous Section is poorly 
efficient, as the brownian diffusive process ignores the shape of the potential. Hence, 
the weight of Eq. fl2.208p suffers of large ffuctuations from step to step, as, for instance, 
there is nothing that prevents two-particles from moving very close to each other even in 
presence of an hard-core repulsive potential. 

The idea of the importance sampling technique consists in using the knowledge of the 
trial wave function '^t{,R) jl35[ I139 j to guide the diffusive process. 

Let us multiply the imaginary time Schrodinger equation fl2.195p by the trial wave 
function ^t{R) and introduce a new distribution f{R,T) = \I't'(-R)^E'(-R, t). We obtain a 
non homogenous Fokker-Plank equation 

- ^f{R, t) = -|^VV(i?, ^) + ■ [vD{R)fiR, r)] + [El{R) - Eo]f{R, r) , (2.212) 



where 
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is the 3A dimensional "drift velocity" and 
is the "local energy". 

The imaginary time evolution for f{R, r) can be conveniently written in terms of a 
modified Green's function 

f{R, r + Ar) = J G{R, R', Ar)/(i?', T)dR' . (2.215) 

Comparing with Eq. fl2.193p it is immediately found that 

G{R, R', At) = G{R, R', Ar)^^ . (2.216) 

It is shown in Appendix [D] that the importance sampling makes the diffusion driven by 
the drift velocity, that carries the walkers along in the direction of increasing 

G4R,R',Ar) = (^-^)-e"-^[^-^'-^-H^ . (2.217) 
The branching factor now contains the local energy instead of the potential energy 

El{R) + El(R') At 



GbiR,R',AT) = e V 2 J . (2.218) 

If the trial wave function is sufficiently accurate, the local energy remains close to the 
ground-state energy throughout the imaginary time evolution. 

As far as the expectation value of the operator O is concerned, from the last term of 
Eq. fl2:2T0|) it follows 

IdR^^^f{R,r) 

Using the central limit theorem, (O) can be computed by sampling the configurations 
from f{R, t) 

2.6.4 Sign problem 

In order to project out the ground state of a given Hamiltonian, the DMC algorithm 
implies a diffusive process, whose starting distribution of walkers is given by the trial wave 
function. Hence, for the diffusion interpretation to be applicable, \E't must be positive 
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definite in tlie wliole configuration space. Tliis is tlie case, for example, of many-bosons 
system, whose ground-state wave function is positive definite. 

The ground state of a fermionic system on the other hand, is described by an antisym- 
metric wave function, to which a probability distribution interpretation cannot be given. 
Let us describe this issue in more details. 

It can be proven that the ground state \E'o(-R) of a regular hamiltonian H is node less. 
Hence, from a strictly mathematical point of view, the search of an antisymmetric ground 
state, '^q{R), corresponds to the search of an excited state of the many-body hamiltonian. 
In terms of the energy eigenvalues, this corresponds to 

Eo < , (2.221) 

where Eq and Eq are the ground-state energies for the bosonic and the fermionic system 
described by H [77j . 

Expanding the trial antisymmetric wave function in terms of eigenstate of the hamil- 
tonian and choosing the energy offset to be Eq, in the limit of large imaginary time we 

get 

lim e-"^^T{R) = lim [ V Cne~^^--^^^^^n{R) + 4^q{R) + • • • 1 (2.222) 

n 

The sum over n runs over the bosonic eigenfunctions of the hamiltonian having a smaller 
energy than Eq] when r — )■ oo these terms diverge. The dots indicate the converging 
term, i.e. the eigenfunctions (both bosonic and fermionic) with energies larger than E^ 
that are exponentially suppressed with respect to '^^{R). 

The exponentially growing component along the symmetric ground state does not 
affect the expectation of the Hamiltonian. Because of the orthogonality between antisym- 
metric and symmetric wave functions, it turns out that 

(^4^^ {Mjym^^y j dRm{R,T)[Hm^]{R) 

^ ' (^ol^o) ™ (*WI^t) JdR^{R,T)^^{R) 
^ ^.^ /^i?[E.Cne-(^"-^o")-vI/„(i?) + c^vl/;f(i?)][gvl/j^](i?) 
/^^[EnCne-(^"-^")-v]>„(i?) + c^^^{R)]^^,^{R) 

E^ . (2.223) 



J dRc^^^{R)[H^^]{R) 



f dRc^^^{R)^^{R) 

However, the variance of the DMC estimate for the energy expectation value 0"^^ = 

|([i^\E'^]^) — ([^f^'^])^! is exponentially diverging. In particular, the bosonic components 
dominates the second term of the variance, as the orthogonality which eliminates the 
symmetric contributions does not apply in the following integral 

- Jr. T^R^sWrnW) • 
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We are left with the contradictory statement that the energy converges to exact eigen- 
value with an exponentially growing statistical error: the signal to noise ratio exponen- 
tially decays. 

In order for the DMC to be used for fermionic systems, it is possible to artificially 
split the configuration space in regions within which the sign of the trial wave function 
does not change. The 3^4 — 1 dimensional subset of the configuration space where the 
trial wave function vanishes is denoted as "nodal surface". In the fixed node approximation 
|14U) . during the diffusion process the walkers crossing the nodal surface are dropped. In 
other words, the nodal surface of the ground-state is imposed on the system, as it defined 
by the constraint "^tIR) = 0. It can be proven fl41] that the energy obtained from a 
fixed node DMC simulation obeys a variational principle. It is important to note that 
this variational principle only applies to ground-state calculations while a much weaker 
variational principle holds for excited-state calculations |142 . 



In the case of nuclear hamiltonian, the overlap between walkers and the trial-wave 
function is complex and the sign problem turns into a phase problem. A generalization 
of the fixed-node approximation, the so-called "constrained-path" |143] was introduced to 
deal with complex wave functions. It amounts in constraining the walkers to diffuse in 
regions where the overlap with the trial wave function is positive. To this aim, a suitable 
choice of the drift terms, used in the earlier AFDMC calculations |144) is 

To avoid the signal to noise ratio exponentially decay, the constrained-path approximation 
is realized by imposing 

^^IMS> 0. (2.226) 

Thus, the walkers having an overlap with the trial wave function that after a diffusive 
step changes sign are dropped. 

Another approach followed to put under control the sign problem is the "fixed-phase" 
approximation, introduced do deal with hamiltonian containing a magnetic field [145]. 
The walkers are forced to have the same phase as the importance function ^t- The drift 
term is given by 

It can be shown that an additional term in the branching has to be considered; in particular 
only the real part of the kinetic energy contribution to the local energy has to be kept 



V^^t{R) „ \V^^t{R) 
^ — - — )• Re ^ — - 

^t{R) L mT{R) J ■ 



(2.228) 



Both for the constrained-path and the fixed-phase approximations an accurate trial 
wave functions would be needed. While in GFMC calculations | il46j the full operator 
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structure of the CBF is taken into account, AFDMC trial wave function only encloses 
pure central Jastrow correlations, which is positive definite. Hence for AFDMC, the nodal 
structure of the nuclear matter wave function is entirely given by the Slater determinant 
of plane waves. 

It is relevant to the purpose of the sign problem discussion to remark that using 
constrained-path approximation, the DMC algorithm does not necessarily provide an 
upper bound in the calculation of energy |147| . Moreover, it has not been proved that 
the fixed-phase approximation gives an upper bound to the real energy. 

For further details concerning constrained path and fixed phase approximations the 
reader is referred to the original papers and to the exhaustive discussion reported in the 
PhD Thesis of Paolo Armani [138j . 



2.6.5 Spin-isospin degrees of freedom and auxiliary fields 



(2.229) 



The method we have described needs to be generalized to account for spin-isospin 
degreed of freedom, that are of major importance in nuclear few- and many- body systems. 

Let us start from an example: within the GFMC approach the eight spin configuration 
of the nucleus (we neglect for the moment the isospin) are represented by |148] 

V ^miR) J 

Each coefficient Oq, represent the amplitude of a given many-particles spin configura- 
tion; for instance 

amiR) = (m \'H) . (2.230) 
The many-particles spin configuration space is closed under the action of the operators 
contained in the hamiltonian. For example, applying ai2 = 2P{^ — 1 yields 

(^miR) 



Ol2?H) 



'^o-miR) 



(R) 



\ 



amiR) 
amiR) 



(2.231) 



Since the total charge is conserved, for the isospin of the we have pnn, npn, or 
nnp; thus, the vector describing the whole spin-isospin structure has 24 entries. In the 
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GFMC algorithm the the imaginary time evolution of Eq. fl2.215p is need to applied to 



each of the 2^- 



Z\{A-Z) 

Eq. (I2.193P generalizes to 



spin-isospin configurations and the the imaginary time evolution of 



a„(i?, r + Ar) = ^ j dR'Go,p{R, R', Ar)a^(i?', r) . (2.232) 



The Green's function depends on the spin-isospin configuration 
Gc,^{R,R',At) = {R,a\e~^^-^^^'^^\R',P) . 



(2.233) 



In order to deal with system having a large number of protons and neutrons, like for 
example medium-heavy nuclei or nuclear matter, GFMC does not seem to be a feasible 
approach. The idea of AFDMC consists in using a single-particle wave function, instead 
of the many-particle wave function of GFMC. For comparing the two methods, the spin 
structure of in AFDMC approach is 



= cll t)i + cii ;)i ® 4it)2 + 4i;)2 ® 4lt)3 + 4U) 



^2.234) 



where complex coefficient c" denotes the amplitude for the i — th particle to have spin 
state a. Taking also the isospin degrees of freedom into account, it can be readily shown 
that the dimension of the vector describing a system with A nucleons is AA. 

Already for such a small nucleus like the dimension of the spin-isospin structure 
of ADMC is a factor 2 smaller than the one of GFMC. The gain in computational time of 
AFDMC with respect to GFMC becomes enormous for larger system. However, GFMC 
is still the best or at least one of the best available method for dealing with hard-core 
potentials like the Argonne Vis- The spin-orbit terms and three nucleon forces have not 
been included in AFDMC algorithm yet, with the notable exception of PNM. 

The main concern of the single-particle wave function is that it is not closed with 
respect to the application of a quadratic spin (or isospin) operator. Let us again apply 
the operator ai2, as we did in Eq. (I2.23ip 



ai2\'H) = 2 



41 t)i + 4\ ® \c\\ t)2 + cji 1)2] ® [4i 1)3 + 4\ 1)3 

c[\ t)i + cil Dil ® \cl\ t)2 + 4| 1)21 ® \cl\ Da + 4\ 4)3 



(2.235) 



The resulting sum of two single-particle wave functions cannot be expressed as a single 
particle wave function. Therefore, if using the standard DMC algorithm, the imaginary- 
time propagator generates a sum of single particles wave functions at each time step. This 
would be catastrophic, as the number of single particle wave function would soon become 
enormous. 

The idea of AFDMC is to use the Hubbard-Stratonovich transformation to reduce the 
spin-isospin dependence from quadratic to linear, making the use of single particle wave 
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functions feasible. In the following we will show how this is done for Argonne potentials 
incorporating the first six operators. The inclusion of the spin-orbit term is possible in 
the case of PNM only and it has beed described at length in Ref. [149j and |138) . The 
first six components of the Argonne potential can be rewritten as 

6 

^ = E E = "^si + VsD (2.236) 

i<j p=i 

where the spin independent and spin dependent contribution read 

^SD = 2 X] ^i^^'iajP^iP 2 ^ (^io^^lajpf^ipT'i ' + 2 5Z ^li^^ ' ■ (2.237) 

io,j/3 ia,jl3 i,j 

As usual, Latin indices label nucleons, while the Greek indices stay for cartesian compo- 
nents. From Eqs. (fL20|) and (fLTG]) . the 3 A x 3 A matrices A" and A"'', and the A x A 
matrix A'^ in the case of the Argonne potential are readily seen to be |25 ) .[22 j 

= v''{nj)5o.p + t;*(r,,-)(3f°fj - 5^p) , 
= v'^^in,)6^p + t;*-(r,,)(3f^.rg - 5.^) , 

=^^i^^,)■ (2-238) 

The matrices A are vanishing on the diagonal, as there is no self-interaction in the po- 
tential. Moreover, since they are real and symmetric, they have real eigenvalues and 
orthogonal eigenstates, given by 



(JT 
iOL ' 



J/3 



Y.^Wn,=K^l,- (2.239) 

i 

It is convenient to normalize the eigenstates as follows 

for p = cr, ar, r. Using the last equation, cTia can be expanded as 

= E (E^-J/^^i/')^,-- (2-241) 
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Substituting tlie latter result in the first term of Eq. fl2.237p yields 

^ E { [ E ( E ^-.7^ [ E ( E ^-,^^^") J } 



^ E { [ E ( E [ E ( E V'^.mo^m) 

ia n l,S m 1,6 

^E(E<.M^'^)'^n- 

n l,S 

The expansion of CTjaTj reads 

(T^aT^ = E ( E Cj/^^i/^^j) C^ia ' 

iquently the second term of Eq. (I2.237P can be rewritten as 



(2.242) 



(2.243) 



consei 



^E(Eo^^O'^" 

n 1,5 



(TT 



^ia,jfS 



m l,S 



.err 
ia 



m l,S 



m,ia 



The same steps can be followed for the term with A'^'^\ It is worth introducing 
of operators written in terms of eigenvectors of matrices A 

j 

In terms of these operators, the spin dependent part of the potential reads 

^ 3A 3 3A 3 A 

= I J2{o:rK + 1 E J2io-n7 + \ E E(o;«)'^n 

ra=l a=l ?i=l a=l n=l 



(2.244) 



a new set 



(2.245) 



(2.246) 
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Since the operators in general do not commute, it turns out tliat 



(2.247) 



With the symbol 0„ we denote the O^, the O^^ and the O^. 

It is now possible to use the Hubbard-Stratonovich transformation, that for a generic 
operator O and a parameter A is defined by 



1 



dxe 2 



(2.248) 



The propagator of Eq. f l2.247p can be recasted in the following form 



(2.249) 



We define a walker to be the 3 A spatial coordinates and I'A spinors, cf. Hence, within 
the AFDMC approach, the spin-isospin coordinate 5* has to be added to the spatial 
coordinate R. The imaginary time evolution (compare with Eq. (I2.193P and Eq. (I2.232P ) 
reads 

^(i?, S,t + At) = j dR'dS'GiR, S, R', S', At)^{R', S', t) , (2.250) 

where the AFDMC Green's function, which includes the integration over the auxiliary 
fields, is given by 



GiR, S, R', S', At) 



m 



2'Kh?AT 
15A 



3A 12 



X 



(Sin 



n=l 



S').. 



(2.25r 



An important point to make is that the operator 0„ contains a sum over particle index 
j, as can be seen from Eq. fl2.245p . However, these operators commute and we can 
conveniently represent the exponential of the sum as a product of exponentials, each 
rotating only one single-particle state. Therefore, the application of the operator to 
the spin-isospin state |S"), turns into a product of independent rotation. For the spin 
rotation, generated by O^, one has 



o 



S') 



cli t)i + cii i)] 



Rotating the j — th single particle state amounts in a change of the coefficients 

(2.253) 



c\\t)A + c\\l)A . (2.252) 
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In order to give an explicit expression to and , tlie following identities have to be 
used 



cos(|a|) + i sin(|a|) 
a 



(2.254) 



cosh(|a|) + ^ _^ sinh(|a|) 



(2.255) 



where in our case the vector a is given by (we omit the index a for brevity) 

a = ^/\Xn\ATXn^n,j ■ 

When A„ < 0, exploiting Eq. (I2.254p . we find 



t 



(3 



cosh(|a|) + ^ sinh(|a|) +cjsinh(|a|) " 



cosh(|a|) + 



lb'' ■ 



sinh(|a|) +asinh(|a| 



1^. 



(2.256) 

(2.257) 
(2.258) 



n,J I 



On the other hand, if A > 0, making use of Eq. fl2.255p the transformed coefficients 



read 



cos(|a|)+i S''' sin(|a|) +rfsin(|a|) 



n,J I 



cos(|a|)+^ J^''' sin(|a|) +cjsin(|a|) 




(2.259) 
(2.260) 



Note that if the integral over the auxiliary fields was computed using the standard 
methods, like the Simpson rule, we would be left with a sum of rotated spinors, one for 
each sampled value of the auxiliary fields x„. 

In the first realizations of the AFDMC algorithm a discrete version of the Hubbard- 
Stratonovich transformation, due to Kooning, was implemented. It essentially consists in 
replacing the Gaussian by a three point weighted sum. With a probability depending on 
the weight, only one over these three values has to be used for rotating the spinors. 

In more recent works, following the spirit of the Monte Carlo algorithm, the Gaussian 
has been considered as a probability distribution. One value is sampled directly from the 
Gaussian and used to rotate the spin-isospin degrees of freedom of the walkers. 

A physical interpretation can be (and has been) given to the auxiliary fields. As can be 
consistently explained in chiral perturbation theory, nuclear interactions can be explained 
in terms of pion exchanges. Within the Born-Oppenheimer approximation, the light pions 
are the fast degrees of freedom, coupled with the slow and more massive nucleons. If the 
fast meson field is integrated out to give a potential, keeping the nucleonic coordinates 
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fixed, and solving for its ground-state energy, then meson coordinate corresponds to the 
Hubbard-Stratonovich auxiliary fields. 

Importance sampling can be implemented to the integral over the auxiliary fields. 
The overlap of the walker with \E't is not generally picked around x„ = 0. Hence, instead 
of sampling from the Gaussian, it is more efficient to sample values of Xn where the 
trial wave-function is thought to be large. One way consists in shifting the Gaussian, 
introducing a drift term analogous to the one used for the spatial coordinates. For the 
detailed calculations of the drift term for the Hubbard-Stratonovich variables, the reader 
is referred to Refs. |25l [I38l ITI9] . 

The agreement between Green Function Monte Carlo (GFMC) and AFDMC energies 
of neutron drops, obtained using the Argonne fg plus UIX hamiltonian, discussed in 
Ref. |150| . supports the validity of PNM calculations carried out within AFDMC with the 
Argonne fg model. The highly accurate GFMC method has been used to study neutron 
matter properties in both the normal |134j and superffuid [151] phases. 

Moreover, using a fixed-phase like approximation, AFDMC also yields results in very 
good agreement with those obtained from Green Function Monte Carlo (GFMC) calcu- 
lations for light nuclei |22] . 



Chapter 3 

Three-body potential in nuclear matter 



3.1 UIX potential within FHNC/SOC approach 

Within CBF approach, the expectation value of a three-body potential, e.g. the UIX 
model, reads 

^ ^ (A-3)!3! {^>l\PP\^>,) ■ ^ ■ ^ 

Let us write V123 as a sum of spin-isospin three-body operators multiplied by scalar 
functions, depending on the relative distances only 

K23 = ^123^^f23 ■ (3.2) 

P 

As for the case of the two-body distribution functions (7^2; is useful to define three- 
body distribution functions g\2Zi reflecting the operatorial structure of V123 

- 3)! p3 J *|,(X)^t^^o(x) ■ ^ ■ ^ 

Hence, analogously to Eq. (I2.113P relative to the two-body potential, the expectation 
value of V123 can be written as 



A 3! 

p 



(3.4) 



Neglecting the Abe diagrams, as in Eq. f l2.102p . the three-body distribution function 
can be approximated by a product of the partial two-body distribution function of Eq. 
f imrj) and (EHZD, denoted by Z^y in Ref. yij, 

9i2z=Y^ 5Z ^^l'(^i2)<y(n3)^r"'V23)x 

ex p' ,p" 

CTr,23^ [dUotAOti + OtiiOt Of;"} + Ot^iOt O?;'})] (3.5) 
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It can be noted from the previous equation, that the FHNC/SOC approximation has 
been exploited as at most two-operators arrive at a given point. Moreover, for the sake 
of brevity, vertex corrections arising from separable diagrams have not been reported. 





Figure 3.1: Cluster diagrams contributing to the expectation value of V'^'". 



The diagrams involved in the FHNC/SOC calculation of the expectation values of the 
V'^'^ term of the UIX potential are depicted in Fig. 13.11 The thick lines represent the 
potential, while dashed and wavy lines correspond to the partial two-body distribution 
functions; vertex corrections, although included in the calculations, are not shown. Be- 
cause of the symmetry properties of the wave function, we can restrict our analysis to the 
permutation (3 : 12) . The other three permutations are accounted for in the symmetry 
factor appearing in front of each of the following expressions 



ex p' 



ex p" ,p" 



nP n'^ nP PTr 

' 13 yxy,12 yx'y',13 yx"y",23^ 



V''{i:12){01^,0l,} (3.8) 



ex p' ,p"' 
p2 r 

=uJ2Yl J ^^^12C^^^13 9xy,12 9x'y', 13 f x'V',23 >< 
ex p' ,p" ,p"' 

CTr[v^^-(3 : 12)raO?;, o£} + 0?;{df2, Og} + OgjOf^, O?;})^ 



(3.6) 



ra^^) =jY.J2 I ^'^i2^'^i3^?:'VM2<,i3^?S'.23CTr[r2-(3 : 12){d{„di}] (3.7) 



(3.9) 
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The computation of all of diagrams (13. 21 a). (13.21 b) and (13.21 c) and all diagrams of Fig. 
13.21 is outlined in Ref. [T7], while the contribution of diagram (3.d), involving three non 
central correlations was first taken into account by the authors of Ref. |19) . 

In addition to the what listed in Eq. (13. Sp . SOR contributions to the three-body 
distribution function, represented by double wavy lines in Fig. 13.21 and denoted by Y^y, 
are included in the calculation of the repulsive term of the UIX potential 



3 3 




3 3 




Oa) =y E / d'n2d'ns9:.y.^,2 gly.iz 9l'y',2,v\^ 12) (3.10) 

ex 

=y E / d'ri2d'ns 9:y,,s 9:>y>,23V''iS : 12) (3.11) 

ex 

Dc) =p' f d'rud'n, gly^,^ Y^,y.^,, gl"y'' ,2,V\?> : 12) (3.12) 

ex 

=y E E / d'n2d'n, 4,,, ,3 eS>'A^v^(3 : 12) . (3.13) 



ex p,p',p" 



3.2 Density dependent effective potential 

As shown in Fig. II. 5^ the inclusion of the UIX three-body potential in the hamilto- 
nian considerably improves the theoretical estimates of the energies of the ground and 
low-lying excited states of nuclei with A < 12 . However, for nuclei heavier than ^H, 
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some discrepancies witli experimental data persist; moreover the empirical equilibrium 
properties of nuclear matter are not correctly reproduced. This problem can be largely 
ascribed to the uncertainties associated with the description of three-nucleon interactions, 
whose contribution turns out to be significant. 

3.2.1 Derivation of the effective potential 

Our work is aimed at obtaining a two-body density-dependent potential vulp) 
that mimics the three-body potential. Hence, our starting point is the requirement that 
the expectation values of V123 and of Vulp) be the same: 

(3.14) 

implying in turn (compare to Eqs. fl2.113p and (13. 4p ) 

E f / ^^3^4 9^23 = E <2ip) 9I2 ■ (3.15) 

A diagrammatic representation of the above equation, which should be regarded as the 
definition of the Vuip), is shown in Fig. 13.31 The graph on the left-hand side represents 
the three-body potential times the three-body correlation function, integrated over the 
coordinates of particle 3. Correlation and exchange lines are schematically depicted with 
a line having a bubble in the middle, while the thick solid lines represent the three-body 
potential. The diagram in the right-hand side represents the density-dependent two-body 
potential, dressed with the two-body distribution function. Obviously, has to include 
not only the three-body potential, but also the effects of correlation and exchange lines. 

In Section [XT] we have examined the left-hand side of Eq. (l3.15p . that has been evalu- 
ated in [1^ within the FHNC/SOC scheme. Here we discuss the derivation of the explicit 
expression of the two-body density- dependent potential appearing in the right-hand side 
of the equation. The procedure consists of three different step, each corresponding to a 
different dressing of the diagrams involved in the calculation 

For each of these steps the final result is a density- dependent two-body potential of 
the form 

vM = J2^''ip,rM, (3.16) 

p 

where, depending on the step, the vP{p,ri2) = Vi2{p) can be expressed in terms of the 
functions appearing in the definition of the UIX potential, the correlation functions and 
of the Slater functions. 

Step I. Bare approximation 

As a first step in the derivation of the density-dependent potential one integrates the 
three-body potential over the coordinate of the third particle 

V}^{P) = ^ I dX3V,23. (3.17) 
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Figure 3.3: Diagrammatic representation of Eq. fl3.15p : tlie two-body density-dependent 
potential includes the effects of both the bare three-body potential and the correlation 
and exchange lines. While g2 dresses the line joining particles 1 and 2, the dressing being 
depicted by a line with a big bubble in the middle, dresses the lines 1 — 2, 1 — 3, and 



Diagrammatically the above equation implies that neither interaction nor exchange lines 
linking particle 3 with particles 1 and 2 are included. Only the two-body distribution 
function is taken into account in the calculation of the expectation value of Vus 



Note that only the scalar repulsive term and one permutation of the anticommutator 
term of the three-body potential provide non vanishing contributions, once the trace in 
the spin-isospin space of the third particle is performed. 

As shown in Fig l3.8[ the contribution of the density-dependent potential to the energy 
per particle of SNM and PNM {v^2\p)) more repulsive than the one obtained from 
the genuine three-body potential UIX. Thus, the scalar repulsive term is dominant when 
the three-body potential is integrated over particle 3. 

Step II. Inclusion of statistical correlations 

As a second step we have considered the exchange lines that are present both in (7123 
and gi2. Their treatment is somewhat complex, and needs to be analyzed in detail. 

Consider, for example, the diagram associated with the exchange loop involving par- 
ticles 1, 2 and 3, depicted in Fig. 13. 4[ Its inclusion in the calculation of the density- 
dependent two-body potential would lead to double counting of exchange lines connecting 
particles 1 and 2, due to the presence of the exchange operator P12 in gi2. This problem 
can be circumvented by noting that the antisymmetrization operator acting on particles 



2-3. 




(3.18) 



94 



Three-body potential in nuclear matter 



1, 2 and 3 can be written in tlie form 

1 ~ Pi2 — Pi3 ~ P23 + -P12-P13 + P13P23 = 

(1 - As - P23) X (1 - A2) , (3.19) 

in wliicli tlie excliange operators contributing to the density-dependent potential only 
appear in the first term of the right-hand side. 

On the other hand, the second term in the right-hand side of Eq. fl3.19p only involves 
the exchange operators P12, whose contribution is included in gi2 and must not be taken 
into account in the calculation of vuip)- 

Two features of the above procedure need to be clarified. First, it has to be pointed out 
that it is exact only within the SOC approximation that allows one to avoid the calculation 
of commutators between the exchange operators P13 and P23 and the correlation operators 
acting on particles 1 and 2. The second issue is related to the treatment of the radial 
part of the exchange operators. Although it is certainly true that one can isolate the 
trace over the spin-isospin degrees of freedom of particle 3, arising from P13 and P23, 
extracting the Slater functions from these operators is only possible in the absence of 
functions depending on the position of particle 3 [_128j. 




Figure 3.4: Three particle exchange loop. 



However, this restriction does not apply to the case under consideration, as both the 
potential and the correlations depend on ri3 and r23. As a consequence, retaining only 
the Pi3 and P23 exchange operators involves an approximation in the treatment of the the 
Slater functions, whose validity has been tested by carrying out a numerical calculation. 

By singling out the radial dependence of the exchange operators, and by computing 
the inverse of the operator (1 — P{^^), where P^'^ denotes the spin-isospin part of Pij, it 
is possible to find a "Slater Exact" density-dependent potential ff2^ (p) whose calculation 
does not involve any approximations concerning the Slater functions. It can be easily 
verified that 

(1 - P-ij^r = (3-20) 
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Scalar contribution 




0.1 0.2 0.3 0.4 0.5 



Anticommutator contribution 




0.1 0.2 0.3 0.4 0.5 

P (fm-^) 



Figure 3.5: Contributions of tlie density-dependent potential to the energy per particle 
(see Eqs. (I3.2ip and (I3.22p ). arising from the scalar term of UIX (upper panel) and from 
the anticommutator term (lower panel). 



thus, the "Slater Exact" density- dependent potential can be rewritten in the form 



1 — Z"' 

^12 



At (4^13^3 - 4) + At (4^13^23 - 4)+ 



ATAT(^12^13^23 



'12^13, 



12^13^23 



'12^23 J 



(3.21) 



Note that in the above equation we have omitted all correlations functions, whose 
presence is irrelevant to the purpose of our discussion. The density-dependent potential 
obtained from Eq. (l3.2ip must be compared to the one resulting from the approximation 
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discussed above, which (again neglecting correlations) leads to the expression 

v',2^ip) = dxsVusil - P^7^l3 - , (3.22) 

where "S. A." stands for Slater Approximation. We have computed (f f2^'(p)) and (f f2"^ (p)) 
for SNM within the FHNC/SOC scheme, for both the scalar and the anticommutator 
terms of the UIX potential. 

The results, plotted in Fig. 13.51 clearly show that Eq. fl3.22p provides an excellent ap- 
proximation to the exact result for the exchanges of Eq. fl3.2ip . Hence it has been possible 
to use Eq. fl3.22p also to compute the contribution coming from the commutator of the 
UIX potential, avoiding the difficulties that would have arisen from an exact calculation 
of the exchanges. 

The second step in the construction of the density-dependent potential is then 

vuip)^v'u''ip) (3.23) 

which is a generalization of the bare potential of Eq. f l3.17p . 

Figure 13.81 shows that taking exchanges into account slightly improves the approxi- 
mation of the density- dependent potential. However the differences remain large because 
correlations have not been taken into account. 



Step III. Inclusion of dynamical correlations 

The third step in the construction of the density-dependent potential amounts to 
bringing correlations into the game. We have found that the most relevant diagrams are 
those of Fig. 13.61 




Figure 3.6: Diagrams contributing to the density- dependent potential. The dashed lines 
with diamonds represent the first order approximation to gb^osF i^ij) ^ discussed in the text. 
Only diagrams with at most one operator attached to a given point are taken into account. 
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Note that, in order to simplify the pictures, all interaction lines are omitted. However, 
it is understood that the three-body potential is acting on particles 1, 2 and 3. Correlation 
and exchange lines involving these particles are depicted as if they were passive interaction 
lines. Moreover, in order to include higher order cluster terms, we have replaced the scalar 
correlation line /^^-^ with the Next to Leading Order (NLO) approximation to the bosonic 
two-body correlation function: 

f^' ^ 9Zse''{r^J) = f-' [l + pj dr,hsh23) . (3.24) 

The full bosonic gtoseif^ij) ot gdd{fij) niight be used instead of the NLO approximation. 
However, including higher order terms would have broken our cluster expansion. The 
correction to of Eq. f l3.24p . whose diagrammatic representation is displayed in Fig. 
13. 7[ can indeed be considered to be of the same order as the operatorial correlations. 

Figure [3l6] shows that the vertices corresponding to particles 1 and 2 are not connected 
by either correlation or exchange lines. All connections allowed by the diagrammatic 
rules are taken into account multiplying the density-dependent potential by the two-body 
distribution function, according to the definition of Eq. (l3.15p . 

We have already discussed the exchange lines issue, coming to the conclusion that 
only the exchanges P13 and P23 have to be taken into account. This is represented by the 
second diagram, where the factor 2 is due to the symmetry of the three-body potential, 
that takes into account both P13 and P23. 
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Figure 3.7: NLO approximation to the bosonic two-body correlation function. 

The explicit expression of f 12^^"* (p) obtained including the diagrams of Fig. 13.61 can be 
cast in the form 

v[2'\p) = U d^sVi2s ki;^(ri3)^7L^e^(r23) 



3 _ 

X (1 - 2A74) + ^9Zse''{ri3)fc{r23)Kr2s)] , (3.25) 
where f{r2z) denotes the sum of non central correlations 

/>23) = X^/^(r23)Of,-. (3.26) 
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Figure 3.8: Contributions of tlie density-dependent potential to the energy per particle 
of SNM (a) and PNM (b), compared to the expectation value of the genuine three-body 
potential UIX: (1^123) j A. 



Note that, in principle, an additional term involving the anticommutator between 
the potential and the correlation function should appear in the second line of the above 
equation. However, due to the structure of the potential it turns out that 

y rfX3{V^i23, /(r23)} = 2 y rfX3V^i23/(r23) • (3.27) 

The calculation of the right-hand side of of Eq. fl3.25p requires the evaluation of the 
traces of commutators and anticommutators of spin-isospin operators, as well as the use 
of suitable angular functions needed to carry out the integration over r3. 

As for the previous steps, we have computed the contribution of the density-dependent 
potential vl-l \P) to the energy per particle. The results of Fig. 13.81 demonstrate that 
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the density-dependent potential including correlations is able to reproduce the results 
obtained using genuine three-body UIX to remarkable accuracy. 
To simplify the notation, at this point it is convenient to identify 

vM^vg"\p). (3.28) 

Note that the above potential exhibits important differences when acting in PNM 
and in SNM. For example, in SNM v^{p,ri2) 7^ for p = 1, cr 12T12, S12T12, while in PNM 
v^ip, ^12) 7^ for p = 1, (T12, 5'i2, as shown in Figs. 13.91 
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Figure 3.9: Different channels of the effective density dependent potentials calculated for 
SNM (a) and PNM (b) for p = 0.16 fm^l 
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Table 3.1: Uncorrelated kinetic energy per particle of 14, 38, 66 and 114 neutrons at 
p = 0.32 MeV compared with the limit A — oo. 

A = 14 A = 38 A = 66 A = lU oo 
E/A{MeV) 56.51 53.50 55.43 56.58 55.71 



3.2.2 Numerical Calculations 

A constrained simulated annealing optimization, described in Section 12. 5. 4^ has been 
performed, by imposing the sum rules for the kinetic energy and for the scalar two-body 
distribution function. In particular the difference between the Pandharipande-Bethe (PB) 
and the Jackson- Feenberg (JF) kinetic energies has been forced to be less than 10 % of the 
Fermi Energy Tp of Eq. f l2.174p . while the sum rule fl2.39p for g'^iru) has been satisfied 
with a precision of 3 %. 

In our calculations we have optimized the variational paremeters for four different 
Hamiltonians, each corresponding to different potential terms: Argonne Wg, Argonne 
v'g + UIX, Argonne fg, and Argonne Vg + UIX. 

The energy per particle of SNM and PNM computed adding to the two-body potentials 
Argonne fg and Argonne v'q the density-dependent potential of Eq. (13.161) . have been 
compared to the results obtained using an hamiltonian with the same two-body potentials 
and the Urbana IX three-body interaction model. In order to show how much the density 
dependent potential differs from the original UIX, we compute the expectation values 
of these potentials with the same correlation functions, i. e. those resulting from the 
calculation with the genuine three-body potential; no optimization procedure has been 
performed for the density-dependent potentials. 

Both calculations have been consistently carried out within the FHNC/SOC scheme. 

It is worth noting that our simulated annealing constrained optimization allows us to: 

i) reduce the violation of the variational principle due to the FHNC/SOC approximation; 

ii) perform an accurate scan of the parameter space. As a consequence, our FHNC/SOC 
calculations provide very close results to those obtained via Monte Carlo calculations, as 
shown in Figs. I3.10l and 13.11^ to be compared with those of Ref. [22j where the agreement 
between FHNC and Monte Carlo methods were not nearly as good. 

Auxiliary Field Diffusion Monte Carlo (AFDMC) approach 

In order to check the validity of our variational FHNC/SOC calculations, we carried 
out AFDMC simulations for both PNM and SNM. 

We have computed the equation of state of PNM and SNM using the AFDMC method 
with the fixed-phase like approximation. We simulated PNM with A = 66 and SNM 
with y4 = 28 nucleons in a periodic box, as described in |152] and |153] . The finite- 
size errors in PNM simulations have been investigated in [T52J by comparing the Twist 
Averaged Boundary Conditions (TABC) with the Periodic Box Condition (PBC). It is 
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remarkable that the energies of 66 neutrons computed using either twist averaging or 
periodic boundary conditions turn out to be almost the same. This essentially follows 
from the fact that the kinetic energy of 66 fermions approaches the thermodynamic limit 
very well, as can be seen in Table 13.11 The finite-size corrections due to the interaction 
are correctly estimated by including the contributions given by neighboring cells to the 
simulation box|i25j. From the above results for PNM we can estimate that the finite-size 
errors in the present AFDMC calculations do not exceed 2% of the asymptotic value of 
the energy calculated by using TABC. The finite-size effects of SNM calculations can be 
estimated from the difference of the energies of PNM obtained with 14 neutrons and the 
TABC asymptotic value, which is of the order of 7%. Although a calculation with 132 
nucleons would be more accurate, performing such a heavy simulation does not appear 
to be justified in the context of a preliminary model of density dependent potential, and 
also in view of the fact that, at present, we are not able to simulate SNM with spin-orbit 
interactions. 

The statistical errors, on the other hand, are very small and in the Figures are always 
hidden by the squares, the triangles and the circles representing the AFDMC energies. 

PNM equations of state 

In the PNM case (see Fig. I3.10p . the EoS obtained with the three-body potential UIX 
and using the density-dependent two-body potential are very close to each other. For 
comparison, in Fig. 13.101 we also report the results of calculations carried out including 
the two-body potential only. In our approximation, with the exception of the line with 
diamonds of Fig. 13. 6[ we have neglected the cluster contributions proportional to . One 
could then have guessed that the curves corresponding to the UIX and density-dependent 
potential would have slightly moved away from each other at higher densities because, as 
the density increases, the contributions of higher order diagrams become more important. 
Probably, in this case a compensation among these second and higher order terms takes 
place. 

The density-dependent potential obtained in the FHNC/SOC framework has been 
also employed in AFDMC calculations. As can be plainly seen in Fig. I3.10[ the triangles 
representing the results of this calculation are very close, when not superimposed, to the 
circles corresponding to the UIX three-body potential AFDMC results. 

SNM equation of state 

In the EoS of symmetric nuclear matter, the above compensation does not appear to 
occur, as can be seen in Fig. 13.111 At densities lower than p = 0.32 fm~^, the curves 
resulting from UIX and the density-dependent potential are very close to one other, while 
for p > 0.32 fm~^ a gap between them begins to develop. 

The gap is smaller when the two-body potential Argonne v'g is used, but the reason 
for this is not completely clear. 

We have computed the saturation density po, the binding energy per particle -E(po) 
and the compressibility K = 9po{dE{p)/dpy for all the EoS of Fig. 13.111 The variational 
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Figure 3.10: Energy per particle for PNM, obtained using the density-dependent potential 
of Eq. f l3.17p added to the Argonne f g (a) and to Argonne f g (b) potentials. The energies 
are compared to those obtained from the genuine three-body potential and from the two- 
body potentials alone. 



FHNC/SOC results are listed in Table [3l2| while those coming from the AFDMC calcu- 
lation with v'q + Vuip) potential are: po = 0.17fm"^ Eq = -10.9 MeV and K= 201 MeV. 

The saturation densities are quite close to the empirical value po = 0.16 fm~^. For the 
genuine three-body potential this is not surprising, since the parameter Uq is chosen to 
fit the saturation density, as discussed in Section ll.2.2[ On the other hand, the fact that 
the density-dependent potential also reproduces this value is remarkable and needs to be 
emphasized. 

The binding energies obtained with vuip) are very close to those coming from UIX 
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Figure 3.11: Same as in Fig. IXTUl but for SNM. 



potential, but they are larger than the empirical value Eq = — 16MeV. 

As for the compressibility, the experimental value K ^ 240 MeV suffers of sizable 
uncertainties. However, also in this case the result obtained with the density-dependent 
potential differs from that obtained with the UIX potential by less than 5%. 



3.3 Chiral inspired three-nucleon potentials in nuclear 
matter 

The work described in this Section, based on Ref. is aimed at testing in nu- 

clear matter the different parametrization of the chiral inspired potentials of Ref. [45j . 
introduced in Section I1.2.2I 
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Table 3.2: Values for the saturation densities, the binding energy per particle, and the 
compressibility of SNM obtained from the variational FHNC/SOC EoS of Fig. 13.111 







< + v{p) 


^^8 + ^123 


v'^ + v{p) 


Po (fm-3) 


0.17 


0.16 


0.16 


0.15 


Eo (MeV) 


-11.3 


-11.2 


-10.3 


-10.3 


K (MeV) 


205 


192 


189 


198 



3.3.1 NNLOL contact term issue 

While the NNLOL chiral interactions provide a fully consistent description of the 
binding energies of and ^He, as well as of the scattering length ^a„rf, some ambiguities 
emerge when these interactions are used to calculate the nuclear matter EoS. 

For our purposes, it is convenient to rewrite the NNLOL chiral contact term of Eq, 
(I1.37P in the form 

Vi(3 : 12) = Vo^ri2Zo(ri3)Zo(r23) . (3.29) 

where the superscript r has a meaning that will be soon clarified. The radial function 
Zo{t) = / {4it) Zo{r) approaches the Dirac 5-function in the limit of infinite cutoff. 
Strictly speaking, the local version of Ve is a genuine "contact term" in this limit only, 
while for finite values of the cutoff it acquires a finite range. 

In addition to of Eq. f l3.29p . the chiral expansion leads to the appearance of six 
spin-isopin structures in the contact term. For example, the scalar contribution is 

Vi{3 : 12) = \/o^Zo(ri3)^o(r23) • (3.30) 

Within this context, the superscripts r and / identify the Tu and scalar contact terms. 

In Ref. [39j it has been shown that, once the sum over all cyclic permutation is per- 
formed, all contributions to the product between the potential and the antisymmetrization 
operator ^123 have the same spin-isospin structure. Therefore it is convenient to take into 
account just one of the contact terms. This result was obtained in momentum space, with- 
out the cutoff functions Fa. As a consequence, in coordinate space it only holds true in 
the limit of infinite cutoff. In particular, for V^(3 : 12) and V^J(3 : 12), it turns out that 

J2 Vo''S{r,3)S{r23)TuAi23 = -Yl Vo''6ir,,)6{r23)Ai23 , (3.31) 

cycl cycl 

making this two terms equivalent. The limit of infinite cutoff is crucial, because the radial 
part of the exchange operator, when multiplied by the Dirac 5-functions, is nothing but 
the identity 

e^'^--'--5(r,,) = 5(r,,) . (3.32) 

After the regularization, i.e. with the 5-function replaced by Zq, the proof is spoiled and 
the six different structures are no longer equivalent. 
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In PNM contact terms involving three or more neutrons vanish because of Pauli prin- 
ciple. On the other hand, the expectation value of the contact terms of the NNLOL 
potential can be different from zero. 

Let us assume that reproducing the binding energies of light nuclei and ^a„rf require a 
repulsive Ve ■ Then, one has to choose either c^^ < or > 0. In PNM, as 



PNM 



(3.33) 



it turns out that is attractive and repulsive. This means that fitting the binding 
energies and the n — d scattering length with either or alone leads to an ambiguity 
in the expectation value of the potential. 
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Figure 3.12: (Color online) Radial dependence of the Function Zo{r), appearing in 
Eq. fl3.29p . plotted for different values of the cutoff A. 



By expanding the cutoff function 

^A(g') = e 

one finds 



A4 



t>^(3 : 12) = K)''^i2[(5(ri3)(5(r23) + O 
fi(3:12) = K,^[5(ri3)5(r23) + 0(|^ 



implying that in PNM 



(V^i'-(3:12))p;vM = 0(|j). 



(3.34) 



(3.35) 
(3.36) 
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From the above equation it becomes apparent that the expectation value of the three- 
nucleon potential, as well as its sign ambiguity, is nothing but a a cutoff effect. Hence, it 
should be regarded as a theoretical uncertainty. Note that, since ~ A, then {Ve)pnm 
is of the same order of the next term in chiral expansion. 

To clarify this issue, let us consider a simple system: a Fermi gas of neutrons, in 
which correlations among particles are not present. The expectation value of the contact 
interaction reads 



IPNM _P 



A 



^V^cf j c^ri2rfri3Zo(ri2)^(ri3) x 



where A is the number of neutrons. The factor 1/2 includes the 1/3! arising from the 
unrestricted sum over particle indices 123, multiplied by a factor 3 from the cyclic per- 
mutations of the potential, all giving the same contribution. The Slater function i{rij) is 
given in Eq. f|2.59p . It can be easily seen that, if V^'^(l : 23) oc (5(ri2)5(ri3), then 

MpM = o. (3.38) 
A ^ ' 

Consider now a Fermi gas with equal numbers of protons and neutrons, where 



SNM 



J^o'' J rfri2rfri3Zo(ri2)Zo(ri3)x 
- ^^'(r23) + ^^(ri2)^(ri3)^(r23)) (3.39) 



and 



\^E/SNM 



'^^'^ j '^''i2^^i3^o(n2)^o(n3)x 



^(ri2)^ iir,,y i{runri3nr23) \ . . 

4 4 4 8 / ■ ^ ' ' 



In the limit of infinite cutoff the above equations imply 

(yTi2\FG o 

A "16^ ^° 
MpM = lpVo^. (3.41) 

As expected from Eq. fl3.3ip . the two contributions have opposite sign. 

We have computed the expectation values of Eqs. fl3.37p . fl3.39p and f l3.40p for different 
values of the cutoff A and density p = 0.16 fm~^. The results listed in Table 13.31 show 
that for PNM the larger the cutoff the smaller the expectation value of the three nucleon 
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Table 3.3: Cutoff dependence of the expectation values of the three body contact term of 
the NNLOL potential in noninteracting PNM. 



A (MeV) 


(^i"")^^MM(MeV) 


300 


9.15 


400 


5.95 


500 


3.60 


600 


2.15 


700 


1.30 


800 


0.81 


oo 






Table 3.4: Same as in Table IX^ but for SNM. 



A (MeV) 


(^;")f^A./M(MeV) 


{VM^M/AiUeV) 


300 


-2.61 


10.21 


400 


-3.61 


8.15 


500 


-4.37 


6.93 


600 


-4.87 


6.30 


700 


-5.15 


5.98 


800 


-5.30 


5.81 


oo 


-5.55 


5.55 



contact term. Note that for A = 500 MeV, the expectation value is still sizably different 
from the asymptotic limit. 

As far as SNM is concerned (see Table [33]), as the cutoff increases the possible choices 
of the three nucleon contact term tend to the asymptotic values of Eq. f l3.4ip . As in the 
case of PNM, the results corresponding to A = 500 MeV, are significantly different from 
the asymptotic values. 

We emphasize that the parameter Ce has not been included in this analysis, even 
though it is itself cutoff dependent. Unfortunately, the authors of Ref. [45j kept A fixed 
to 500 MeV. Had this not been the case, their fit to the experimental data would have 
resulted in a set of different constants Ce, corresponding to different values of A. It would 
have been interesting to extrapolate the expectation value of Ve to the limit of infinite 
A, where the cutoff effects associated with the regularization procedure are expected to 
vanish. 
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3.3.2 FHNC/SOC calculations 

Using the relations for the constants and the radial functions given in Eqs. fll.47p and 
( ll.48p . the computation of the diagrams of Fig. 13.11 with the V3 and V4 terms of both the 
TM' and NNLOL potentials is the same as that of V"^^ reported in Section 13.11 and in 
Ref. dZ]. 

Thanks to the identity 

(ai ■ fi3){a2 ■ r23)(ri3 ■ f23) = 7^{cri3 + 'S'la, 0-23 + S23} , (3.42) 

io 

the term Vi of Eq. f ll.37p . appearing in both the TM' and the NNLOL potentials, can be 
written in the form 

Wq Ti3r23 

Vi{3 : 12) = — — ^l/(ri3)t/(r23){ri3, r23}{(Ti3 + 5*13, (T23 + S23} ■ (3.43) 

60 [ri3 ■ r 23) 

Aside from the radial function, Vi is completely equivalent to V3, the anticommutator 
term of the UIX potential. Therefore, we were allowed to use again the results of Section 
Oand Ref. H?]. 

Furthermore, exploiting the identities 

(air23)(a2 ■ r23) = ^{5*23 + cr23, 0-13} 


(aa ■ ri3)(ai ■ ri3) = ^{^23, ^13 + ^13} , (3.44) 


we can rewrite the Vd term in a form that has again the same spin-isospin structure as 
the anticommutator contribution of the UIX potential 

Vd{3 : 12) = -^{n3,T23}[{ai3,(T23}V^^{n3,r23) + 

{Sl3, (T23}VS''{n3, r23) + {cri3, S23}V^^{r,3, r23)] , (3.45) 

where 

VD^irr3,r23) = Y{ri3)zo{r23) + Zoiri3)Y{r23) 
VD^iri3,r23) = Zo{ri3)T{r23) 

VS''{ri3,r23)=T{n3)zo{r23). (3.46) 

In conclusion, including Vd amounts to properly adding the above radial functions to 
those already appearing in V3. 

The Ve term of TM' is completely equivalent to Vr (see Eq. fll.49p ). This allowed 
us to use the results of Section 13.11 and Ref. fl7j for the diagrams of Fig. 13.21 The same 
holds true for the chiral contact term Ve in PNM, as {Tij)pNM = 1, while in SNM the 
calculation of Ve requires the evaluation of the diagrams of Fig. 13.11 
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Let us start from diagram fl3.1l a). whose analytic expression is given in Eq. (13.61) . As 
pointed out in Ref. [T^, integrating and tracing over the radial and spin-isospin variables 
of particle 3 leads to the appearance of an effective density dependent interaction 

E Kyy',i2iP)0{2 = pY.j ^^13 gly^r, Q^'y^a^GTi, [^^(3 : 12)] , (3.47) 



such that 



ma) =fpY.llj druA^al^'vAy^M ■ (3-48) 



ex p 



The subscripts xy label exchange patterns at the ends of the generalized correlation 
lines. In particular, dd correspond stands for direct-direct, de for direct-exchange, ee 
for exchange-exchange, and cc for incomplete circular exchange. The matrix = 
1, 3, 3, 9, 6, 18 has been defined in Eq. (I2.125|) . 

It turns out that the only nonvanishing term of the density dependent potential is 

VEyy',12ip) = ^0 pYI / ^^^^3 9xy,l3 9x'y' ,23^o{ri3)Zo{r23) ■ (3.49) 

ex 

As far as diagram f l3.1[ b) is concerned, from Eq. fl3.7p . it clearly follows that only 
r-type generalized correlation lines contribute. Hence 



y",12 



2:0(^3)2:0(^23). (3.50) 



Diagram fl3.1[ c) does not contribute to Ve, while the calculation of the spin-isospin 
trace of diagram (13. li d), with three generalized operatorial correlation, yields 

dSiy) = yW^o^'p' / drudr^,[~2gl, gl, g^, + 9g'^^ gl^ g^, + 9g'[, g^ g^- 

^9l2 9ll 923 + 18^231 5^12 9l3 923 " 12^231 5^12 9l3 923 + 18^231 9l2 9l3 923 + 
18C23l5'l2 5^13 5^23 ~ 12^23l5'l2 9l3 923 + 9^23lfi'l2 9l3 923 + 9^231 9l2 9l3 923^ 
6^231 fl'12'^ 913 923 - 12^231 5'l2 5'l3 5'23 + 18^231 5'l2 9^3 9%, + 18^231 5'l2 9x3 923 + 

m23i9%9\3 923Mn3)zoir23) ■ (3.51) 

The matrices ^231; depending on the angles formed by the vectors ri, r2 and r^, are 

defined in Eq. (I2.130p . Following Ref. [19j we have considered only the direct term 

of the generalized operatorial correlations. As a consequence, in the previous equation 

p p 

9ij - 9dd,ij- 

In order to find the optimal values of the variational parameters, we have employed 
a procedure similar to simulated annealing, the details of which have been explained in 
Section 12.5.41 of this Thesis. 
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For the calculation of the density dependent potential of Section I3.2[ we constrained 
the difference between the Pandharipande-Bethe (PB) and the Jackson-Feenberg (JF) 
kinetic energies to be less than 10 % of the Fermi Energy Tp and required the sum rule 
involving the scalar two-body distribution function, g'^{ri2), to be fulfilled with a precision 



For this comparative study of chiral inspired potential, in variational calculations of 
SNM we have imposed the further condition, firstly considered in Ref. [12], that the sum 
rule of the isospin component of the two-body distribution function 



be also satisfied to the same accuracy. 

Using also the sum rules for the spin and spin-isospin two body distribution functions 
leads to a sizable increase of the variational energies, which turn out to be much higher 
than those obtained releasing the additional constraints, as well as the AFDMC results. 
The same pattern is observed in the results of variational calculations not including TNF. 

For this reason, we have enforced the fulfillment of the sum rules for g'^{r 12) and g'^{r 12) 
only. 

For potentials other than UIX, it turns out that the variational energies of PNM 
resulting from our optimization procedure are lower than the AFDMC values at p > po- 
By carefully analyzing the contributions of the cluster expansion diagrams, we realized 
that the value of diagram f l3.1[ a) was unnaturally large. In particular, we have found 
that a small change in the variational parameters leads to a huge variation of the value 
of the diagram. Moreover, the minimum of the energy in parameter space was reached in 
a region where the kinetic energy difference was very close to the allowed limit. 

To cure this pathology, we have constrained the PB-JF kinetic energy difference to be 
less than 1 MeV, regardless of density. The variational energies obtained imposing this 
new constraint are always larger than the corresponding AFDMC values and the value of 
diagram (13. lb ) is brought under control. For the sake of consistency, the same constraint 
on the kinetic energies has been also applied to SNM. In addition, the variational energy 
minimum does not correspond to the maximum allowed violation of the constraints. As 
a consequence, it would be largely unaffected by a slight modification of the constraints. 

3.3.3 AFDMC calculations 

We have computed the EoS of PNM using the AFDMC approach with the TM' and 
NNLOL chiral potentials combined with the Argonne fg NN interaction. 

An efficient procedure to perform AFDMC calculations with three-body potentials is 
described in Ref. [154J . Since V3 is equivalent to the anticommutator term of the UIX 
model (while the commutator, V4, is zero in PNM), and in PNM the Ve terms of both 
the TM' and NNLOL potentials do not show any formal difference with respect to the 
repulsive term of UIX, the inclusion of these terms reduces to a replacement of constants 
and radial functions. The authors of Ref. [154J also described how to handle the Vi for 
the TM model, and no further difficulties arise in the case of the NNLOL potential. 



of 3% 




(3.52) 
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As the Vd term has never been encompassed in AFDMC, it is worthwhile showing how 
the calculation of this term reduces to a matrix multiplication. The expectation value of 
Vd is given by 

(Vd) = Yl [^^(' • ^■^) + ■■ tk) + VD{k : ij)] (3.53) 

i<j<k 

with Voii '■ jk) = Voii '■ kj) (otherwise all six permutations need to be summed). Thanks 
to this property one can write 

(Vd) = J2 ^^^3 • • (3-54) 

i<k,j 

It is possible to write Vd^J : ik) of Eq. f l3.45p in terms of cartesian components operators 
VdU : tk) =(F„, l3jZjj-5k + Zai;i3jT.,j-Sk)WT'^j , ^]<^k} ) (3.55) 

where 

Ya^■,f,, = F(r,,)r^ 

T„,;3,=r(r,,)(3f,^4-r^). (3.56) 

The anticommutation relation {erf, cTj} = 25°'^ makes the expectation value of Vd a sum 
of 2)N X 3iV matrix multiplications 

(1/^) = 2^(F,, ■,/3jZi3j-Sk + Zai-iSjY/sj-sk + Tai-pjZpj-sk + Zai;i3jTi3j.Sk)o'iO'k 
i<k,j 

= 2 ^({r, Z} + {T, Z}U,sk cjtai , (3.57) 

i<k 

analogous to those of Ref. |154] . In order to compute the expectation value of Vd the 
former expression has been added to the cartesian matrices associated with the two-body 
potential. 

Finally, as the free-gas value obtained with 66 neutrons turns out to be very close to 
the thermodynamic limit of 73.00 MeV, the finite size corrections for 66 neutrons tend to 
be small. 

For the same reasons adduced for the AFDMC calculations of the density dependent 
potential, we simulated PNM with A = 66 neutrons in a periodic box, using the fixed- 
phase approximation. Finite-size effects are expected to be larger when the density is 
bigger, as the dimension of the box decreases. In order to check the validity of our 
calculations, at p = 0.48 fm~'^ we have repeated the calculation with 114 neutrons. For all 
the potentials, the energies per particle obtained with 114 neutrons are higher than those 
obtained with 66 neutrons. The authors of Ref. [152] found a similar behavior for PNM 
at p < 0.32 fm"^ in the case of the v'g plus UIX hamiltonian, and ascribed part of this 
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difference to the Fermi gas energy, amounting to 72.63 MeV and at 74.15 MeV for 66 and 
114 neutrons, respectively. However, the difference of the energy per particle obtained 
with 66 and 114 neutrons is always within 4 MeV. It is worth noting that once finite-size 
effects on the Fermi gas energy are accounted for, the residual finite-size effects do not 
exceed 4% of the energy per particle. 

3.4 Nuclear Matter EoS 

3.4.1 TM' potential 

The results of Fig. I3.13[ showing the density dependence of the energy per nucleon 
in PNM, indicate that, once the new constraint on the difference between PB and JF 
kinetic energies is imposed, the agreement between FHNC/SOC (solid line) and AFDMC 
(triangles) results is very good. 

The most striking feature of the results displayed in Fig. 13.141 is that, despite the 
parameters of the three body potentials being different, all SNM EoS obtained from the 
TM' potential turn out to be very close to each other. This is probably due to the fact 
that these potentials are designed to reproduce not only the binding energies of and 
^He, but also the n-d doublet scattering length "^a^ui- 

It is remarkable that although the parameters of TM' potentials were not adjusted to 
reproduce nuclear matter properties, the EoS saturates at densities only slightly lower than 
Po = 0.16fm~'^, and the compressibilities are in agreement with the experimental value 
K ^ 240 MeV. On the other hand, the binding energies are larger than the empirical value 
Eq = —16 MeV and rather close to the one obtained from the UIX potential, ~ 10 MeV 
and shown in Section 13.21 The numerical values of all these quantities are listed in Table 

Table 3.5: Saturation density, binding energy per particle and compressibility of SNM 
corresponding to the TM' EoS displayed in Fig. 13.141 





tm; 


TM'2 


TM'3 


Po (fm ^) 


0.12 


0.13 


0.14 


Eo (MeV) 


-9.0 


-8.8 


-9.4 


K (MeV) 


266 


243 


249 
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Figure 3.13: (Color online) Equation of state of PNM obtained using the AFDMC (trian- 
gles) and FHNC/SOC (solid lines) approaches with the TM; (a) and TM'2 (b), TM'3 (c) 
plus Vg hamiltonian. 
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Figure 3.14: (Color online) Equation of state of SNM resulting from FHNC/SOC varia- 
tional calculations with the TM' plus v'g hamiltonian. 



3.4.2 NNLOL chiral potentials 

The results displayed in Fig. 13. 151 show that, as in the case of the TM' potentials, the 
EoS of PNM computed within the AFDMC and FHNC/SOC schemes are very close to 
each other over the entire density range. 

The EoS of Fig. 13.151 are softer than those obtained from both the TM' (compare to 
Fig. I3.13p . and UIX (see Fig. I3.1UP potentials. This is due to the ambiguity in the term 
Ve, discussed in Section [3.3.11 

In the NNLOL2, NNLOL3, and NNLOL4 models the constant ce is negative. There- 
fore, the contribution of Ve is attractive, making the EoS very soft. When Ve is repulsive 
(i.e. Ce is positive), as in the NNLOLi potential, its contribution is very small and the re- 
sulting EoS, while being stiffer than those corresponding to the other NNLOL potentials, 
remains very soft. 

The recent astrophysical data of Ref. [16j suggest that the EoS of PNM be at least 
as stiff as the one obtained with a readjusted version of the effective density-dependent 
potential of Lagaris and Pandharipande in combination with the Argonne v'q two-body 
interaction [31]. Therefore, the EoS resulting from chiral NNLOL potentials are not likely 
to describe a neutron star of mass around 2Mq. 

The SNM EoS corresponding to the NNLOL potentials are displayed in Fig. 13.161 
The fact that the NNLOL4 potential provides the stiffest EoS, while in PNM provided 
the softest, is not surprising. As discussed in Section I3.3.1[ when the contact term is 
attractive in PNM, it is repulsive in SNM, and viceversa. A large cancellation between the 
repulsive core of the Argonne fg and the strong attractive contact term contribution of the 
NNLOL4 potential is observed. This could influence the variational results, which for this 
particular three-body force could be less accurate than for the other interactions. As the 
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corresponding AFDMC results do not show a similar behavior, giving a simple physical 
interpretation to the inflection point at p ~ 0.24 fm~^ resulting from the FHNC/SOC 
calculations turns out to be difficult. 

The results listed in Table 13.61 show that none of the chiral NNLOL potentials fulfills 
the empirical constraints on the SNM EoS. All potentials overestimate the saturation den- 
sity, while the compressibility is compatible with the empirical value only for the NNLOL2 
and NNLOL3 models. As for the binding energies, they are closer to the experimental 
value than those obtained using both the UIX and TM' potentials. 

As a final remark, it has to be noticed that using the scalar repulsive term instead 
of provides more repulsion, resulting a stiffer EoS. As stressed in Section 13.3. 1[ this 
issue needs to be addressed, taking into account all terms that become equivalent in the 
limit of infinite cutoff only. Moreover, since the discrepancies among these terms are of 
the same order as the NNNLO term of the chiral expansion, other contact terms have to 
be included [155]. 




Figure 3.15: (Color online) Same as Fig. 13331 but for NNLOLi (a), for NNLOL2 (b), for 
NNLOL3 (c) and for NNLOL4 (d) plus hamiltonian. 
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Table 3.6: Saturation density, the binding energy per particle, and the compressibility 
related to the NNLOL Eos displayed in Fig. IXTBl 





NNLOLi 


NNLOL2 


NNLOL3 


NNLOL4 


Po (fm-3) 


0.21 


0.20 


0.19 


0.17 


Eo (MeV) 


-15.2 


-14.6 


-14.6 


-12.9 


K (MeV) 


198 


252 


220 


310 




Figure 3.16: (Color online) Same as Fig. I3.14[ but for NNLOL plus fg hamiltonian. 



Chapter 4 

Response of nuclear matter 



The understanding of neutrino-nucleons interactions is required in both astrophysics, 
mainly in the context of supernovae explosions and of the cooling of neutron stars, and high 
energy physics, in particular to reduce the systematic uncertainties of the data analysis. 

In the low-momentum transfer regime (|q| of the order of tens MeV), the non relativis- 
tic limit of the weak current matrix element is expected to be applicable. The nuclear 
response to weak probes delivering energy u and momentum q at leading order in |q|/m 
reads 

^(q, ^) = \Y. I{^/I0q|^o}r5(c^ + Eo- E^) , (4.1) 

/ 

In the above equation A is the particle number and Oq the one-body weak operator 
that induces a transition from the ground state to the excited state I^P/}, which 

are eigenvalues of the nuclear hamiltonian with energies Eq and Ej, respectively (see Eq. 

(E2ZD). 

The non relativistic Fermi (F) and Gamow- Teller (GT) operators describing low energy 
weak interactions are 

i i 

Of = E^rw = ^^E^^'"^^-^' (4-3) 

i i 

where gv = 1-00 and qa = 1-26 are the form factors at zero momentum transfer, while 
is the isospin-rising operator acting on the i-th nucleon. 

The spin longitudinal and spin transverse components of the Gamow- Teller response 
functions are defined as the components parallel to q and orthogonal to q, respectively. 
They can differ significantly at larger values of |q| and, in principle, have to be calculated 
separately. However, when not otherwise specified, with "Gamow- Teller response" we refer 
to the total response, given by the sum over the cartesian components: 

5^^(q,a;)= fl^n^mmu; + E, - E^) , (4.4) 

a=x,y,z f 
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The differences between the integrated spin-longitudinal and spin-transverse responses 
will be discussed in Section 14. 5. 31 

Calculations of the charged current weak response have been recently performed in 
Refs. |17t 1156] . In the high energy domain the short range correlations of the CBF play 
a major role, while for low neutrino energy < lOMeV, the response is mainly affected 
by long range correlations that lead to excitation of collective modes, taken into account 
within the correlated Tamm-Dancoff approximation (CTDA). 

Even including the spin-isospin dependent correlations of Eq. (I2.30p . generally the 
CBF states are not eigenstates of H. However, it has been argued by Feenberg that the 
Hamiltonian matrix has smaller non diagonal elements in CBF than in the non interacting 
FG basis. Hence, perturbative many-body calculations are expected to converge more 
rapidly in CBF. Neglecting orthogonality corrections, the CBF weak operators matrix 
elements read 

{*,i6,i*„}^ (iic:^^^ (4.5) 

Following Refs. |17t 1156) . in this Thesis we only consider transitions between the 
correlated ground-state and correlated Iparticle-lhole (Ip — Ih) excited states. The n > 
2 particle-hole correlated states give a smaller contribution, mainly at large excitation 
energy, the size of which will be estimated at a later stage, studying the sum rules of the 
weak response. 

The CBF matrix element between the ground-state and Ip — Ih excitation reads 

where Pm and hi denote the whole set of quantum numbers of the single nucleon state, 
namely the momentum, the spin and the isospin projections along the axis. 

This quantity, entering all our calculations of the response function, will allow us to 
define the effective weak operators, as discussed in the next Section. 



4.1 Effective weak operators 

We define the effective weak operators O^^^ through the relation 

As for the calculation of the hamiltonian expectation value, a cluster expansion of 
the weak operator correlated matrix element can be performed [157]. The smallness 
parameters in this case are f^j — 1 and ff-. Following ^33j, we rewrite the matrix element 
of Eq. (|12D as 



[^,^;H.\Pd^:F\m,) ^WWW^ -K-Rb. (4.8) 
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A cluster expansion of botli Ra and Rb needs to be performed, so tliat tlie order n of tlie 
perturbative series in terms of ffj — 1 and ffj is given by 

i%^;k.\OM^-'> = ■i?^)<5,+^=.. (4.9) 



Lm 



4.1.1 Cluster expansion of R}y 

It is convenient to carry out the cluster expansion of R^ first, since it does not involves 
the transition operator. The product of correlation operators can be written as in Eq. 
I2.68[ associated with the denominator of the two-body distribution function 

jrtjr= i + ^X(2)(xi,Xj-) + ^ X^^\xuXj,Xk) + ... (4.10) 

i<j i<j<k 



Because of the symmetry of the ground state wave functions, one finds 

{A-N)\N\ 



J2 (^o|X(^)(x„ . . .,xa)\^o) = ^' ,^, (^o|X(^H^i, • • ■,x^)\^o) (4.11) 



i<j<... 



and an analogous relation holds for the Ip — Ih state. 

Using the results of Appendix [B], it is possible to extract particles from the Slater 
determinant of the ground state, to obtain 



.xn)x 



Ail/jnA^l) ■■■ ^nj,{xN)]. (4.12) 

In order to make the simplification among disconnected diagrams manifest, it is worth 
introducing a new index, Ui = 1, . . . , hi — 1, hi + 1, . . . , A, labeling the A — 1 states of 
a system lacking both single-particle states Pm and hi. Thus, adding and removing the 
contribution of the hole state, the A^-body term of the numerator reads 



:vl/o|^t^|vI/o) = $^^ J dx,,.. 



i^nA^l) . . . 'lp*nj,{xN)X^^\xi, . . . ,XN)A['^nAxi) . . . llJnr,{xN)] + 

1p*h^{xi) . . .■^^^{xn)X^^\xi, . . . ,XN)A['^hXxi) . ■ .'^nr.ixN)] + ■■■ + 

ijl^{xi) . . .tpl^{xN)X'^^\xi, . . . ,XN)A[i'ni{xi) . . . iJh{xN)] • (4.13) 

Although the sum of the cluster terms precisely gives the denominator of the two- 
body distribution function, the diagrammatic rules stemming from Eq. fl4.13p slightly 
differ from those of g^{r 12). 

• Wavy lines represent both scalar and operator correlations, f^ - — 1 and ff- and two 
kinds of vertex and exchange lines are present. 
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• The bare vertex, arising from ipniixiyipniixi), does not involve the hole state hj; it 
tends to the "standard" vertex in the thermodynamic limit only. On the other hand, 
the direct term of the hole state, iphi{xi)*iphi{xi), is represented by an hj-veitex: a 
closed loop around the internal point i. 

• The hole state is not exchanged in the bare exchange lines, but there is an additional 
/ij-exchange line, coming from (xj), which connects points i and j. At 
most one hi vertex or hi exchange line can appear in a given diagram. 

It it possible to factorize the sum of cluster diagrams in two subsets, as shown in Fig. 
14.11 The first subset, denoted with 1 + Ylci^i)^ addition to the unity, contains all 
the connected diagrams having one hi vertex or one /ij-exchange line. All the remaining 
diagrams, both connected and disconnected, belong to the second subset. 




Figure 4.1: Factorizations of the sum of cluster diagrams contributing to the numerator 
of Rb- The square brackets on the left enclose 1 + Ylci^i) fro^^ which the denominator, 
1 + YlciPm)^ can be obtained replacing hj with Pm- 

The diagrammatic cluster expansion rules for the denominator can be readily obtained 
from those of the numerator by replacing the hole state hi with the particle state pm- The 
sum of the denominator's cluster diagrams can be factorized in two subsets. The first, 
1 + Ylc(P'm)^ made of the unity and of connected diagrams with one Pm vertex or one 
Pm-exchange line, while the second cancel the corresponding one from the numerator. 

Therefore, the following analytic expression for Rh is obtained 



At first order in ^ — 1, the explicit expression of Eq. (14. lip , diagrammatically repre- 
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sented in Fig. 14. 2[ reads 



Rh 



^ 1 + ^ / dru Ea,(«i«P^ 1^12(1 - A242e*P--^i^)|«i«p™) 



(4.15) 



where \ai) denotes the spin-isospin state of particle i. 



1 + 2x 



1 + 2x 



Pm Pm 

. + 2x<^i> 



Figure 4.2: First order diagrams in terms of ^„/^ — 1 contributing to i?^. 



4.1.2 Expansion of Ra 

Apart form the square root, the denominator of Ra is identical to the numerator of Rf,, 
hence the same diagrammatic rules apply. In the numerator, on the other hand, the weak 
transition operator Oq appears. Using the symmetry of the Ip ~ Ih and ground-state 
wave functions, we can write 

{^,^.,,^T^d^T\^o) = A{^,^.,,^\T^d^il)T\^o) , (4.16) 

with Oq(l), defined in Eq. (14. 3p . acts on particle 1 only. 
The cluster expansion of J-'"'"Oq(l)J-' reads 

.Ft6q(l)^ = 6q(l) + ^xf)(xi;x,)+ J2 x['\xl■,x^,XJ) + ■■■ (4-17) 

Ki l<«<j 

Using again the results of Appendix [HI the orthogonality of the Slater minors intro- 
duced therein and the properties of the antisymmetrization operator A, one obtains 



^(*P™;/.J-^^Oq(l)^|vl/o) = J]— ^ J] / dr,, 

V'p„(^l) • • • i^*nA^N)X[^\xi] ■ ■ .,XN)A[i)hXxi) ■ ■ ■ i^n^ixN)] + • • • + 

V>a^(a;i) . . . %S^N)]xf\x,- XA)A[^nA^i) . . . ^h,M]\ ■ (4.18) 

where the index has been defined just above Eq. fl4.13p . The diagrammatic rules for 
Rfy need to be modified in order to include the weak transition operator acting on particle 
1 and to deal with the p — h exchange line. 
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• The weak transition operator Oq(l), carrying momentum q and a spin-isospin op- 
erator, is attached at point 1. It is represented by a dashed hne and an arrow 
indicating the flow of the momentum q. 

• All internal points but point 1 must be reached by at least one wavy line. 

• Like in the diagrammatic expansion of Rh, bare vertices are due to ip^,{xi)ipn.{xi) 
and are lacking of the hole state; in the bare exchange line the hole state is absent as 
well. If the diagram is finite in the thermodynamic limit, both vertices and exchange 
lines tend to the "standard" ones of the FR cluster expansion. There are neither h- 
and p- vertices, nor h- and p- exchange lines. 

• There is one single p/i-exchange line connecting points i and j, arising from tlj*^{xi)ilJhi{xj), 
or one single ph vertex, coming from ip*^{xi)'4)hi{^i)- In order for a given diagram 

not to vanish, either the p/i-exchange line or the ph-veitex need to be connected 
with particle 1. 

As in the case of Rh, the sum of cluster diagrams can be factorized in two subsets, 
shown in Fig. 14.31 The first, indicated with 'Yl,c^QiPm-,hi), consists of the connected 
diagrams with particle 1 and one p^-exchange line or one ph- vertex. The second one 
includes both connected and disconnected diagrams without p/i-exchange and p/i-exchange 
lines. 




Figure 4.3: Factorization of the diagrams contributing to the numerator of Ra] in the 
square brackets on the left Y^Q^q.Pmhi) is contained. 



As in the calculation of diagrams not containing the weak transition operator, the 
p/i-exchange line and the ph- vertex cancel with the corresponding ones arising from the 
denominator. Hence, the following expression for Ra is obtained 

The first order diagrams in / — 1 contributing to Ra are displayed in Fig. 14.41 
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qA qA 



+ 



1 + 2x 



Figure 4.4: Diagrammatic representation of tlie first order in / — 1 of 



4.1.3 Cluster expansion of the matrix element 

Substituting Eqs. fl4.14p and fl4.19p in Eq. fH^ . tlie CBF matrix element of tlie weak 
transition operator takes tlie form 

(^..;.JO,|0) = .^^^ff.fl , , , (4.20) 

the diagrammatic representation of which can be easily deduced from Figs. 14.11 and 14.31 
For the sake of giving a unified description of matrix elements associated with both the 
Fermi and Gamow- Teller transitions, it is convenient to distinguish the common radial 
part from the specific zero momentum form factors and spin-isospin operators, defining 

Oq(l) = e'^'^-O^^iX) ■ (4.21) 

From Eq (gS]) it follows that O^^(l) = gyr^ and O^'^iX) = QA^irt- 

Zero-th order 

The weak operator's matrix element at zeroth order in / — 1, displayed in Fig. 14. 5[ 
corresponds to the non interacting Fermi gas result 

(^P.;.JOJO)W = A{^>,^,,^d^{l)\^>,) ^ ^ of . (4.22) 

Using Eq. fl4.18p . the right hand side of the above equation can be written in the form 

of = cirie-^P--^e^''^^e^'^^''M«pJa.(l)|a..) 

= 5q,p™-h.(«p„|ar(l)|«h.), (4.23) 



124 



Response of nuclear matter 



where the discretized momentum conservation is expressed by the Kronecker delta func- 
tion. 

The SNM spin-isospin operator matrix elements (ap^|Oo-T(l)|a/i,); along with all the 
other matrix elements appearing in the rest of this Section, are given in Appendix [E] for 
both Fermi and Gamow- Teller transitions. 

First order 

At first order in / — 1, diagrams arising from both the numerator and the denominator 
of Eq. ( I4.20p contribute to the matrix element. In order to include all the first order terms, 
three-body cluster contributions have to be taken into account. These contributions, 
neglected in |157[ W7\ 1156] bring about some inconsistencies, to be discussed later. 

At first order in / — 1, the two-body cluster term of Eq. fl4.17p is given by 

X^'\x^;x2) = {f -1,0^(1)}. (4.24) 

The analytic expressions of the two-body diagrams of Fig. 14. 6[ corresponding to the 
numerator of Eq. (l4.2Up . can be derived by substituting the latter relation in Eq. fl4.17p . 
Since all of these diagrams are converging in the thermodynamic limit, the index can 
be safely replaced by the standard n^, including the hole state hi. Thus 

oSi = ^^q,P,„-h, / rfri2e"''i-'--5^(aia,„|{/-l,a.(l)}|aia,,), (4.25) 
OnI = ^Kpm-h^ J dru^{ap^a2\{fi2 - 1, 0<xr(l)}|ah,tt2) , (4.27) 

0.2 

O^nI = -^'5q,P.-h. 1 rfri2e-^'^-^-£i2 5^(a,^a2|{/i2 - 1, a.(l)}AYI«/..«2) • • (4.28) 

The first order diagrams of the denominator of Eq. fl4.20p are displayed if Fig. 14.71 
By noting that the two-body cluster term of Eq. fl4.10p is given by 

X(=^)(xi,X2) = 2(/i2-l). (4.29) 



I 



Figure 4.5: Zero-th order matrix element O 
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Pmhj p^hi Pmhj p„h 

. . 

qA qA^^ ^ qA 

^Na ^Nb ^Nc '^Nd 



Figure 4.6: Two-body diagrams of tlie first order term in / — 1 coming from tlie numerator 
ofEq. flCTj) . 



Pm Pm hi hi 

2.. 2x^^ 2X. 2x 

OS 0« og> otli 



Figure 4.7: First order diagrams stemming from tlie denominator of Eq. (l4.2Up 
the following analytic expressions of the diagrams of Fig. 14.71 can be derived 

O^Dc = — I rfri2^(aia/ij/i2 - MaiOh,) 

= -V / ^ri2e*"-'^^£i2 5^(«ia/.J(/i2 - • (4-30) 



"1 



Performing a Taylor expansion of the denominator of Eq. ( I4.20p . it can be shown that, 
up to an error of the order of (/ — 1)^, once multiplied by O^"*, and cancel O^^, 
while O^^l and OgJ can be interpreted clS cl correction to Oj^^ and O^^. 

At first order in / — 1 three body cluster diagrams have to be taken into account. 
They are associated with the numerator of Eq. fl4.20p and arise from the the three body 
cluster term of Eq. (14.171) 

X(3)(xi;X2,X3) = {/23 - l,Oq(l)} = 26q(l)(/23 - 1) . (4.31) 
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Because of the property {xi;x2,Xs) = {xi;x3,X2), when the latter expression for 
the three-body cluster term is inserted in Eq. (14. 181) . one finds 

^ / dxidx2dx3 2?/';^(xi)?/';2(x2)V'p„(a;3)C'q(l)(/23 - l)-A[^/'ni(a;i)?/'^^2(x2)^/'h,(a;3)] + 



(4.32) 



i^*p^{XlWn2{x2)lp*n,{x3)d^{l){f23 " (xO^/-^^ (a;2)^n3 (a^s)] 



The first order three-body cluster diagrams can be grouped in two subgroups; those 
coming from the first line of Eq. fl4.32p are depicted in Fig. 14. 8[ the others in Fig. 14.91 



Pmhj 




o 



(1) 

Ne 




Pmhj 




Pmh 




Pmh 



Nh 





Pmhj 



Figure 4.8: Three body diagrams emerging from the first line of Eq. fl4.32p . 



Diagrams and vanish, as the part of the diagram containing the weak tran- 
sition operator is connected neither to the ph-vertex nor the the p/i-exchange line. In the 
thermodynamic limit diagram 0^^\ reads 



O 



(1) 

Ng 



2p2 



(aia2ap„|0<xr(l)(/23 - 1) A'2^l"i"2ahJ • 



(4.33) 



Note that in the limit of zero momentum transfer, q — )■ 0, cancels the contribution 

of diagram O^^. This is an indication that three-body diagrams need to be taken into 
account. 
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Figure 4.9: Three body diagrams emerging from the second hne of Eq. fl4.32p . 



For the diagram 0^]^ we obtain 

)(!)_ V 



on 

{aia2apjdar{l){f23 - l)Pi7l"i"2ah,) ■ (4.34) 
The latter expression vanishes, because 

J (iri3^i3e^P™-"" = . (4.35) 

For the same reason we find 

(«l«2«p„|ar(l)(/23 - 1)A7A7I«1«2«/..) = . (4.36) 

The last diagram coming from the first line of Eq. f l4.32p gives a non zero contribution 

(«l«2«p„|ar(l)(/23 - 1)A7A7I«1"2«/^.) • (4-37) 
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The first two diagrams of Fig. 14. 9[ being disconnected, simplify with the denominator. 
All the remaining diagrams of Fig. 14.91 vanish. This can be shown by noting that the bare 
exchange line is represented by the modified Slater function that does not contain the 
hole state 

J rfr,,-4e*'--- = , (4.38) 

and using the matrix elements 

(ap„a2a3|0^r(l)(/23 

^ (ap„a2a3|O^T(l)(/23 

^ («p„«2a3|0^r(l)(/23 
a2037^0hi 

Second order 

As in the first order case, diagrams coming from both numerator and denominator of 
Eq. (I4.20p have to be taken into account. In order to consistently include all the second 
order terms, in principle up to five-body cluster diagrams need to be computed. In our 
calculations |158] . we have considered second order two-body cluster terms only, as in 
Refs. |157t HZl 1156] . Second order three-body cluster diagrams are in fact much smaller 
than the corresponding two-body ones, as confirmed by a direct computation of some of 
them. 

Considering the second order term, to the two-body cluster of Eq. fl4.24p becomes 

X^^\x,;x2) = {/i2 - l,Oq(l)} + (/i2 - l)0q(l)(/i2 - 1) . (4.40) 

The additional term (/12 — l)0q(l)(/i2 — 1) brings about new diagrams for the nu- 
merator of Eq. fl4.20p which are analogous to those of Fig. 14.61 As can be seen in Fig. 
I4.10[ the only difference consists in the replacement of the single wavy lines with doubly 
wave lines, accounting for the second order correlations (/12 — l)(/i2 — !)• 

The second order two-body numerator diagrams of Fig. 14.101 are given by 

02 =^'^q,P.-h. / dme-'^^^- 5^(«i«p„|(/i2 - l)a.(l)(/i2 - l)|«i«h.) , (4.41) 

O^Ji = - ^ W.-h. 1 rfri2e^P'"-^-£i2X 

5^(«i«p,J(/i2 - l)a.(l)(/i2 - l)A7l"i«/..) , (4.42) 

Ol 



- l)Pi7l«h,«2«3) = 
-l)P,7A7l«/^.«2«3) =0 

- l)Pr2"A7l«/..«2«3) = . (4.39) 



4.1 Effective weak operators 



129 



Pmhj 



o 



(2) 
Na 




Pmhi 



o 



(2) 




Figure 4.10: Two-body diagrams of the second order term in / — 1 coming from the 
numerator of Eq. fl4.20p . 



O 



(2) _ P 



Kpm-h, j rfri2^(ap„a2|(/i2 - l)0^r(l)(/i2 - l)|ah,a2) , (4.43) 



O 



(2) 
Nd 



q,Pm-hi 



dri2e 



5^(ap„a2|(/i2 - l)ar(l)(/i2 - l)AY|a/..a2) 



(4.44) 



0^2 



As far as the denominator of Eq. fl4.20p is concerned, the two-body second order 
diagrams, shown in Fig. I4.20[ can be obtained from those of Fig. (14. 7p by again replacing 
the single wavy lines with doubly wavy lines. As a matter of fact, once the second order 
term in considered, the two-body cluster term of Eq. (14.101) reads 



X(2)(a;i,X2)=2(/i2-l) + (/i2-l)' 



(4.45) 
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Figure 4.11: Diagrammatic representation of the second order two-body cluster term of 
the denominator of Eq. (I4.20p . 



The analytic expressions of the second order two-body diagrams of Fig. 14.111 are very 
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similar to those of Eq. f l4.30p 

Ogl = -^j drue'^--^^Hu $^(aiap.|(/i2 - 1)^A7I«i«p„.) • (4-46) 



The results of the second order two-body cluster terms coincide with those obtained 
by the authors of Refs. [TFTl E]. 



4.2 Effective interaction 



Using the formalism of CBF and the cluster expansion technique, the authors of Ref. 
|157[ HZ], were able to develop an effective interaction, obtained from the bare Argonne 
Vg potential, which incorporates the effects of the short-range correlations. In Ref. |37] . 
the two-body effective interaction of Ref. |157[ H7] was improved with the inclusion of 
the purely phenomenological density dependent potential of Ref. [lllj, accounting for the 
effects of interactions involving more than two nucleons. The CBF effective interaction, 
vll^ , suitable for use in Hartree-Fock calculations, is defined through the matrix elements 
of the hamiltonian in the correlated ground-state 

{^o\PhP\^o) =Tf + (^ol^^il^l^o) • (4.47) 



As suggested by the above equation, the effective interaction allows one to calculate 
any nuclear matter observables using perturbation theory in the orthonormal FG basis. 
However, in general, extracting the effective interaction is a very challenging task, in- 
volving difficulties even more severe than those associated with the calculation of the 
expectation value of the hamiltonian in the correlated ground state. 

The procedure developed in Ref. [H] consists in carrying out a cluster expansion of 
the Ihs of Eq. f l4.47p and keeping only the two-body cluster contribution. The sum of the 
two-body cluster contribution of the potential and kinetic energies of Eqs. f l2.149p and 
f l2J[53|) reads 



^ I rfri2CTr 



2b 



F12V12F12 (ViFi2)(ViFi2; 

m 



;i - A24) 



(4.48) 
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On the other hand, the expectation value of the effective potential is given by 



= ^ / rfrisCTr 
2b 2 ' 



hi 



(4.49) 



Therefore, the effective potential at two-body cluster level turns out to be 



hi 



FvuF 



2b 



— (ViFi2)(ViFi2) 
m 



(4.50) 



The effective potential of the above equation slightly differs from the one reported 
in the literature. The authors of Refs. [171 1156] have not done the integration by parts 
leading to Eq. f l2.153p . As a consequence they have neglected the terms in which the 
gradient operates on both the correlation function and on the plane waves. In our effective 
potential, these terms, although small compared to the other contributions, are fully taken 
into account. 

We have improved the effective potential by adding the tree-body cluster contributions 
of Eqs. (12. 164112. 173p . This allowed us to consistently include the UIX potential, whose 
leading order terms emerge at three-body cluster level. 

As in the construction of the density dependent potential from UIX, the issue of the 
exchange pattern has to be carefully analyzed. The distinctive feature of the present 
calculation is that v\{^ contains the correlation between particles 1 and 2 making possible 
to implement the inversion of P"^ of Eq. (I3.20p in a straightforward way. To be definite, 
consider the three-body cluster contribution of the ground-state expectation value of the 
two-body potential 



(^12) 



3b Z 



1 _ fy^T p'^ 

-'12 ^12 



9 p'^T p'i 
'^-'13 ^13 



2i^l7^l7^12^13^23) 



(4.51) 



Note that the above equation summarizes the terms corresponding to Eqs. (I2.16ip . 
(I2.165p . (I2.166P and (I2.167p . A comparison with Eq. (I4.49P immediately leads to 



^12 



3b 



p / (iraCTrs X{xi, X2; X3){l 



-'12 ^'12) 



P12 



^2 
12 



12 -'13 ^12^13^23 



X 



(4.52) 



A similar argument holds for the three-body cluster term of the kinetic energy and the 
three-body potential contributions to the effective potential. 

In Fig. I4.12l the central and tr components of the effective potentials at two-body and 
three-body cluster level are compared. The starting bare NN interaction is the Argonne 
fg; for the three-body cluster results the UIX three-body potential has been included in 
the calculations. 

Starting from a bare hamiltonian whose only interaction is the Argonne fg NN po- 
tential we have computed the EoS of SNM for the low-density regime using both the 
new three-body cluster effective interaction and the older one with only two-body cluster 
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Figure 4.12: Central (a) and the tr (b) channels of the effective potentials at two-body 
and three-body cluster level calculated for SNM at p = 0.16 fm~^. 



diagrams. The results have been compared with the corresponding FHNC/SOC calcula- 
tions, displayed as a shaded region Fig. l4.13l to account for the PB and JF kinetic energy 

difference. The curve corresponding to vll"^ is much closer to the FHNC/SOC results 



than the one obtained with the older v^^^ 



3b 



3b 



In the lower panel of Fig. I4.13[ the EoS of SNM are shown for an hamiltonian con- 
taining the Argonne fg NN potential along with the UIX three-body interaction model. 
Again the curve obtained from the three-body cluster effective potential is close to the 
full calculation and, it exhibits the saturation, which is a remarkable feature. 
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Figure 4.13: EoS of SNM in the low density regime. In tlie upper panel the bare hamilto- 
nian only contains the Argonne Vq potential, while in the lower panel the UIX three-body 
interaction model is added to it. The dotted and the solid lines display the three-body 
and the two-body cluster effective potential results, respectively. The FHNC/SOC calcu- 
lations are represented by the shaded region, accounting for the PB and JF kinetic energy 
difference. 



4.3 Correlated Fermi gas, Hartree-Fock approximations 

In both the correlated Fermi gas (CFG) and correlated Hartree-Fock (CHF) approxi- 
mations, the weak response of cold SNM, defined in Eq. (14. ip . is given by |17] 

S^g{ci,uj) = jY1 \{^Pr.-M\dl"\'^o)\'S{u + e,^ - e,J . (4.53) 

Pm.hi 
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Within the CGF approximation, the single particle energies are obtained from Eq. (12. 5p 
in the case of a non-interacting hamiltonian 




(4.54) 



The single particle approximation is retained in the CHF; however the potential enter 
the calculation of the single particle energies. As mentioned in Section 12.11 the Hartree- 
Fock approximation is not suitable for nuclear potentials having a repulsive core, like the 
Argonne models, because it does not encompass the correlations between nucleons. We 
use instead the effective potential described in the previous Section, which is appropriate 
for mean field calculations. Thus, Eq. (12. 5 p becomes 

^n, = -^+Yl I ^^i^n,(a;i)V^n,(a;iX/^[V^n.(a;i)^„,(Xj) - i^n.iXi)^)^^^ j)] . (4.55) 



4.4 Correlated Tamm Dancoff approximation (CTDA) 

The Tamm-Dancoff approximation amounts to expanding the final state of Eq. (14. ip 
in series of Ip — l/i excitations 

mi 

Because the hamiltonian is translationally invariant, the total momentum q of the 
state is conserved, and the momenta of the particle, p^, and the hole, hj, satisfy the 
relation p^, — hj = q. 

The eigenvalue equation for the effective hamiltonian defines the excitation energy Uf 

^e//|^^^™^ = ( 5^ + J] f)f ) |v,,}™^ = ^E, + . (4.57) 

i i<j 

Multiplying from the left the previous equation by (\E'p^^./j | and using the orthonormality 
of the Ip — Ih states yields 

Y,{^,„,,jH^ff\^,^.,,^)CL, = Cl^iEo + Uf) . (4.58) 

mi 

Hence, the coefficients Cl^^ can be interpreted as the eigenvectors of the hamiltonian 
Hnj-mi = {"^ p„;hj\H\'^ p^-hi) , while the excitation energies Uf are the associated eigenvalues. 

In Appendix [B] it is shown that singling out particle 1 from the Ip — Ih wave function 
^Pm/ii yields the result 

^p^k. = 4f E(-ir^v^.i(^i)^rA,(^2 ... xa). (4.59) 
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The indexes rii and hi denote the single particle states of the Fermi gas and of the Ip— Ih 
excited state, respectively. Note that, in the above equation, hi must be regarded as a 
function of rii] in particular it turns out that hi{ni) = rii except for hi{hi) = p^. 

The minor "^^J^^^ has been obtained from "^p^hi by removing the column corresponding 
to the state hi and the first row, associated with particle 1. 

Extracting particle 2 from the Slater determinants leads to 

%..k. = }^ _ Yl (-l)"^^"^^'^[^^i(^i)^'^.(^2)]*ri?„..(^3 ...xa), (4.60) 

V V / n.i<n2 

where "^^^^^ the minor of is lacking particle 2 and the state 

Being a one-body operator, the kinetic energy, on account of Eq. fl4.59p . contributes 
only to the diagonal part of Hnj^mi 



A(^p^;hi\- ^^Pm\hi) = / dxxi\)l^[xx)y- ^^i)n^\ 



Xl 



^ + (4.61) 
2m 2m 2m 

|k|<A:F 



The diagonal matrix elements of the effective potential can be computed employing 
Eq. IKm\i 

2 {Pm'i hi\Vi2^ \pm] hi) 

= \Yj d^l,2i'lSxi)^l^{x2)vli^A[i)nAxi)i'n2{x2)] 
rii 

= Y dxi,2'ipl^{xi)iJn2i^2)vll^A[ipm{xi)ip„2ix2)] 

nin2 

+ Y ^^h2%JXi)^*^^{x2)vf2^A[lljpJXi)^n2{x2)] 
ni 

~Y dXi,2^Pn,iXl)ipl^{x2)vll^A['llJnA^l)^hAx2)] 



ni 



dxi,2'ippJxi)i^l^{x2)vf2^A['ipp^{xi)'^hXx2)] ■ (4.62) 

The sum of the first term of Eq. fl4.6ip and the first of Eq. fl4.62p corresponds to the 
ground state energy of the system 



— + -2^ / dXi,2i^*^^{Xl)tp*^^{x2)vti^ A[4jn,{Xl)4jn2{x2)] ■ (4.63) 



|k|<fci? nin2 
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By recognizing in the sum of Eq. f l4.6ip and the Eq. fl4.62p the single particle energy of 
Eq. (I2.5p one gets for the diagonal part of the Ip — Ih matrix elements of the hamiltonian 

i'^Pm-MlH^^l'^p^-M) = Eo + ep„ - eh, + {pm hilvll^ \h, p^) . (4.64) 

where we denote with (p„ hi\0i2\hj pm) the two-body matrix element of the operator Ou 

{Pnhi\di2\hjPm) = / dXi,2'lp*p^iXi)lpl.{x2)di2A[lphjiXi)lpp^{x2)]. (4.65) 



The off-diagonal matrix elements of the hamiltonian come from the two-body potential; 
making again use of Eq. f l4.60p it turns out that 

^^^^^(*P„;h>i2^l*Pm;h,) = {Pnh,\vli^h,p,^) . (4.66) 

Collecting Eqs. (I4.64p and f l4.66p . one finds the following compact expression [17] for 
the matrix elements of the hamiltonian 

{^p,^.h^\H'f^^p^,h,) = {Eo + ep„ - ehJS^Jij + (p„ hi\vll^\hj p^) . (4.67) 

Finally, substitution in the eigenvalue equation fl4.58p gives 

[(^P™ ~ ^hJ^mnSij + {Pn hi\vl{^ \hj p^) C^- = Uf ■ (4.68) 

mi 

The nuclear hamiltonian commutes with the total isospin, T with the total isospin 
projection along z, and with the total spin, S*, but not with the total spin projection 
along Sz because of the tensor term of the potential. Solving Eq. fl4.58p would then 
lead to coefficients Cl^^ such that the final states of Eq. fl4.56p are eigenstates of S*, T and 
T,: I/) ^ \f)TTzS- 

The combinations of particle hole pairs that are eigenstates of the total spin S and its 
projection along the z-axis Sz-, that define the particle-hole Clebsch-Gordan coefficients, 
are shown in Table UT] The differences between the total spin states of the particle particle 
pairs, also given in Table 14. H are due to the phase factor appearing in the canonical 
transformations to particles and holes [52]. For a more detailed discussion, see Appendix 
lEl The treatment of the total isospin can be done in complete analogy, replacing the up 
and the down single particle spin states with the proton and the neutron isospin states, 
respectively. 

It is possible to reduce the computational complexity of the eigenvalue equation by 
carrying out the summation over the spin-isospin projections along the ^-axis of Eq. 
f l4.56p . grouping the combinations of \^p^-h,) with deffnite T, Tz, S and Sz 

|^'/)tt.5= J2 C;jhf''l^p™;h.)TT.5S., (4.69) 
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Table 4.1: Spin configurations for a particle particle pair and a particle hole pair for 
spin-1/2 particles. 



Total spin state particle particle particle hole 



s = 




= 1 


tt 




-ti 


s = 


1,5. 


= 




it) 


;^(tt-ii) 


s = 


1,5. 


= -1 


u 




;t 


s = 


0,5, 


= 




it) 


75(tt + ii) 



A further simplification arises by noting that the final states of both the Fermi and the 
Gamow- Teller transitions are characterized by having T = 1 and = 1. To simplify the 
notation, the isospin indexes may then be omitted 

1^/)^= E C'^Sl*P.;h.)55., (4.70) 

and it is understood that T = 1 and = 1. 

Unlike in Eq. fl4.56p the sum is now restricted to the momentum of the particle and 

f TT S S f 

of the hole and to Sz, furthermore the coefficients C*p^h^ ^^'^ ^Li related through 
the Clebsch-Gordan coefficients of Table 14.11 

The eigenvalue equations for the coefficient C*p^^^ can be obtained from Eq. fl4.58p 

ssA^r>.-MM%^-M,)ssi <S = ClX'iEo + ujf) . (4.71) 

Cf ss 
i, amounts in diagonalizing the block diagonal hamilto- 

nian for the two subsets of the Ip — Ih basis having T = 1, = 1 corresponding to 5 = 

and to 5 = 1 . This is a much less expensive computational task than diagonalizing the 

hamiltonian in the full Ip — Ih basis, as it was in Eq. fl4.58p . As a consequence, this 

approach allows for considering a larger number of momentum states. 

Since the single particle energies and of Eq. f l4.64p depend neither on the spin 

nor on the isospin of the states pm and hi, Eq. (I4.68P can be recasted in the following 

form 

E [(^p-^ - ^hj5p.p„5h.h/5.5^ + (Pn h,\viif\h, pjssis] ciS: = ciT; ■ (4-72) 

In the matrix element of the potential {Pn'^i\'Vi2'^\hjPm)ss'^s^ the spin and isospin pro- 
jections along the z axis of the particle hole pairs p^ — hi and p„ — hj are combined as 
in Table 14.11 to have definite S, and total spin projections along z equal to S'^ and 5*,, 
respectively. We recall that the total isospin and its ^-projection are T = 1 and = 1. 
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The direct and exchange terms for the case S = 0, relevant to the Fermi transition, 
for SNM read 



■J Hm/oOO 



2V 



= / dri2e-''»"'i^[ (p t ^ t \vt2^\n t P t) + (p t ^ 4 1-^12-^1^ t P i) 



eff\ 



4 
V 



(4.73) 



(p„hj|^;i2|hjp, 
1 

~ 2V 



to/ 000 



dYi2e'^'^'^^^[{p t n t \vl2P\2\n t P t) + (P t ^ i \vl2P12\n t P 



+ (pint 1-012-^^21^1 ipt) + {pini\vli^P^2\nip\)] 
I dri2e*'^--^-(^;j2 - ^^2 + 3<2 - S^^iJ ) , 



(4.74) 



where = hj — h^. For the sake of simphcity, the superscript "e//" has been omitted 
where the channels of the effective potential are specified. 

For the Gamow Teller transition, the final state has S* = 1; hence it is necessary 
to compute the nine matrix elements (Pn h^l 1)^2 '^|hj p„) 55^5^ corresponding to Sz,S'^ ~ 
-1,0,1 
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= I y" drise-^^"-- 



2<J + <2(l 



'-'^12 



(4.75) 



(Pn hi|'0i2'''|hj Pm)5=l5^,=l5^ = l 



V 
1 

V 



Vi2 - 



<2 + <2 +(^^i2-<2) 



(4.76) 



12' 



_ 1 

6^2 



[ dr,2e— ,^r^ (^i2-iyi2).i2 

J n2 



(4.77) 



(Pn hi|'Oi('^|hj Pm)5=15^=05^=l 
1 



V2V 



rfri2e^^-"--[ (p t t |^il^A2|n i P t) - (p t 4 1^12'^ ^2!^ 4 p 4)] 



(4.78) 
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(Pn hi|'0i2''^|hj Pm)s=lS'=-lSz=l 



V 



(iri2e-''i"'i2 (p t ^ t l^i^-^l^ ipi) 



'^12" 



^3^12 - ^^12)' 

^2 
'12 



(4.79) 



(Pn hi|'0i2'^|hj p^)|^i5,^_]^5.^^i 

^ y (iri2e*'''^""(p t n t \vl2 Pi2\n ip\.) 

Jkij-ri2 („,t 



V 



(4.80) 



12 
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(4.81) 



= ^ j rfri2e^'^-'-^^[ (j9 t n t |^)il^A2|r2 1 P t) - (j' t ^ i \vl2 Pi2\n \ p \) 



-(pint |t'l'2-^A2|^ i P t) + (P i ^ i |t'l2-^A2|^ i P i)] 
1 



V 



2 ^ 

ikij-ri2 



vl^ - vl^ - vl^ + - 2(t;^, - v%) (1-3 



-'12 



-'12 ^12^ 



^12 



12 



(4.82) 



Replacing the final state of Eq. (I4.70p in the definition of the response, Eq. (14. ip 
yields 



f s 



(4.83) 



/ S Pmhi5z 



4.5 Numerical calculation of the response 

We model the infinite system using a cubic box of side L with periodic boundary 
conditions. Hence, the single particle wave functions are the plane waves of Eq. (12. 9p 
with the discrete momenta of Eq. (12.100 . For zero temperature SNM, all single-particle 
states with |kj| < kp are occupied in the ground state. The momenta of the Ip — Ih 
excitations are such that \hj\ < kp and |pi = hj + q| > kp- For the hole and particle 
momentum to be on the lattice of allowed momentum states in the box, the momentum 
transfer must be such that 

q= ^K.,^.„,^gJ- ng^ = 0,±l,±2,... . (4.84) 

Given the magnitude and the direction of the momentum transfer, the side of the box 
can be determined by inverting the latter equation 

L = ^J<+<+<- (4-85) 



4.5 Numerical calculation of the response 



141 



The size of the basis, that can be increased by increasing ^jn^^ + ri^^ + , has been 



1y 



determined requiring that the response of a system of noninteracting nucleons computed 
on the lattice agreed with the analytical result of the FG model [17 tll56) . The FG response 
is obtained replacing the effective operator with the bare operator in Eq. (I4.53P and using 
the single particle energy of Eq. (I4.54p . The analytical calculations can be performed using 
a continuum of momentum states ^52j, while the numerical result consists of collection of 
discrete delta function peaked at the values of the single particle energies. For a better 
representation of the results, as well as for fitting purposes, a gaussian representation of 
the energy conserving delta function has been adopted 



5{x — xq) — 7- 



o■^/7^ 



exp 



a 



(4.86) 



For sufficiently small values of the gaussian width a, the results become insensitive to it. 

All the results that will be shown in this section refer to SNM at equilibrium density 
p = 0.16 fm~^. 




Figure 4.14: Fermi response functions calculated at g = 0.3 fm~^. The left panel shows the 
response calculated at two-body cluster level for both O^-^-^ and vll^ for different choices 
of correlation functions. In the right panel the three body cluster has been included. The 
plus marks show the response for a non interacting FG. 



4.5.1 CFG and CHF 

FHNC/SOC calculations and their associated minimization procedure, explained in 
Section I2.5.4[ provide a set of correlations function, corresponding to the minimum of 
the hamiltonian expectation value. We have found the best correlation functions for 
the Argonne fg, fg two-body potentials, and for comparison we have also considered 
the correlations of Ref. corresponding to Argonne fig. With these correlations, 

the Fermi and Gamow- Teller response functions have been evaluated in CFG and CHF 
approximations. 
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When only two-body cluster diagrams are considered, as in Refs. 1156] . the CFG 
response, suppressed by a 20 — 25% with respect to the FG case, exhibits a sizable de- 
pendence on the choice of correlations, as shown in the left panels of Figs. 14.141 and 14.151 
These figures refer to a transfer momentum q = |q|(4x + Ay + Az)/\/48 with |q| = 0.3 
fm~^. The folding gaussian function has a width of 0.25 MeV. 

This unphysical effect is removed once the effective weak transition operator is com- 
puted at three-body cluster level. As a matter of fact, the CFG curves in the right panels 
of the aforementioned figures are very close, when not superimposed, to each other. There- 
fore our results appear to be more robust than those of Refs. |47[ [156] , as the physical 
quantities should not be sensitive to the details of the short range behavior of the corre- 
lation functions. 




Figure 4.15: Same as Fig. 14. 141 but for Gamow- Teller response functions. 

A small dependence on the choice of the correlation is observed in CHF calculations, 
also shown in Figs. 14.141 and 14.151 The reason for this lies in the single particle energies, 
which depend on the first four components of the NN potential (see Eq. (12. lip ). It turns 
out that the first four component of the fg potential are very similar to those of fg, while 
the same behavior is not observed for the full fig. 

We observe that the shift of the strength to higher u is slightly enhanced by three-body 
clusters, expecially for Fermi transition. 

4.5.2 CTDA results 

The nuclear matter response calculated in CTDA for |q| = 0.30 MeV is displayed 
in Figs. 14.161 and 14.171 for Fermi and Gamow- Teller transitions, respectively. The peak 
corresponding to the collective mode is shifted to lower energies when the tree-body cluster 
is included. This effect, ascribed to the change of the single particle energies, is mitigated 
when the UIX potential is included in the hamiltonian. 

When the three-body cluster is included, it produces in a depletion of the Gamow- 
Teller resonance at |q| = 0.30, particularly apparent when the nuclear hamiltonian is 
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Figure 4.16: Fermi response functions calculated at g = 0.3 fm~^ using vq' and Vq/ + UIX 
interaction models. The upper panel shows the full response across all the values of cu, 
the lower is a magnification of the small u region. The plus marks and the cross marks 
display the response for FG and CFG respectively. The remaining lines refers to CHF 
approximation and CTDA (see text for details). Responses are folded with a Gaussian of 
width 0.25 MeV. 



lacking of the three-body potential. This feature can be explained by looking at Fig. 
I4.18[ where the Fermi and Gamow- Teller responses are plotted for different values of |q| 
ranging from 0.10 fm~^ to 0.50 fm^^. The maximum of the peak of the Fermi response is 
observed at |q| ~ 0.40 fm^^, like in the two-body cluster approximation of Ref. [37]. For 
the Gamow- Teller case, on the other hand, the position of the resonance is shifted from 
|q| ~ 0.30 fm~^, corresponding to the two-body cluster calculation, to |q| ~ 0.25 fm~^. 




Figure 4.17: The same as Fig. 14. 161 but for Gamow- Teller response functions. 
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Figure 4.18: Fermi (upper panel) Gamow- Teller (lower panel) response functions calcu- 
lated at g = 0.10,0.15,0.20,0.25,0.30,0.40,0.50 fm"! using vq, + UIX potential and 
correlations. 

4.5.3 Sum rules 

The set of final states in Eq. ( 14. ip is not exhausted by Ip — Ih excitations. In 
principle, transitions to more complex multi p — h states should also be considered. So 
far, the contribution of these states have been neglected; however, an estimate of their 
importance can be obtained computing the sum rules. The static structure function is 
defined by 

S{q) = J dcoS{q,co) 

= j{^o\dldM- (4.87) 

While a direct integration of the CTDA response functions allows for the evaluation of 
S{q), from the last line of the latter equation it is clear that S{q) can be computed by using 
the variational ground state (VGS) resulting from FHNC-SOC calculations. In particular, 
the knowledge of the two-body operatorial distribution functions of Eq. (I2.11ip is needed 
to compute the structure function. 

While the VGS calculations include all the multi p — h excitations in the CBF basis, 
in the CTDA only the correlated Ip — Ih states are taken into account. Therefore, multi 
p — h contributions can be estimated from the difference S^'-^^{q) — S'-''^^^{q). Note, 
however, that an interplay between many-body correlations and multi p — h excitation 
could in principle take place. In fact, while VGS includes many-body correlations through 
the chain summations, in the CTDA of Ref. [37] only two-body cluster terms have been 
considered. 
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Figure 4.19: Static response function for Fermi (upper panel) and Gamow- Teller (lower 
panel) transitions. Results are shown for a noninteracting FG (dashed lines), for the 
integral of the CTDA S'(q, u) at two-body (squared) and three-body (stars) cluster level 
and for the VGS (solid lines). 



The static structure functions for the Fermi and Gamow- Teller transitions are dis- 
played in the left and in right panels of Fig. I4.19[ respectively. The dashed lines refer 
to the noninteracting FG, while the squares and the stars represent the two-body and 
three-body cluster results, respectively. The NN potential of the effective hamiltonian is 
the Argonne v'q, while in the tree-body cluster calculation the UIX interaction has been 
used. Note that, in the sum rules' calculation, the authors of Ref. [47j have set to one 
the factor qa of the Gamow- Teller transition. For a better comparison with their results, 
we did the same choice of normalization. 

As noted in Ref. |17], the variational calculations of Ref. |159j . because of the ap- 



proximations involved in FHNC/SOC calculations, do not show the correct behavior, 
S{q = 0) = 0, required by baryon number conservation. On the other hand, the static 
structure function obtained within CTDA does exhibit the appropriate low-momentum 
limit. 

As far as the multi p — h excitations are concerned, the two-body cluster results show 
that their contribution is smaller than the dominant Ip — Ih excitation, but it is not 
negligible. When three-body cluster is accounted for, the VGS and the CTDA curves get 
closer. The shift turns out to be detectable for the Fermi transition case, while for the 
total Gamow- Teller response is very small. This is a clear indication that the difference 
between FHNC/SOC and correlated Tamm-Dancoff results has largely to be ascribed to 
the multi p — h excitations. 

There are experimental and theoretical indications that the spin longitudinal and spin 
transverse response functions can differ significantly due to tensor forces. Thus, we studied 
how the UIX three-body force affects these quantities computing the spin longitudinal and 
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spin transverse static structure functions, defined as 

S\D = \ I dujY,\{^f\q-df\^o}?5{u + E,-Er,) (4.88) 

S^il) = J I rfc^$^|{%|gAdJ^|^o}r5(c^ + i?o-i^n) (4.89) 

(4.90) 

As shown in Fig. I4.20[ the inclusion of the UIX potential brings the CTDA curves 
for S^{q) slightly closer to those of VGS across all the values of |q|. As far as the 
transverse static response function is concerned, the position of the maximum of the 
CTDA calculations including UIX potential almost coincide with the VGS results. At 
small momentum transfer however, the three-body cluster contributions move away the 
S'"^(g) obtained from CTDA from the variational results. 

1.5 I , , , , , , n 1.5 I , , , , , , n 
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Figure 4.20: Spin longitudinal (left panel) and spin transverse (right panel) static response 
functions, calculated using VGS (solid lines), two-body CTDA (squares), two-body CTDA 
(stars), and noninteracting FG (dotted lines). 

Recently |160] . the sum rules relations (additional sum rules with an increasing power 
of u in the integrand can be defined) has been inverted to obtain the spin-response function 
of PNM at zero momentum transfer. The authors of Ref. |160] used the formalism 
of AFDMC to compute the operatorial distribution functions, obtaining very promising 
results. As a follow up of the present work, we are planning to compute the response 
function of PNM and compare it with their findings. 



Conclusions 



The main body of this Thesis has been devoted to the discussion of a novel approach, 
developed in Ref. f27\ fIE\ . allowing one to obtain an effective density-dependent NN 
potential taking into account the effects of three- and many- nucleon interactions. 

The resulting effective potential can be easily used in calculations of nuclear properties 
within many-body approaches based on phenomenological hamiltonians, including the 
effects of strong NN correlations, which can not be treated in standard perturbation 
theory in the Fermi gas basis. Moreover, the derivation of the density-dependent NN 
potential is fully consistent with the treatment of correlations underlying the FHNC and 
AFDMC approaches. 

While the reduction of n-hodj potentials to a two-body density-dependent potential is 
reminiscent of the approach of Refs. [2^ [5U] . our scheme significantly improves upon the 
TNI model, in that i) it is based on a microscopic model of the three nucleon interaction 
providing a quantitative description of the properties of few nucleon systems and ii) allows 
for a consistent inclusion of higher order terms in the density expansion, associated with 
four- and more-nucleon interactions. 

As shown in Chapter [21 the results of calculations of the PNM and SNM equation 
of state carried out using the density-dependent potential turn out to be very close to 
those obtained with the UIX three-body potential. In this context, a critical role is 
played by the treatment of both dynamical and statistical correlations, whose inclusion 
brings the expectation value of the effective potential into agreement with that of the 
UIX potential (see Fig. 13. 8p . This is a distinctive feature of our approach, as compared 
to different reduction schemes based on effective interactions, suitable for use in standard 
perturbation theory [441 1161] . 

Using the density- dependent potential we have been able to carry out, for the first 
time, a AFDMC calculation of the equation of state of SNM consistently including the 
effects of three nucleon forces. The results of this calculation show that the fg-|-UIX 
hamiltonian, or equivalently the one including the effective potential, fails to reproduce 
the empirical data. 

The FHNC results obtained using the v'g potential indicate that the 5-6 MeV un- 
derbinding at equilibrium density can not be accounted for replacing the f g with a more 
refined model, such as fig. Hence, the discrepancy has to be ascribed either to deficiencies 
of the UIX model or to the effect of interactions involving more than three nucleons. 

In order to improve on the UIX three-nucleon interaction, we have performed nuclear 
matter calculations using the new generation of chiral inspired three-nucleon potentials 



148 



Conclusions 



in coordinate space [46j. We liave carried out a comparative analysis of the EoS of PNM 
and SNM obtained using the different parametrizations of the NNLOL potential, as well 
as the improved versions of the TM model discussed in Ref. |45) . 

The calculation of the SNM EoS has been been performed within the variational 
FHNC/SOC approach. In the case of PNM we have also used the AFDMC computa- 
tional scheme, the results of which turn out to be in close agreement with the variational 
FHNC/SOC estimates. 

Our analysis shows that the transformation from momentum to coordinate space 
brings about a cutoff dependence, leading to sizable effects in nuclear matter. As dis- 
cussed in Section 13.3. 1[ the contribution of the contact term, which in PNM would vanish 
in the A oo limit, can not be fully determined fitting the low energy observables. 
Moreover, the NNN contact terms of the NNLOL2 and NNLOL3 models turn out to be 
attractive in PNM, leading to a strong softening of the EoS. 

An illustrative example of the uncertainty associated with the local form of the NNN 
contact term is provided by the results of Fig. 13.161 and Table 13.61 The NNLOL4 model 
largely overestimates the empirical value of the compressibility modulus of SNM, thus 
yielding a stiff EoS. On the other hand, as pointed out in Section I3.4.2[ it predicts a soft 
EoS of PNM. 

The impact of this is ambiguity is large, since compressibility is a most important 
property of the EoS. The recent discovery of a ~ 2 Mq neutron star appears in fact to 
rule out dynamical models yielding a soft EoS of /3-stable matter. 

None of the considered three-nucleon potential models simultaneously explains the 
empirical equilibrium density and binding energy of SNM. However, among the different 
parametrization that we have analyzed, the NNLOL4 and TM3 provide reasonable values 
of po- It has to be emphasized that this is a remarkable result, as, unlike the UIX model, 
these potential do not involve any parameter adjusted to reproduce po- 

In order to resolve the inconsistencies involved in the contact term, one should include 
all contributions to this term arising from the chiral expansion at NNLO. Moreover, as 
pointed out by the authors of Ref. |155] . due to the choice of the regulator function (see 
Eq. fll.36p ). a fully consistent treatment should also take into account NNNNLO contact 
contributions. 

In future works we are planning to include the NNNLO contributions to the three- 
nuclear interactions of Refs. |l0 t [68 ] . as well as the three-body NNNNLO contact terms of 
Ref. |155) . For consistency, at that point also the NNNLO long-range terms calculated in 
[69] should be accounted for. However, since we are not using the chiral potential in the 
two-body sector, we need to deal with the issue of determining the low-energy constants 
entering the three-body interaction. 

The last Chapter of the Thesis reports the results of a calculation of the weak response 
of nuclear matter, carried out using effective operators and effective interactions derived 
within the framework of the CBF approach. This calculation significantly improves on 
those of Ref. \47 \ 1156] . as it explicitly includes the contributions of three-nucleon clusters 
and three-nucleon forces. 

As far as the effective weak operators are concerned, we have shown that the sizable 
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dependence of the response function on tlie clioice of tlie correlation function is in fact 
unphysical. When three-body clusters diagrams are considered, the transition matrix 
elements, in which correlations enter only through the effective weak operators, become 
nearly independent of the correlation functions. 

The three-body cluster contributions described in Section 12.5.31 have been consistently 
included in the construction of the effective interaction. As a first step we have computed 
the EoS of SNM for an hamiltonian including the two-body potential only. In this case, 
the three-body cluster effective interaction provides a EoS of SNM much closer to to the 
one resulting from full FHNC/SOC calculations, compared to the one obtained using 
the two-body cluster effective interaction of Ref. [17| l37] . The leading contributions of 
the UIX potential, emerging at three-body cluster level, have been also included in the 
effective potential. As a result the EoS of SNM exhibits saturation at p 0.18 fm~^. 

Inclusion of the three-nucleon interactions also affects the single particle spectrum, 
leading in turn to a shift of the CHF response as a function of energy transfer. 

The main effect of the three-nucleon force on the TDA response originates from a 
change of the off-diagonal elements of the effective interaction. As a result, the collective 
mode associated with the Fermi transition at |q| = 0.3 MeV turns out to be shifted 
to lower energy, although its magnitude is nearly unaffected by the three-body cluster 
contributions. On the other hand, a depletion of the peak is observed for the Gamow- 
Teller transition. The analysis of the TDA response for different momentum transfer also 
reveals a sizable effect of three-nucleon cluster contributions. 

The sum rules for the Fermi transition comes closer to the variational results once 
the three-body cluster cluster is taken into account, thus confirming the importance of 
many body effects, which are included in the variational calculations through the chain 
summations. The residual discrepancy cannot be accounted for by the n > 4— body 
cluster contributions, and is likely to be ascribable to the effect of multi p — h excitations, 
which are taken into account in variational calculations. This effect appears to be even 
larger even larger in the structure function obtained from the Gamow- Teller response. In 
this case the change due to the three-body cluster is in fact very small. 

Some improvements are observed in the sum rules of the longitudinal and the trans- 
verse response. As shown in Fig. I4.20[ the results of the three-body cluster calculations 
of S^{q) carried out within TDA are slightly closer to the variational ones for all the 
values of |q|, compared to the two-body cluster case. Moreover, the position of the max- 
imum of the TDA calculations of the static transverse response is almost coincident with 
that obtained from variational calculations. At small momentum transfer however, the 
three-body cluster contributions move the TDA S'^{q) away from the variational result. 
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Appendix A 
Hartree-Fock equations 



The starting point to derive the Hartree-Fock equation is the variational principle 

< = (^|//|^) , (A.l) 

valid when the Slater determinant has been normalized to unity. At this aim we require 
that single particle wave functions be orthonormal 

dxi}n,{x)lpnAx) = Sn,n, ■ (A. 2) 



It is straightforward to show |162] that the expectation value of the Hamiltonian is 
made of three contributions 

A 1 ^ 

E[^] = J] 4^ + - ^ [J„^„^, - . (A.3) 

ni=l ni,nj=l 

The first contribution comes from the kinetic energy, as no genuine one-body potential 
appears in the nuclear Hamiltonian of Eq. (II. lip 

r 

In^ = / dXitlj*^^{Xi)^ilJnXXi) , (A.4) 

while the direct term J„.„ and the exchange term Kn^n arise from the two-body potential 



Jn^n, = J dXidXjij^Xi)ij*^^{Xj)Vij'^nXXi)'^nj{Xj) 
Kn.nj = / dXidXjtlj*^^{Xi)tlj* {Xj)Vij1pnj{Xi)tpnA^j) ■ (^-5) 



The variational principle is applied by requiring that the functional E[$] being sta- 
tionary upon variations of the A occupied spin-orbitals ipm, with the constraint that the 
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spin-orbitals must remain ortliogonal. To acliieve tliis, a set of Lagrange multipliers 
e„. „ . is introduced in the variational equation 



6Em- J2 ^n.,n, 
ni,nj=l 



(A.6) 



Since the energy is a real quantity, the Lagrange multipliers are elements of an Hermitian 



matrix, e„.„j. — eJi^,„^. 



The matrix of the Lagrange multipliers can therefore be diagonalized performing a 
unitary transformation on the spin-orbitals 



(A.7) 



The new Slater determinant differs by a phase factor from the previous one, = det{U)'^ , 
while the functional E[%p\ is not affected by this unitary transformation. 

To simplify the notation, instead of working with the primed indeces, we assume that 
the diagonalization has been made from the beginning. Therefore Eq. (1A.6I) can be 
written as 



0, 



(A.8) 



being {e„J the eigenvalues of e^^^.. 

Carrying out the variation leads to the canonical set of integro-differential 

Hartree Fock equations 



(A.9) 



The direct and exchange operators, J and respectively, result from the differentiation 
of the direct and of the exchange term defined in Eq. (lA.5p . While the direct operator is 
local 



(A.IO) 



as in order to evaluate their action on a given function at a point in the configurational 
space, we need to know just the value of the function at the same point. On the other 
hand, the exchange operator 



(A.ll) 



is not local, since we need to know the value of the function ■?/'„. in all the configurational 
space. 



Comparing Eq. (1A.9|) and (12. 2p . it is readily seen that the Hartree-Fock potential is 
given by 



Y.[Jn,-Kn,]. (A.12) 



The solutions of the eigenvalue equation enter the definition of the Fock hamiltonian, 
through the direct and exchange operators. For this reason the Hartree-Fock equations of 
Eq. (1A.9I) are often referred to as pseudo-eigenvalue equations and are solved iteratively. 
The starting point is building the Fock hamiltonian with a guess of approximate spin- 
orbitals. Solving the resulting eigenvalue equation provides new spin-orbitals that can be 
used for the Fock hamiltonian. The procedure must be repeated until the solutions are 
close enough to the spin-orbitals used for the construction of the operator. This method 
of solution is often called self- consistent field. 

Using the orthonormality condition (lA.2p . e„. can be obtained by taking the scalar 
product of Eq. (lA.QO with ipn^ 

^n, = In^ + ^[-^n^nj " -^ninj , (A. 13) 



or, in Dirac notation, 



^I'^i) + y^\ninj\vij\ninj - njUi) . (A.14) 
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Appendix B 
Slater determinant 



The groud-state \I/o of A non interacting fermions is described by the Slater determi- 
nant of single-particle states ipmi^i) 



Isolating the first particle from the determinant amounts in 

^0 = 4t E(-l)"^'^-(^lX^n.(^2 ...Xa) 



(B.l) 



(B.2) 



where '^^^^^^ {x2 . . . x a) denotes the Slater determinant of system with A — 1 particles 
and one hole in the state rii. In other words it is the minor of the matrix defined in Eq. 
f lB.ip from which the first row and the riith column have been removed. 

Isolating a second particle is a slightly more involves issue, as the sign depends on the 
ordering of ni and n2. It turns out that if n2 < rii the multiplicative factor, (— 1)"^"'"^, 
is identical to the one of Eq. ( ]B.2p . When n2 > fii, the fact that the column associated 
with the sate ni had been already removed has to be taken into account because the 
column index of \E'[^_^„^(x2 . . . xa) does not correspond to the state index. In this case, 
the multiplicative factor arising when isolating the n2-th column is then (— 1)"^ 



V V / ni>n2 



(B.3) 



ni <n2 

With "^^^ni 712(^3 ■ ■ ■ •^^) denote the Slater determinant for a system of A — 2 
particles lacking of the states ni and n2. It corresponds to the minor of the matrix of 
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Eq. (IB.ll) from which the first two rows and the colums rti and n2 have been removed. A 
more compact expression can be given for the latter result 

^ - -1-1) ^ (-l)"^+"^+U[^n.(xi)^„,(x2)]vl/^^„^,„,(a:3 ...xa). (B.4) 



ni<?i2 



The many-body state '^p-h is obtained replacing the h—th column of \E'o with another 
column with state index p, having energy larger than the Fermi energy. The procedure 
for isolating particles from '^ph follows the same step as the one of \E'o 

"^vh = E (-l)"^'-"^-''^[V^^i(^i)V^^.(^2)]vl^L'^^,.,(x3 ...xa). (B.5) 

^N{N - 1) „f^^ 



The indexes rii have the values 1, . . . , /i, . . . , A, while 

hi = p if i = h, 
hi = rii otherwise. 

The generalization to the three-particles case is straightforward 



(B.6) 



ni<?i2<"3 



ni<?i2<"3 

(B.7) 
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C.l Four-body reducible diagrams 

In this appendix we report the exphcit calculations of the four-body reducible diagrams 
involved in the calculations of the FR three-body cluster contribution of (fi2)- 
The analytic expression of diagram (a) drawn in Fig. 12.191 is 



V 



dir (a) 



EL 

2A 



rfri234^'(ri3)CTri234 ^12^^12^12(^14 - 1)A3 



(C.l) 



Because of the tracing over the spin-isospin variables of particle 4, only the central part 

^2 
34 



of F^^ does not vanish; hence, only the scalar part of the exchange operator survives 



V 



dir (a) 
4b^3b 



p3 1 

/ (iri2(ir34CTri234 



/ ciri234^2^ri3)CTri234 [Fl2t^l2Fl2(F34-l) 

c/ri2rfri3«'(ri3)CTri234 [^12^^12^2(^24 - 1) 



Fi2f 12-^12 (-^; 



34 



where we have used 



rfri2£2(ri2) 



P 



Renaming the integration variables yields to the expression entering Eq. fl2.16ip 

o2 



V 



dir (a) 
46-^36 



PL 
2 



(iri2(iri3CTr 



123 



Fi2Vi2Fu{F^3 - 1) 



(C.2) 
(C.3) 

(C.4) 



The calculation of diagram (b) of Fig. 12.191 follows analogously. 

The contribution of the diagram depicted in Fig. 12.201 is given by 



nbi3b = ^ I rfri234^(n2)^(ri3)^(r23)CTri234 ^12^^12^12(^34- 1)^12^13 



(C.5) 
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As before, since only the central part of F^^ survives, it turns out that 



(C.6) 



Thanks to the property 



drui{n^)£{\Tu-ru\) = -e{n2) 

P 



(C.7) 



and relabeling the integration variables, one arrives to the following expression 



<^3b = / dr,2drUi{ruyCTTi23 F,2Vi2Fu{F^^ - l)Pi2 



(C. 



Notice that, to recover the contribution to Eq. fl2.165p . one the symmetry 1 f-)- 2 needs 
to be exploited. 

Let us now consider the diagram of Fig. 12.211 



pl3 
^46-5- 3fe 



1^ / rfri234^(ri3)£(ri4)£(r34)CTri234 FuVi2Fi2{F^, - l)PisP3^ 



(C.9) 



Carrying out the trace over spin-isospin variables of particle 4 makes only the central part 
of (^3^4 — 1)P34 not to vanish. Hence only the scalar part of P13 contributes 



<^36=|^^ / C?ri234^(ri3)£(rM)^(r34)CTri234 F,2Vi2Fl2{F,\-l)P3, 

= / t/ri2rfrl4t/r34^(ri3)£(ri4)£(r34)CTri234 ^12^12^12(^14 " 1)^34 

= ^ / rfri2rfr34^(r34)'CTri234 Fi2Vl2Fi2{F'i^ - 1)^34 



(C.IO) 



where in the last line we have used Eq. (lC.7p . A change of integration variables leads to 



pl3 



Mb^sb - y / rfri2rfri3^(ri3)'CTri23 ^12^^12^12(^13 " l)^i3 



(C.ll) 



The diagram of Fig. 12.221 can be written as 



' dri234^(ri2)^(ri3)^(r24)^(r34)CTri234 [^12^^12^12(^34 " 1) A2A4 



P_ 

/ 

A V 



y- / rfri2rfri3rfr34£(ri2)£(ri3)^(r24)^(r34)CTri234 ^12^^12^12(^14 - 1)^2^34 



(C.12) 
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having used the fact that only the central term of (^3^4 — 1)-P34 contributes. Since i{r24) 
^(|r34 — ri2 + rial), exploiting the convolution property of Eq. (lC.7p . it turns out that 



4U3b = -P" / rfri2rfr34^(ri2)£(|r34 - ri2|)^(r34)CTri234 F^2Vl2Fl2{Fl - 1)^12^34 



V / c?ri2rfr34^(ri2)^(|r34-ri2|)£(r34)CTri234 ^12^^12^12(^24 - 1)^4^, 



12 



p2 / rfri2t/ri3£(ri2)£(ri3)£(r23)CTri23 F^2Vi2Fi2{Fl^ - l)PnPi2 ■ (C.13) 



Appendix D 

Green's function for importance 
sampling 



In this appendix we report the derivation, originally due to Kalos, of the expression for 
the importance sampling wave function. For the sake of simplicity we will limit ourselves 
to an accuracy of Ar^, that can be achieved by using the Trotter Fromula of Eq. (12. 206 p . 
Hence, the "original" Green's function is 



G{R, R',At) 



m 



3A ,2 



Let us start from Eq. fl2.216p . that it is worth rewriting 

MR) 



G{R,R',Ar) = G{R,R',At) 



MR') 



(D.l) 



(D.2) 



We want to include the ratio of the trial wave functions in the diffusive part of the Green's 
function. Assuming that for small Ar the displacement \R — R'\ be small, a first order 
Taylor series expansion in i? — i?' yields 



GdiR,R\AT) ~ A/'exp 
~ A/'exp 
~ AAexp 



m 



m 

2ffAT 

m 

2h^AT 



{R-R 



/\2 



1 I ^^MR') 



{R-R'f + vd{R!) ■ {R-R') 



m 



where 



Af 



m 



M 

2 



(D.3) 



(D.4) 



.2TTh^AT. 

As usual, the walkers distribution is represented by a sum of discrete delta functions, 
centered in the position of each walker 



^l;{R,To) = J2^iR-Rk) 

k 



(D.5) 
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Green's function for importance sampling 



In general, the imaginary time evolution does not preserve the normalization. The 
branching factor is nothing but the multiplicity of the walkers that at r + At generated 
starting from R'. Hence an integration over the possible final positions has to be performed 

MR) 



N{R') 



dRG{R,R',T) 



MR') 



dRG{R,R',T) = 

To compute the integral, as before we perform a Taylor series expansion 

N(R') ~ f dRGiR! R, r) [l + ^^^M^ " (R ' R') 
J L Vt{R) 



(D.6) 



+ 



(D.7) 

i,j=l a,/3=l 

where denotes the a-th cartesian component of the coordinates of the i-th particle. The 
Green's function depends on the integration variable only through the gaussian diffusive 
term, which is symmetric in the relative coordinate R — R'. Therefore, only even terms 
survive 

A 3 



N{R') 



1 



2MR 



i=l a=l 

X exp[-AT{V{R') - Et)]. 
Using the following result for the Gaussian integral 

m{R - R! 



m{R - R 



i/^2 



Af J dRexp 
for the multiplicity we obtain 

AtV^MR' 



I \2 



At. 



(D.8) 
(D.9) 



N{R1) 



1 + 
exp 

exp 
exp 



2 MR') 
1 v^MR' 



At 



2 MR' 

^fHMRn_^v 

\ MR') 

Ar(Ei(i?')-^T) 



X exp[-AT{V{R') - Et)] 
[ + Et- V{R!) 



Therefore, the branching factor reads 



Gb{R,R',AT) = e 
hence the full Green's function turns out to be 



Ar El(R')-Et 



G{R, R',At) 



m 



3A 
2 — - 

e ' 



R~R' 



.2Tlh?AT 

which differs from the one of Eq. (I2.217P by terms of the order Ar^ 



(D.IO) 

(D.ll) 
(D.12) 



Appendix E 



Matrix elements of Fermi and 
Gamow- Teller transitions operator 



In this appendix we give the results of the Fermi and Gamow- Teller matrix elements 
of the diagrams computed in section 14.11 As far as the two-body clusters are concerned, 
our findings coincide with those of Refs. [157J, while the three-body contributions have 
not been computed so far. 

To simplify the notation the subscripts i for the particle and j for the hole states 
have been omitted. Hence, the spin and isospin states of the particle and of the hole are 
denoted by Sh, th, Sp and tp. 

Note that the Kronecker deltas St^^p and Stf^^n are common to all the matrix elements, 
showing that in the transition a neutron decays into a proton. The spin-structure, on 
the other hand, is in general more involved; particularly for the Gamow- Teller transition 
where we need to distinguish the three cartesian components. 



E.l Zeroth order 



{ap\Oar{l)\ah) 



(E.l) 



Fermi 



gv{(yp\T^ \ah) = gvSsp,sJtp,pSti,,n 



(E.2) 



Gamow- Teller 



gAiaplA '(Ji\ah) 
I I (+) yi \ 

5'A(ap|ri^Vi|a/i) 



(E.3) 
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Matrix elements of Fermi and Gamow- Teller transitions operator 



As a general property of the Gamow- Teller transition, the matrix element between the spin 
states Sp —l,Sh — i is equal to minus the matrix element between Sp Sh =t- Moreover 
the matrix element between Sp —-IjSh —t is the complex conjugate of the one between 
Sp =^,Sh = i- For brevity we report the results for Sp Sh —t and Sp — "[,Sh— i, while 
the ". . . " symbolize the other contributions. 



E.2 First order 



O 



(1) 

Na 



^(aiap|{/i2 - l,C'aT(l)}|aitt/i) 



ai 



Fermi 



(E.4) 



9v^{aiSp\{fi2 - l,Ti'^^}\aiSh) = 9vSt^,pStf^,nSsj„sh^fl2 



(E.5) 



Gamow- Teller 



2^12^12 



ai 

,t[24/: 

L \ ' 12 ^ 

...} 



12 



(E.6) 



gA^{aiSp\{fi2 - l,r^^Vf}|Q;iSft) = gASt^,pStf„n[ 



ai 



[24/, 



^12^12 

12 „2 
'12 



+ 



...} 



3|/i2(a:i2 -iyi2) 



)] 



+ 



(E.7) 



E.2 First order 
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o 



(1) 

Nb 



ai 



...}. 



{Xi2 - iyi2)zi2 



J.1 

'12 



+ 



)] 



+ 



^(aiap|{/i2 - l,6(,r(l)}AYl"i«ft) 



Fermi 

J](aiSp|{/i2 - l,r(+)}A7l«i5;.) = 

Oil 

gvK,pK,n5s,,snmt2 - 1) + fl2 + 3(/r2 + AT)] 

Gamow- Teller 



ai 



3^12^12 



[6(/i*2 - /S) 

(/r2 - 1) + fi2 + /r2 + + (/i*2 - ( - 1 + 3 



(E.8) 



(E.9) 



(E.IO) 



2^12 (a;i2 - %2)~ 



...} 



12 



)]- 

(E.ll) 



gA^{aiSp\{h2 - l,T['^^al]P^i\aiSh) = gA5tp,pK,n{^ 



ai 



tr 

12^ „2 
'12 



+ 



- i{n^ - 1) - ^/[^ - ^A^2 - ^^fl2 - {fi2 - fll) - ^ - 3 



,yi2(a;i2 -«?/i2)y 



r2 

'12 



(E.12) 
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Matrix elements of Fermi and Gamow- Teller transitions operator 



Ss,,Ss,,2 [(/r, - 1) + fi, + A + 5/- + ifi, - fii) (i - + - "'^^ 
<^..A,^[6(/i*2-/i*2: 

...} 



^2 
'12 



^i2(a;i2 - iyu) 



j,2 
'12 



+ 



o 



(1) 

Nc 



Fermi 



OC2 



Gamow- Teller 



9A^{Spa2\{fl2 - 1,t}^Vi}|s,jQ;2) = 9AStp,pSt^,n\^ 

'^..A,48[A-i]- 

...} 



+ 



+ 



(E.13) 



(E.14) 



9v ^{Spa2\{fi2 - 1, Ti^^}\sha2) = 9vK,pK,n5s^,snKf 12 - 1) (E.15) 



(E.16) 



9A^{Spa2\{fl2 - l,T[^^al]\Sha2) = 9AStj„pSt^,n[ 



02 



+ 



...} 



(E.17) 



9A^{spa2\{fi2 - l,ri^^Vi}|sfta2) = 9ASt^,pK,n[ 

...} 



+ 



(E.18) 



E.2 First order 
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The spin-isospin matrix elements of are identical to those of 0)^^ for both Fermi 
and Gamow-Teller transitions. 



o 



(1) 

Da 



^(Q;iap|/i2 - l\aiap) = 4(/i2 - 1) 



(E.19) 



o 



(1) 

Db 



5^(«l«p|(/l2 - l)A7l«l«p) = ifl2 - 1) + HfL + fl2 + 3/lT) 



(E.20) 



ai 



We do not report the results for the spin-isospin matrix elements of O^^ and O^^^ as 
they are identical to those of O^a ^^nd O^^, respectively. 



Fermi 



Gamow-Teller 

^A^(QilQi2Sp|T^'^Vi(/23 - l)Pi2\0iia2Sh) = gAStp,pSt^,n\^ 
<^Sp,t<^5h,t [12/23 1^2 



+ 



23 



+ 



...} 

<^5p,t^^/.,tfl2/23^^5^ + 
' 23 

A A r ^„-/<xr , ,ftrf, , 3|/23(a^23 -^I/23) M , 

...} 



(E.21) 



9V ^{aia2Sp\T['^\f23 - l)A'7|ai«2Sfe) = 9vK,pK,nK,^h'^f23 (E-22) 



(E.23) 



(E.24) 



168 



Matrix elements of Fermi and Gamow- Teller transitions operator 



o 



(1) 

Nl 



L \ ' 23 ^ 



[12/, 



2:23 (a;23 - W23) 



...} 



23 



+ 



(E.25) 



(E.26) 



Fermi 

gV^{aia2Sp\ri'^\f23 - l)Pi7A7l«l«2S^) = 9vStp,p5tf^,nSsp,SH4=f23 

Gamow- Teller 

2^23^23 



(E.27) 



<^..A,t[6(/23 -/23) 
K^^s,,, \{r23 - 1) - /2''3 - /23 + 123 + 2(/^3 ' Hi) ( " 1 + 3 ^Hii^^^^li^ ) 1 + 

. . . } (E.28) 



9A^{aia2Sp\Ti'^^a^{f23 - l)-PiTA7l«l«2Sfe) = yA^t^.p^t^.nj 
^..,t^s„t[6(/2*3-/23) 

<^.p,t5..,4 [ - ^((/23 - 1) - /23 - f23 + AY) + 2(/2*3 ' ftl) ( ^ + 3 

...} 



, 2/23^23 



23 



^23(3:^23 - ^y23) y 



23 



(E.29) 



E.3 Second order 
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(/23 - 1) - /Ja - /23 + /23^ + 2(/23 " /23) ( 2 " 3 3 ) 

^ '23 ^ - 

'23 - 

...} 



+ 



(E.30) 



E.3 Second order 

In order to write in a compact form tlie Gamow- Teller matrix elements, it is worth 
introducing the following quantities 

M =24[/^2(/i*2 + 2/i*D - A7(/i*2 + + A*2((/r2 - 1) - K2 + 2(/L + /i*2)] 

^2 =8[/^2/f2 + (/l2 - + "^fufll + '^fufll + 2/lT^ - /r2/l2 + AT/L" 

(A2 - 1) A*2 - 2A"2A*2 + A'^2A*2 + 2ATA*2 - 4A*2A*2 - 4A*2'] 

Bl ^Q[fUfu - 1 + 3A"2 + A"2 - 5AT + 2A*2) + (-(A2 - 1) + A"2 + 7A'^2 + fll+ 

B2 =(A2 - 1)' + 2(A2 - 1) A"2 + A"2' + 2(A2 - 1) A'^2 - 6A"2A'^2 + A'^2' + io(A2 - 1) AT+ 

2/12/^2^ + "^fuf 12 + fl2 ~ 2(/{^2 ~ l)/l2 ~ 6/^2/12 ~ 2/^2/12 + 10/i7/i2 — 

8A*2' + 2(A2 - 1) A*2 - 2A2A*2 - i4A^2A*2 - 2A7A*2 - i6A2A*2 - 8A*2' 

C^ =24[2/f2A*2 + A*2' - A*2(2A7 + A*2)] 

C2 =4[(A2 - 1)' - A"2' - A"2' + AT' - 4A"2A*2 - 4A*2' + 4ATA*2 + 4A2'] • (E.ai) 



J](aiap|(/i2 - l)a.(l)(/i2 - l)\aiah) (E.32) 

ai 

Fermi 

ai 

gvK,vK,n^s,,sAfl2{{ft2 - 1) + A"2) + 3AT(A"2 + A'^2^) + 6A2(A*2 + A2)] (E.33) 

Gamow- Teller 
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Matrix elements of Fermi and Gamow- Teller transitions operator 



9A^{aiSp\{fi2 - 1)tJ^Vj^(/i2 - 1)|q;is^) = gASt^,pdt^,n[ 



12 



+ 



...} 



'12 



+ 



(E.34) 



ai 



[- 



?/i22;i2 



^?2 



+ 



^s,,^^suA - ^^2 + >ll 



?/i2(a;i2 - iyx'i) 



^2 
'12 



+ 



(E.35) 



o 



(2) 
Nb 



Fermi 



^12 
12-" 



+ 



...}. 



12 



+ 



^(ai«p|(/i2 - l)0..(l)(/i2 - l)P[2^\a,ah) 



gvY.{a,SpXh2 - 1)t^H/i2 - l)Pf2l«i%) 



ai 



2(/r2 - i)(/r2 + 3(/r2 + at)) + i2(/i*2 + f%?\ 



(E.36) 



(E.37) 



(E.38) 



E.3 Second order 
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Gamow- Teller 

3^12^12' 



oil 



12 



+ 



...} 



12 



+ 



(E.39) 



ai 



[- 



yi2Zl2 



12 



^12(2^12 - iyu) 



...} 



12 



+ 



(E.40) 



o 



(2) 
Nc 



^^^^("l^PiKAs - l)-rJ^^0-^(/l2 - 1) At I ) = 9AStp,pStf^,n[ 



^Sp,^Ss^,^ \B2 + Bi-^ 
'- '12-' 



+ 



...}. 



^12(2:12 - ^2/12) 



^2 
'12 



+ 



X](ttpa2|(/i2 - l)C'<,^(l)(/i2 - l)|afta2) 



Fermi 

a2 



(E.41) 



(E.42) 



(E.43) 
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Matrix elements of Fermi and Gamow- Teller transitions operator 



Gamow- Teller 

2^12^12' 



02 



12 



...} 



+ 



2^12(2^12 - iy 12) 



12 



+ 



9A^{Spa2\{fl2 - l)'r}^''o-l(/l2 - i)\Sha2) = gASt^,p5t^,n\^ 



0:2 



1/12^12 



12 



+ 



|/i2(a;i2 - %2) 



12 



+ 



...} 



gA^{Spa2\{fl2 - 1)tJ^V[(/i2 - l)|s/iQi2) = ^A(^tp,p(^t;„n{ 



a2 



<^sp,t<^sa,t 1^2 + Ci 



-12 

rl2^ 



+ 



Ci 



2:12 - ^1/12 J 

J.2 
'12 



+ 



Like in the first order case, the spin-isospin matrix elements of are identical to 



(E.44) 



(E.45) 



(E.46) 



(2) 

those of 0]y^ for both Fermi and Gamow- Teller transitions. 



o 



(2) 
Da 



- l|«i«P> = 4[(/r2 - if + 3{ru + fu + 3/iT' + 2/1*2' + Qfif)] (E.47) 



o 



(2) 
Db 



^(aiap|(/i2 - l)Pf2'l"i"p) 



ai 



{{fi2 - 1) ' + 6(/r2 - i)(/r2 + fi2 + 3/i7) - 3[/r2' - 6/r2/r2 + /r2'+ 



QKIifL + K2) - 3/iT' - ^{fl2 + 6/1*2/1*2 - 3A*r)]} 



1 2 



(E.48) 



E.3 Second order 
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We do not report the results for the spin-isospin matrix elements of O^^ and 0[^^ as 
they are identical to those of O^]^ and O^^, respectively. 



Appendix F 

Total spin of a particle-hole pair 



We have shown that the Hartree-Fock ground state of infinite nuclear matter is a Slater 
determinant of plane waves. The single particle wave function can be characterized by the 
quantum number \a) = |k, s, s^, t, t^). with s = 1/2 and t = 1/2. Since our Hilbert space 
is a direct product of spin and isospin spaces, for the sake of simplicity we may consider 
only the spin, as the treatment of the isospin degree of freedom proceeds analogously. 

The general creation and distraction operators are denoted by c|^^ and Ck,s^, where 
it is understood that s = 1/2. The particle and hole operators can be defined by the 
relations 

bls^=S-^,-s^C_^,-s. \k\<kF. (F.l) 

The following phase convention has been adopted 

^k,.. = (-1)^/2-^^ . (F.2) 
This is a canonical transformation, for it does not affect the anticommutation rules 

{ak,s,, Ok/^^j = {bk,s,,bl'^s'J = Skk'Ss.s'z , (F.3) 

while all the the other anticommutators vanish. It can be shown |i52j that the s^-dependent 
phase operator foj^^^ creates a hole of spin projection along the z-axis Sz- 

Using the canonical transformation we may now compute the total spin of the particle 
hole pair \pm', hi). We start from following the well known relation for the total spin of 



the particle particle states 

4/-h,tlO) =^ ^ = l,^, = l (F.4) 

;7f«/-h,, + </-h,t)|o) ^ ^ = i,^. = o (F.5) 

<A,JO) =^ S = l,Sz = -l (F.6) 

4«/-h,4-</-h,t)|0) =^ S = 0,Sz = 0. (F.7) 
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Total spin of a particle-hole pair 



Through the canonical transformation we can relate these states with the correspond- 
ing particle-hole states 

(F.8) 

cl,^c^,m (F.9) 

(F.IO) 

4,;Ch,4)|0). (F.ll) 

By comparing Eq. (IF.7P and (IF. lip , we immediately find the results listed in Table 



t /,t 



^-h,t)|0) 
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Spostandomi di 700 Km a Nord-Est, ringrazio gli amici del tennis, tra cui Giacomo, 
Andrea, Piero e Gabriele: sono stato accolto come uno di loro. Gh allenamenti, le partite 
e la palestra in loro compagnia sono stati fondamentali, soprattutto quando ero nero per 
dei conti che non tornavano, o per delle soluzioni che non si trovavano. 

Nella citta adottiva ho avuto il privilegio di conoscere Daniel, un amico sincero, con 
cui ho condiviso momenti "epici ed eleganti" durante questi quattro anni. Se mai saro 



nonno, egli sara presente nei racconti che faro ai miei nipoti per via dei sei giorni in 
bici da Trieste a Roma, i quattordici da Milano a Santiago di Compostela (con annessi 
i cinque giorni aggiuntivi nella sua natia Galizia, terra di Breogan), e delle varie gite 
verso Tarvisio, lo Zoncolan, Venezia e passo Pramollo. Lo ammiro non tanto per essere 
un ciclista decisamente migliore di me, come dimostrano i minuti di ritardo che accuso 
appena la strada inizia a salire, ma per le sue qualita di uomo, veramente di rara eleganza 
(esula da questo discorso la scelta delle scarpe); Epic and Glory! 

Per la sua laconica e saggia ironia, che da giusta misura a difficolta a prima vista 
enormi, ma soprattutto perche "altrimenti si arrabbia", sono riconoscente ad Andrea, 
altresi noto come Jens. Come non nominare anche I'altro compagno di cadute, Castillo 
alias Edu, non solo perche mi abbia citato nei ringraziamenti della sua Tesi, ma perche 
da sempre il massimo, anche quando non strettamente necessario (e per questo detiene il 
record di scalata a piedi dello Zoncolan). 

Dalla Galizia alle Asturie, per aver costantemente allenato la mia eloquenza, perche 
«La cultura o ce I'hai o non ce l'hai», trova giustamente posto in questi ringraziamenti 
Irene. Non e da tutti perseverare nell'errore di frequentarci e di fingere di ridere alle nostre 
battute, cara Isabelita. A proposito "S'z Evita estuviera viva, Isabel seria soltera... " sono 
certamente in debito con Ignazio per questa canzone e per il suo umorismo sorprendente. 

Bene, mi sto dilungando forse troppo, ma dopo 190 e passa pagine, aggiungerne due 
o tre non e che cambi poi molto la situazione. Certo, avrei potuto cavarmela con un " mi 
trovo nella fortunata situazione di avere tante persone da ringraziare e poco spazio per 
farlo |20] "; tuttavia non posso prendere spunto dalla Tesi di Akmal pure per i ringrazia- 
menti! Vado avanti. 

Nicola, per I'impegno profuso come rappresentante, per le discussioni alia lavagna e 
per I'ottimo manzo all'olio. Giorgio e Maurizio, per aver condiviso i primi tempi sotto 
lo stesso tetto, per la calma e per la stimolante anfetaminicita, rispettivamente. Matteo, 
Kansas city e di per se un motivo bastante per essere qui presente. Poi Ambra, Daniele, 
Alberto, Aurora (anche per il polpo tramite Michele!), Umberto... 

Ancora Daniele, un altro pero, precisamente il mio attuale coinquilino, nonche com- 
pagno di giuochi universitari. Grazie per le cene, per i dialoghi notturni, per non aver 
mai fatto caso alle spigolosita del mio carattere. Se verrai nella mia nuova casa una bella 
pulizia del piano cottura non te la leva nessuno! Dovrai portare anche Giorgia (che puo 
venire anche da sola, eh), solerte compagna di colazioni mattutine, nonche abile autrice 
della foto sul visto, ottenuta con espedienti fotografici d'avanguardia. Non ti invito sola- 
mente per migliorare la danielesca qualita della pulizia del suddetto piano cottura, non ti 
preoccupare! 

Da Trieste a Trento. Colgo I'occasione per rendere edotti gli sparuti lettori del fatto 
che queste due citta, che nell'immaginario collettivo sono quasi sovrapposte, forse perche 
italiane entrambe dal '18, sono assai distanti; metaforicamente, come il mare e la mon- 
tagna. Paolo "Ormoni", senza di te I'AFDMC avrebbe molti piu segreti di quanti ancora 
non ne abbia, cosi come il monte Bondone e il Trentino tutto. Ammiro Diego in partico- 
lare per I'ostinazione e I'intuizione; lo ringrazio per le cose fatte e per quelle che fara per 
risolvere alcuni "problemucci" del codice, per non parlare della minuziosa rilettura della 



presente opera. Elia, per le sciate e per aver ricambiato la visita. Enrico, senza il tuo 
aiuto (e quelle di AURORA che non e una persona) dovrei ancora rispondere al referee. 
Giorgia, per essere la barra di grafite nel reattorc della saletta. Tra I'altro vorrei sotto- 
lineare I'economicita della pensione "Roberto" a Pergine: un buono-cafFe a notte, "extra" 
inclusi (e che "extra"). 

Giusto per essere il Trip Advisor de noantri, consiglio anche la pensione "Arianna" a 
Barcellona. Squisita ospitalita al costo di un piatto di polpette (vanno bene anche sciape). 

Tempo di serieta. 

Ringrazio Stefano per avermi iniziato all'AFDMC, per I'infinita pazienza dimostratami 
e per avermi procacciato un posto di lavoro. Non cambiando famiglia, sono in debito con 
Serena e le birbette, per I'ospitalita in quel di Los Alamos e per il supporto morale 
seguente. Vorrei calcare i vostri passi. 

Senza Francesco Pederiva e Kevin Schmidt il Capitolo suU'AFDMC non avrebbe mai 
visto la luce, tra I'altro perche I'AFDMC stesso non sarebbe mai nato. Resto sempre 
stupito dalla chiarezza, dalla competenza e dalla pazienza con cui cercate di far entrare 
nel mio cervello quanto i vostri abbiano prodotto. 

E' d'uopo ringraziare Stefano Fantoni, le cui idee mi limito a sviluppare; spero che 
I'operaio possa imparare dall'architetto. Mi stupisco ancora come abbia potuto, a ragione, 
crcdcrc che Ic cquazioni di FHNC, in apparenza molto intricate e soggette a millemila 
possibili errori, potessero portare a risultati cosi accurati. Credere nelle idee a lungo 
termine e un insegnamento prezioso, anche al di fuori della Fisica. 

Omar, dopo quella della mia laurea specialistica, ti ripromisi che mai piu avresti 
corretto una Tesi ad Agosto. Mi sento un poco responsabile di averti fatto rimangiare 
tal solenne promessa. Ogni volta che discutiamo di Fisica, o di politica (notare I'uso 
improprio delle maiuscole), c'e sempre qualcosa che dovrei sapere, di "ovvio" che non 
so e che tu invece sai. Ti ringrazio per lasciar trasparire quello che pensi al riguardo, 
piu dal tuo viso che dalla tua voce. Mi sto anche accorgendo che alcune tue espressioni 
idiomatiche stanno entrando nel mio parlato. Grazie di cuore. 

E veniamo alia Veronese piu terrona della storia; braccio quando c'e bisogno di una 
mente e mente quando c'e bisogno di un braccio: Cristina. L'ultimo Capitolo e largamente 
frutto delle discussioni avute a Natale, a Pasqua e a Ferragosto, piu in quasi tutti gli 
altri giorni feriali e festivi nel mezzo. Magari avrei da dissentire sui tuoi gusti culinari, 
forse guidati dal figlio che porti in grembo, a cui sono gia affezionato, ma sei veramente 
insostituibile, non solo come fisica. Ecco, te I'ho detto, mo' accontentati per i prossimi 
cinque anni. 

Dulcis in fundo, Barbara. Qui mi limito ad un sintetico, ma maiuscolo, "GRAZIE", 
riservandomi il diritto di approfondire in private. 



